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Means and Examples of Aeronautical 


Research in France at ONERA 


The Twenty-Second Wright Brothers Lecture 


MAURICE ROY* 
Office National d' Etudes et de Recherches déronautiques 


OPENING REMARKS 


On this anniversary may I suggest that we turn our minds back 
fifty years to the end of 1908. 

The exploit at Kitty Hawk was already 5 years old. An 
unfortunate accident has occurred on September 17, at Fort 
Myer, in which Orville Wright had been seriously injured. 

Meanwhile, his brother Wilbur was in France, in the region of 
Le Mans. Since August 4, he had been making repeated flights 
of ever-increasing duration which had aroused immense enthu- 
siasm and were to bring world-wide fame to the Wright brothers. 

Fig. 1 shows the famous biplane flying in the twilight above the 
military camp at Auvours before landing after nightfall. In 
France on that day, October 3, 1908, Wilbur Wright made the 
first night flight in the history of aeronautics. 

Fig. 2 dates back to October 15, 1908. On that day, flying with 
his passenger at an average height of 33 ft. above a racetrack and 
its surrounding meadows, Wright remained airborne for 1 hr. 9 
min. 45 sec. Officially, he had covered a distance of 34 miles, 
but, in fact, the figure lay between 45 and 50. In so doing, 
Wright simultaneously smashed the world distance and duration 
records for flights made with or without a passenger, and this feat 
was to secure for him a contract worth $100,000. 

For me today, this particular picture, which recalls that 
triumphant exploit, has a very special significance, for Wright’s 
passenger on that flight was none other than Paul Painlevé, 
who was both a Member of the French Academy of Sciences and 
Professor of Mechanics at our ‘‘Ecole Polytechnique’’—two titles 
which, it so happens, are also held by the Frenchman whom the 
Institute of the Aeronautical Sciences has called upon this year 
to deliver the Wright Brothers Lecture. 

Perhaps then, when I tell you here how grateful I am for this 
great honor conferred upon me, which I deeply appreciate, you 
will bear with me if I let my thoughts go back fifty years to 
Wilbur Wright's flights in my country and to the man who was 
with him when he established that record—the man whose pupil 
I was at the Ecole Polytechnique 30 years before I returned there 
to occupy the same Chair. : 

Before starting my lecture I would like to give well-merited 
tribute and thanks to all of my collaborators at ONERA who 
helped me in preparing this paper. 


Presented before the Institute of the Aeronautical Sciences in 
the Natural History Building Auditorium, Smithsonian Institu- 
tion, Washington, D.C., December 17, 1958. 

* Director. Also, Member of the Academy of Science. 
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INTRODUCTION 


OSMONAUTICS is currently very much to the fore- 


front in the news. 


It embraces and extends aero- 
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nautics, and I would like to propose including both, at 
least on certain occasions, under a general denomina- 
tion of ‘“‘aerocosmonautics.”’ 

In your country, the sciences and technology of 
space are subjects which have been backed by initial 
advances and abundantly treated. 

Since France has not yet launched any artificial 
satellite or built any circumlunar space vehicle, I pro- 
pose to confine myself here to the field of aeronautics, 
where there is still so much progress of manifest utility 
to accomplish. 

I shall accordingly content myself with presenting 
some examples of aeronautical research and expeti- 
ments undertaken in my country by ONERA, a body 
whose mission is akin to that of the illustrious NACA, 
now NASA, but bearing in mind the considerable dif- 
ference between the scales of the respective resources. 


MATERIALS 


The progress in materials always influences construc- 
tion and sometimes even determines design. The his- 
tory of aeronautics, as will that of ‘‘cosmonautics’’ in 
the future, bears striking testimony to this. 

First, the characteristics of metallic materials have 
evolved to a great extent; their behavior has assumed 
ever-increasing importance under conditions of high 
temperature—initially in piston engines, then in jet en- 
gines and rocket motors, and, finally, in every form of 
structure subjected to aerodynamic heating due to very 
high speeds. 

Second, materials of a very different nature, rather 
improperly called plastics, have made their appearance. 
Obtained by a solidifying and hardening polymeriza- 
tion process, these materials were at first usable only at 
low temperatures, but in the future they too will be sub- 
jected to more and more stringent requirements under 
heat conditions. 

In the wide field of metal and plastic research, the 
work undertaken by ONERA has been directed mainly 
at acquiring a deeper knowledge of the relationship 
which exists between the structure of such materials 
and their behavior when subjected to mechanical and 


thermal stresses. In all this work, emphasis has been 
laid on the time-history of the phenomena, that is to 
say, on the influence of the time factor and, hence, of 
velocities and accelerations. 

This objective, in itself highly scientific, is a very dif- 
ficult one to attain. It falls sufficiently outside the 
usual targets sought by industry to preclude any over- 
lapping in research work, though at the same time 
fostering intelligent and understanding cooperation be- 
tween research and industry. 

I propose, first of all, to cite two means of investiga- 
tion developed on these lines. 

The first example relates to the ONERA electronic 
probe (see references 2, 3, 4,5, and 6 and Fig. 3) which, 
by spectrography of the X-rays emitted and diffracted 
by a test sample when struck by a beam of elec- 
trons, enables the local composition of the metal to be 
analyzed quantitatively and its mode and state of 
crystallization to be known. 

The quantitative analysis provides a very satisfac- 
tory degree of accuracy (threshold of detection: 0, 
02°) but takes into account only those elements of 
which the atomic numbers are greater than that of so- 
dium. Like the crystallographic analysis, for any se- 
lected point on the surface of the sample, the quantita- 
tive analysis affects only an extraordinarily minute 
area, which is of the order of two-thirds of a billionth of 
a square inch. 

On the scale indicated for 20 microns, Fig. + shows the 
pin-point areas of a complex sulphide analyzed by the 
electronic probe. 

The second example I wish to cite concerns ma- 
chinery designed to investigate fatigue phenomena in 
materials subjected to very-high-frequency stress re- 
versals; the frequency used initially for these tests was 
6,000 cps, but we are currently working with 90,000 eps. 

This particular machine (Fig. 5), which generates al- 
ternating tensile and compressive stresses, uses sta- 
tionary longitudinal elastic waves formed in an assem- 
bly made up of successive bodies of revolution, these 
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being, respectively, a piezoelectric exciter, a mechani- 
cal amplifier, and the test piece itself. 

Amplitude, the peak value of which is about 10 uy, is 
observed by microscope and enables stress magnitude 
to be deduced in the elastic range. The frequency is 
measured by stroboscope. 

As may be seen, this apparatus is attractive in its 
simplicity, yet allows fatigue tests of an entirely new 
order to be made, it being possible to complete a three- 
billion-stress-reversal test program in only nine hours, 
which is three thousand times faster than with conven- 
tional pulsators. 

Among test data obtained in this way, I would like 
to mention only one fact observed, which is that alter- 
nating elastic strain, at any rate up to the maximum 
frequency used by us, does not alter the ageing process 
on an aluminium alloy of 4 per cent copper content, 
probably because fatigue failure occurred without any 
cold working. On the other hand, the effects of cold 
working can be observed if alternating strain exceeds 
the elastic limit. 

A variety of equipment (see references 9 to 19) has 
been evolved to enable one to follow the evolution in 
the structure of materials while under heat treatment 
during manufacture and when subjected to working 
conditions under thermal and mechanical loads. In 
this context, I would like to mention (a) the process of 
chemical attack known as “‘contrasting,’’ whereby the 
surface tinting which takes place, before the precipi- 
tation themselves are actually observable, (b) the use 
of a reactive agent to condition a surface for the taking 
of alumina replicas which can then be observed under 
an electronic microscope and kept as testimony of suc- 
cessive stages in a given evolution process, and (c) ra- 
dio-crystallography equipment to study, by the X-ray 
method, defects in homogeneity betrayed by satellite 
lines of diffraction patterns. 

Another example of intensive study in the field of 
practical application, using fundamentally scientific 
means, is the surface heat chroming of metals belonging 
to the ferrous group and the plasticization of chro- 
mium. For such work, a halogenous reducing atmos- 
phere is used. 


Fig. 6 illustrates how penetration in depth of the 
chromium-vapor diffusion layer varies according to 
the duration of the treatment at 1,100°C. and the na- 
ture of the steel being processed. 

Fig. 7 shows, with varying and highly significant de- 
grees of magnification, the difference in purity between 
an example of trade-supplied chromium and the new 
metallic chromium obtained by sintering the plasticized 
powder product evolved by ONERA. 

In the following moving pictures,* we shall see first 
the cobalt specimen being put into the furnace where 
controlled alternating heating is observed by a cine- 
microscope. Next small martensitic plates breaking 
surface will appear, then disappear. They show the al- 
most reversible modification under study. 

By means of a microphone and an amplifier, we shall 
hear the crackling of cobalt crystals, that is to say, the 
jerky noise which accompanies the allotropic modifica- 
tion of cobalt heated up to 500°C., a phenomenon 
which was revealed at ONERA and is still little known. 


STRUCTURES AND VIBRATION 


Almost the whole of the activity of ONERA’s 
Structural Division is concentrated on investigating 
vibration phenomena in flight so as to constantly im- 
prove apparatus and methods used to detect these vi- 
brations and to determine the critical speeds of aircraft 


* At this point in the Lecture, a short film sequence was shown. 
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and the remedies to be applied for any shortcoming ex- 
perienced. 

During the six years that ONERA has made it its 
business to ascertain the critical speeds of all proto- 
types produced by French constructors before the first 
flight is made there has been no flying accident in France 
ascribable to structural flutter. In 1957, a dozen air- 
craft and missile prototypes were subjected to this in- 
vestigation which includes ground testing of prototypes. 
For this purpose a road convoy, comprising all necessary 
equipment, visited each manufacturer concerned and 
even went to Italy. On one occasion, when tests had 
to be carried out at our Test Center at Modane, at a 
time when rail and road communications had been cut 
by severe floods, a smaller quantity of equipment was 
transported there by helicopter. 

At present, about ten days are sufficient to put an air- 
craft through the full ground-test program, then some 
twenty days or so are needed for mechanical computa- 
tion and interpretation. 

I propose to draw attention to only three aspects of 
the recent evolution in investigation procedures for 
dealing with such problems, which are still so vitally 
jinked with the safety of aircraft and missiles. 

The first objective is to seek to obtain a precise 
knowledge and assessment of the structural damping 
occurring in an aircraft by in-flight experimentation and 
from tests carried out on the machine itself. 

The first method used,”"?! which is known as the im- 
pulse method and which has already been extensively 
developed, has notably enabled the maximum permissi- 
ble speeds of certain recent aircraft to be raised by 30 
knots and even by as much as 80 knots. In this 
method, powder-energized impulsors are used similar 
to those shown in Fig. 8 mounted on a wing which un- 
derwent in-flight torsional tests; the response signal of 
the structure to the disturbance produced when these 
impulsors are energized, is recorded by telemetry on the 
ground and then suitably filtered prior to being ana- 
lyzed by a novel backwards-reading method. The im- 
pulsors are actually small rockets having a very short 
duration impulse and a thrust usually of between 40 and 
1,000 Ibs. 
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Fig. 9 gives an example of the signal recorded by te- 
lemetry and analyzed with respect to three different 
modes, and indicates the respective frequencies m and 
damping values a. 

The second method, now being compared with the 
first,*? uses in-flight harmonic excitation, as is usually 
done for ground vibration tests. 

Vibrators, similar to those shown in Fig. 10 fitted to a 
Caravelle jet aircraft to investigate the feasibility of 
improving its present performance with a view to its 
future development, are fed by variable-frequency gen- 
erators with power amplifiers and transistors. They 
put the structure into a state of harmonic vibration both 
the frequency and amplitude of which can be regulated 
at will, the amplitude being carefully monitored to ob- 
viate any risk of divergence. 

Recent initial tests using the harmonic excitation 
method have given encouraging results; they have no- 
tably shown that the response signal can be picked not 
from noise by analysis, even when the amplitude of the 
random disturbances is several times greater than that 
of the signal. 

Wind-tunnel tests have revealed also that models 
used to study vibration phenomena constitute excellent 
means for detecting parasitic excitation generated by 
the air stream in the wind tunnel. This form of excita- 
tion, which had generally escaped notice up to the pres- 
ent, manifested itself in our big transonic wind-tunnel 


RESPONSE SIGNAL FIRST MODE FLEXURE n ‘ 
(Hz) 
ll 
| 


—> TIME 
3 NODES FLEXURE 


| 
4 NODES FLEXURE 
“pre 
| 29,3 0,019 | 


rmo 10 kg 


Fic. 10. 


FR 


DA 


at 
tion 
para 
A 
ena | 
for 
mod 
TI 
what 
each 
powe 
be ck 
unity 
pates 
other 
Fis 
tions 
is ass 
ing ti 
In 
vides 
gerou 
essed 
proac 
limite 
tural 
come: 
when 
ered, 
Jus 
chani 
conce 
The 
ous ai 
which 
ated r 


196 
‘ 
© 
0 
; 
C 
nig 
| 
| 
Fic. 9. 
Ha 
| “CARAVELLE’ 7 
Airplane 


30 


AERONAUTICAL RESEARCH IN FRANCE AT ONERA 197 


1001 
| 


| 
= 
| 


FLIGHT 
MACH NUMBER 


G2 
0,2 
TRANSFERS 
42 n 
a a [12] 
— e,, q- S21 21 [4—>2] 
0,4 


at Modane notwithstanding the fact that the test sec- 
tion is located more than 850 ft. from the fans. Such 
parasitic phenomena are worthy of attention by aero- 
dynamicists. 

A second example of research into vibration phenom- 
ena involves the introduction of an ‘‘energy method”’ 
for detecting combinations of fundamental vibration 
modes capable of giving rise to a grave risk of flutter. 

This method of discrimination*’ is based on the some- 
what simplified idea that, in an association of p modes, 
each of these modes receives-——as an algebraic sum—a 
power for dissipation from the other modes which can 
be characterized by the sum of the differences between 
unity and the ratio of the power which the mode dissi- 
pates in binary association with each of the (p — 1) 
other modes, to the power it dissipates alone. 

Fig. 11 shows, plotted against air speed, the devia- 
tions in frequency and damping when an isolated mode 
is associated to another; it also shows the correspond- 
ing transfers. 

In its present approximate form, this method pro- 
vides for qualitative and rapid discrimination of dan- 
gerous associations, the latter being subsequently proc- 
essed by much more precise, detailed computation. 

It is possible—and this is a new and interesting ap- 
proach—to release the energy method from the linearity 
limitations so generally accepted in the matter of struc- 
tural vibration, yet which are so debatable when it 
comes to high—-and even transonic—speeds, or again 
when such components as control surfaces are consid- 
ered, owing to hinge clearance and friction. _ 

Just as it has renovated theoretical and applied me- 
chanics in many different fields, here also the energy 
concept can be beneficial. 

The third example of ONERA research into danger- 
ous aircraft and missile vibration concerns the way in 
which wind-tunnel or free-flight testing is being associ- 
ated more and more closely to the testing carried out on 


the actual aircraft or missiles themselves, both on the 
ground and in the air, as already referred to above. 

Every endeavor is made to obtain reasonable aeroelas- 
tic similarity in models, a thing which it is easier to 
achieve in missiles than in aircraft. An overall check 
of critical speeds is then carried out on these models, 
occasionally self-propelled, either dropped from an air- 
craft** or booster-launched from the ground.” 

Fig. 12 shows, in position on its launcher and 
equipped with a booster, a missile model of this type 
featuring twin simulated ramjets connected to the cen- 
tral body by a straight stub wing. 

Such tests, whether carried out by dropping or 
ground-launching, are also useful for checking and im- 
proving calculation methods. 


APPLIED PHYSICS 


The importance of free-flight testing of models or 
missiles has instigated a great deal of research work in 
the field of applied physics in order first to design, de- 
velop, and improve the wide variety of apparatus needed 
to pick up and measure the interesting quantities, then 
to transmit, receive, process, and correct the experi- 
mental data. 

Two typical setups of this kind designed at ONERA 
are: (1) the radio-trajectograph* which, while the test 
vehicle equipped with a radio transmitter describes its 
trajectory, records the variations of differences in phase 
detected by ground receiver stations, taken two at a 
time, and (2) the 5-channel (Fig. 13) ONERA teleme- 
try unit,” or the 19 channel-type with 4 continuous 
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channels, which is associated with a very full range of 
frequency-modulated pickups; the complete equipment 
weighs very little and takes up little space, while the 
overall degree of accuracy is within | per cent. 

If experiments involving model or missile launching 
and the exploitation of such flights are to be con- 
ducted both economically and efficiently, all primary 
tasks of preparation and execution must be closely co- 
ordinated and all operators must observe a rigorous 
discipline in work. 

These requirements were fully met in a few years at 
ONERA with the formation of special “operational” 
teams, and it is indeed a pleasure for me to be afforded 
this opportunity of praising their competence, con- 
scientiousness, and diligence. 

Fig. 14 shows an ONERA experimental ram-jet de- 
signed for research into propulsion at altitudes of over 
65,000 ft. and speeds in excess of 2,500 m.p.h. set up for 
launching by one of these teams. 

Means used for the automatic handling and reduction 
of wind-tunnel and flight data have been considerably 
developed in France*’ as in the United States, but some- 
what later. 

In the case of our trajectograph, the seven receiver 
stations taken two-by-two yield 21 combinations and 
the plotting of about a thousand points of the trajectory 
is automatically accomplished in an hour and a half. 
An automatic tracking board traces this sequence of 
one thousand poiits with a positional accuracy better 
than | ft. in the three space coordinates. 

Another piece of apparatus worthy of mention® con- 
cerns the automatic reading of photographs of wind- 
tunnel manometer banks (Fig. 15); these photographs 
are introduced into a machine in which a scanning pho- 
tocell counts the number of lines not masked by the 
liquid. By means of electronic devices, these readings 
are recorded on teleprinter tape which then passes 
through the data-reducing machine. 

In half an hour, 50 frames of a 50-tube manometer 
are read, processed, corrected, and recorded. 

I would like to give two other examples of novel phys- 
ical methods that find application in wind-tunnel re- 
search works. 

(1) Interferential strioscopy (Fig. 16), which has 


been improved at ONERA, has very great sensitivity 
and can be used for fairly rapid quantitative evalua- 
tions.*! By means of a film sequence, taken at 4,000 
frames per sec., the little-known pulsations of a bow 
shock wave in the process of formation can be quanti- 
tatively analyzed at each point. 

(2) Magnetic model-suspension,** which 
models to be held in a wind-tunnel test section without 
material means of support (Fig. 17), suppresses all dis- 
turbance and interaction. In this method, a group of 
coils creates magnetic fields, the intensities of which are 
controlled by the position assigned to the model, this 
position being defined in the three space coordinates by 
luminous beams directed at photoelectric receivers. A 
notable contribution of this system, which is particu- 
larly advantageous in research wind tunnels, has been 
to provide fundamental data about base and sting ef- 
fects and about their interactions in the case of axisym- 
metrical bodies. 

Concerning noise and vibrations, I should mention 
that recently some work has been undertaken on jet 
engine noise as well as on the vibrations and the mechani- 
cal and thermal effects induced by the hot exhaust gases 
on aircraft structures. ** 

Moreover, research on the sound-proofing of air- 
craft has been under way for a long time. 

Fig. 18 illustrates the progress made in France in this 
latter respect over the period extending from 1947 to 
1957, thanks to the pioneering work of ONERA which 


enables 


NOISE| civil TRANSPORT 1950s 
~... —BOMBER 1947 
(DB) 
90 


“CARAVELLE” 
JET TRANSPORT__-> 
1957 
150 600 2400 
FREQUENCY (Hz) 
Fic. 18. 


9600 


€ 
= ii, 
2 
v : 
Fic. 17. 
— 
4 
| 
it 


200 


FOURNAL OF THE AERO/SPACE SCIENCES—APRIL. 1959 


stall 
4} 
+--+ 


Fic. 19. 


Rotating.» Regular 


Flow 


06 O8 » 1 
Fic. 20. 


| 


(UPSTREAM) 


~ = Flow coefficient 
= Total pressure 
coefficient 


-e-e-Experimental points 


in regular flow 
eee Experimental points 
for rotating stall 
1+¢, 
c= *=Characteristic 


parameter (constant) 


v=1- $=Relative width 


of the separated flow 


| 
: 114 pz 
x 4 a 122 pz 
130 pz 

20 


Deflection 


10 


Coefficient of momentum of the jet 
| 


| 

| 

TOTAL PRESSURE 
| 


0.05 


Fic. 21. 


010 0.15 0,20 


AIR 
100 Lb/s 


Fic. 22. 


Fic. 23. 


LABORATORY 
HEATER 


QUICK OPENING 
VALVE 


VACUUM 


PRESSURE 


ISOLATING VALVE 


TEST CHAMBER 
Fic. 24. 


EXHAUST GASES tic 
do 
(b) 
| + + 0,2 \ 
i HI 2 in 
0 — roc 
~ $=0,67 
stuc 
bus 
. 
frui 
| 7 4 in w 
5 
Fic. 25. melt 
* 


AERONAUTICAL RESEARCH IN FRANCE AT ONERA 201 


has resulted in the very remarkable acoustic comfort of- 
fered by the French medium-range jet liner Caravelle. 

A last example of scientific research work in general 
physics applied to a very practical case of engineering, 
particularly in the aeronautical industry, concerns the 
ONERA three-dimensional milling machine, all opera- 
tions of which are controlled, according to a concept of 
my own, solely by mathematical equations defining the 
shape of the object to be reproduced.” 

This can be applied with advantage to the reproduc- 
tion of wings or blades, the contours of which are rea- 
sonably easy to define with a few terms from a Fourier 
series, easily handled by simple electronic computers. 

By ways of example, Fig. 19 shows one half of a 
stamping die machined by this process and used, in this 
instance, for mass-producing axial compressor blades. 


ENERGY AND PROPULSION 


ONERA research work on energy and propulsion 
problems has been mainly directed—and deliberately 
so—towards fundamental subjects, and this has been 
done in close liaison with industry, which naturally 
tends to concentrate more on progress in technological 
fields. 

Numerous special and rather interesting measuring 
instruments have been developed, and I would like to 
mention some of these as examples: (a) a gas analyzer, 
based on infrared absorption, which can be used in 
flight** to monitor jet-engine operation and rating; 
(b) a very sensitive and accurate flowmeter,” for use 
in connection with propellants tested in our micro- 
rockets*~—" (as a matter of fact, this flowmeter has 
been adopted to measure blood flow in the “‘artificial 
heart’’ now in service in a Paris hospital); (c) a con- 
tinuous-type electro-optical thermometer*! with a re- 
sponse time of 10~° sec., to measure flame temperatures 
between 1,400 and 3,300°K.; and (d) a spectrograph to 
study the distribution of chemical reactions in a com- 
bustion zone.” 

Thanks in particular to the use of such equipment, 
fruitful research techniques have been recently de- 
veloped for the study of combustion in ramjets and the 
chemical kinetics of liquid propellants. 

For example, film sequences have been made to 
study combustion instability in an ONERA micro- 
rocket equipped with transparent tanks which are pho- 
tographed at the rate of a thousand frames per second 
to provide accurate flow measurements, while pressures 
are recorded to within a tenth of a millisecond. 

In studying aerodynamic heating at high and very 
high speeds—and this has a bearing on propulsion as 
well—rocket jets can be used to provide high-tempera- 
ture supersonic flows, as can be shown in film sequences 
in which an AGARD B ogive is placed in a Mach 2.5 
rocket exhaust jet at an absolute temperature of 
2,600°K. These films* show that the nose of the ogive 
melts almost instantaneously, but that the shouldering 


* At this point in the Lecture, a short film sequence of this kind 
was shown. 


thus formed causes separation of the outer flow and de- 
lays downstream heating. 

This example of a simplified experiment is linked with 
the fundamental problem* of heat exchange phenomena 
occurring in the boundary layer, upstream and down- 
stream of separation—a problem which is now being 
studied and which may give new solutions for the shape 
and structure of re-entry noses on hypersonic aircraft 
and missiles. 

Insofar as gas turbines are concerned, I would like to 
mention the research work devoted to the phenome- 
non of rotating stall which is so detrimental to high- 
efficiency axial compressors and whose underlying the- 
ory has been established and confirmed,” as Fig. 20 
shows in summarized form. 

As one of the possible ways of delaying or opposing 
this highly undesirable phenomenon, I suggested air 
blowing at the trailing edges of guide vanes, the theory 
of which has since been established and verified, as 
shown in Fig. 21. 

This method, which transposes jet-flap research on 
isolated wings (research which I will consider in a mo- 
ment) to airfoil cascades, can find useful application in 

various types of gas turbines, for, as can be seen, mod- 
erate blowing is capable of markedly deflecting the rel- 
ative air flow and can therefore replace a cascade of ad- 
justable stationary blades. 

I will also mention a third example of mixed theoreti- 
cal and experimental research which likewise resulted in 
what I will term an “‘elegant’’ theory, in the sense that 
such theories are not only all-encompassing but also 
highly accurate in their application. This concerns the 
theory of ejectors in which the driving flow is super- 
sonic.” * In distinguishing between supersonic, satu- 
rated, and mixed regimes, this theory has linked up a 
very wide and complex variety of experiments. It has 
been utilized to establish the design of the ejectors 
used at ONERA’s aerothermodynamic laboratory at 
Palaiseau to simulate very high altitudes. 

A similar application (Fig. 22) of the theory, though 
on a larger scale, is projected for an intermittent wind- 
tunnel for ramjets and rockets at Modane. A flow 
at 20 lb. sec. of burnt exhaust gases at a pressure of 
0.1 atm would be extracted by a two-stage air ejector 
and then discharged into the open air. 


LARGE-SIZE WIND TUNNELS 

While on the subject, I would like to point out that 
ONERA has built an Aerothermodynamic Test Cen- 
ter** at Modane, at the foot of the Alps (Fig. 23), which 
makes direct use of a waterfall with a height of 3,000 ft. 

By direct action on water turbines, this waterfall 
drives the 50 ft.-diameter fans of an 85 megawatt tran- 
sonic wind-tunnel, called S1 MA, which, until the Arnold 
Engineering Development Center tunnels at Tullahoma 
came into operation, was the most powerful in the 
world. The waterfall also drives the 65 megawatts 
compressor of the closed circuit S2MA supersonic wind 
tunnel now being completed and to which are associ- 
ated two intermittent tunnels, one of which is the hyper- 
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sonic S3MA while the other, the S5MA, is especially in- 
tended for high supersonic tests on ram-jets and rockets. 
The large 26-ft.-diameter transonic wind tunnel has 
been in operation for 5 years and has been used for nu- 
merous tests already. 

Though direct hydraulic drive is not used in the 
United States, our plant at Modane is well up to U. S. 
dimensional standards, as many of our American friends 
have been able to see for themselves. 

Fig. 24 shows the layout of S3MA which is now ready 
for the first testing. 

Fig. 25 gives an example of one of recent tests. The 
Nord-Aviation ‘‘Griffon’’ has undergone very full high- 
subsonic testing with its turbo-ram-jet engine running, 
and under the control of the same pilot who later put the 
aircraft through its flight testing. 

Fig. 26 shows yet another example: here a Dassault 
‘““Mystére”’ aircraft is set up in the test section in readi- 
ness for cockpit jettisoning inside the tunnel at speeds of 
up to 500/550 m.p.h.; these tests, in fact, proved par- 
ticularly instructive. 

Another type of full-scale test, carried out in the big 
open-circuit wind tunnel operating at Chalais- Meudon 
since 1930, is shown in Fig. 27. The Deltaviex model, 
seen here with its pilot aboard, is fully airborne, but a 
funicular suspension—improved at ONERA—reduces 
the number of degrees of freedom. This system now 
allows low-speed stability and maneuverability studies 
to be undertaken—studies in which there is a renewed 
interest because of STOL and VTOL problems. 


AERODYNAMICS 


The science of aerodynamics, with its attendant com- 
plications of compressibility, viscosity, and turbulence 
has been, and continues to be, the object of much activ- 
ity, both theoretical and practical, at ONERA, and has 
given rise to unceasing development of wind tunnels of 
all types as well as of their equipment, which has be- 
come more and more automatic. 

Thus, for example, eight years ago, and at about the 
same time as the NACA, though each was ignorant of 
the secret work the other was doing, ONERA was com- 
pleting a transonic device which was immediately used 
to throw light on numerous problems connected with 
French aircraft, then just beginning to cross the noto- 
rious sound barrier. 

Of the great number of aerodynamic studies and of 
the aerodynamic facilities developed at ONERA over 
the last ten years or so, I can obviously quote only a 
few, though this is inevitably unfair to the rest of the 
work done. 

I shall therefore mention only three of these, each of 
different character. 

The first example relates to the fundamental problem 
of flow separation on a three-dimensional body. By 
making certain judicious assumptions in conjunction 
with a principle known as ‘‘prevalency’’—a generaliza- 
tion of an idea put forward by L. Prandtl—ONERA en- 
gineers”*! have established a separation theory for the 
laminar boundary layer, whereby the line of separation 


is defined as the ultimate pos:tion of a streamline passing 
through a characteristic point of separation on the body. 

The transition problem was analyzed concurrently 
and, by combining two criteria of longitudinal and 
transverse stability,°” it became possible to make a num- 
ber of forecasts later borne out by experience. 

To check the theory of three-dimensional separation 
in a laminar boundary layer, calculations were made, 
particularly in respect to a 10 per cent thickness-ratio, 
symmetrical-airfoil delta wing of elliptic cross section. 

Fig. 28 shows that the calculated line of separation is 
consistent with the actual flow conditions observed in a 
hydrodynamic tunnel which I designed in 1950 and 
which is very widely utilized when complex phenomena 
must be detected and understocd. 

Asa second example of fundamental research work, to 
enable the valuable nature of these detections to be ap- 
preciated, a film sequence will show us the flow around 
a thin delta wing the angle of attack of which is pro- 
gressively increased, then decreased and reduced to 
zero. The first frames concern the flow outside the 
boundary layer of the upper surface seen from above. 
They show the formation of vortex sheets rolling up 
into scrolls from the apex on both leading edges, in ac- 
cordance with my theory on rotational subsonic flow 
for wings of this type.*! 

Then there is a progressive decrease in delivery of the 
colored liquid used to show up the streamlines. We shall 
see the flow pattern changing and turning into the 
boundary-layer streamlines close to the wing surface.* 

The next sequence shows a further example of these 
visualizations, this time relating to the fow along an 
ellipsoid set at zero incidence; with a change from 
outer flow to flow along the body, the theoretical 
“separation line’’ already mentioned is seen to settle 
down fairly well. 

My third example deals with recent theoretical in- 
vestigations on very new effects of hypersonic aerody- 
namics. 

A good intrinsic and physical definition of these ef- 
fects might well be to state that they take up where 
nonlinear phenomena—in particular, the entropy varia- 
tion linked with local dissipation of energy—cease to 
play a secondary role, a condition usually correspond- 
ing to Mach Numbers greater than 5. 

In such flows, if they are stationary, waves of en- 
tropy variation are propagated—in the Hugoniot sense 
of the word—along the streamlines, while expansion 
and compression waves are propagated along the Mach 
lines in the same way (Fig. 29). 

The approach theory which I made known twelve 
years ago” consists in neglecting on the body, when this 
is permissible, the effect of compression waves reflected 
by the front shock wave. In this way very simple cal- 
culations become feasible, and in a recent study one of 
my collaborators was able to specify the second-order 
approximation. 

* At this point in the Lecture, two short film sequences were 
shown, 
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Parallel to this, he made a deep investigation of the 
structure of the ‘“‘shock layer’’ and rigorously demon- 
strated that in this layer the streamline pattern is virtu- 
ally superposable on the geodetic pattern of the surface 
of the body.*® 

This theoretical study has led to a most instructive 
analysis of the ‘‘shock layer’’ in the case of a flat plate 
parallel to the airstream, having a leading edge of finite 
thickness.” 

In the conventional representation given in Fig. 30, 
this leading edge behaves as a “‘source of entropy”’ 
which influences a partial and limiting layer of low 
density and considerable thickness. 

It is around this layer that the fluid flows away un- 
affected by the leading-edge effect, and the flow is so 
to speak ‘‘aware’’ of the obstacle only by virtue of this 
“shock layer,’’ which is generated by an extremely 
narrow strip of upstream flow of given breadth and as- 
tonishing expansion and rarefaction characteristics. 


STOL VTOL 


Despite the great progress made since the flight at 
Kitty Hawk, birdlike zero-length take-off and landing 
have not yet been achieved by conventional aircraft. 
Although the helicopter has been quite successful in 
certain uses, there is no doubt that part of aeronautical 
research must still be devoted to STOL and, ultimately, 
VTOL designs. 

Several different solutions of this problem are being 
experimented with in your country: either the aircraft 
is vertical to the ground, then tilts over to the horizon- 
tal and returns to a vertical attitude while in flight, or 
else the engines—-with or without the wings-—swivel 
through 90° after take-off and again prior to landing. 
That these methods are claiming equal attention indi- 
cates that none is superior to the others. 

In connection with STOL, ONERA has carried out 
and published extensive work on the controlled and 
forced circulation by boundary-layer suction and blow- 
ing at wing leading or trailing edges. The most recent 
of these works have mainly been devoted to trailing- 
edge blowing, associated or not with the action of a 
flap, which is itself blown at its hinge. 

A theory has been elaborated®**-® which suitably 
links up experimental results and which is quite literally 
“lit by flow visualization.* 

Applications leading to a quite notable reduction in 
take-off and landing runs can be deduced from this 
theory, provided that blowing air is supplied in reason- 
able abundance, which may well lead to the use of spe- 
cial turbo-blowers on aircraft. 

I shall only mention an interesting application on 
our Deltaviex experimental flying model. The take-off 
weight of this tiny machine (Fig. 31) is 2,080 Ibs. and 
its lifting area is 36 sq.ft. only; the wing span does 
not exceed 124 in., sweepback at the leading edge be- 
ing 70°, while maximum static thrust of its turbojet is 
940 lb. It is probably the jet plane with both the 
sharpest sweepback and narrowest wing span in the 
world. 
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Piloting this machine proved very difficult, as its 
very low roll inertia gave it a response time shorter in 
many instances than that of its excellent pilot. How- 
ever, by using a blowing system at the wing-flap trail- 
ing edges involving no more than a 2 per cent bleed 
from the engine air flow and by subordinating the action 
of these air jets on each wing to the angle and rate of 
roll, ONERA succeeded in achieving a forced stabili- 
zation which restored normal piloting characteristics to 
the aircraft, as you will be able to see now by watching 
this uncommon aircraft take off (short film sequence). 

To return to the problem of STOL and VTOL, I 
would like to mention the result of a study carried out in 
accordance with a personal concept of mine aiming at 
utilizing in various ways, for lift and thrust, a dual flow 
turbojet, which one might term a turbojet with ‘‘poly- 
valent secondary flow.’’ The flow is given moderate 
mechanical compression, the magnitude of which 
is adjusted to flying speed. Compression, at a maxi- 
mum on take-off, gradually diminishes until it is zero 
at maximum flying speed—envisaged at Mach 2 to 2.5 
at least. This variable compression is achieved by us- 
ing pivoted guide-vanes for the two-stage compressor 
and for the single-stage turbine of the free low-pressure 
rotor, the latter being left to windmill after low-pres- 
sure compression has been cut off. 

The secondary flow, the output of which is of the or- 
der of 1.5 times that of the central flow, can also be used 
(Fig. 32) for blowing on wings and tail units during cer- 
tain stages of take-off and landing. 

When such blowing is resorted to, the secondary flow 
is almost wholly diverted to nozzle combustors used for 
lift, so located and set as to cause their thrust to be al- 
most exclusively transformed into lift and to contribute 
to trimming of the aircraft, their jets being discharged 
at an adequate distance not only from the engine air 
intakes to preclude the possibility of exhaust gas being 
inadvertently recycled, but also from any sensitive 
part of the aircraft. 

I shall not dwell here on the advantages nor on cer- 
tain technological difficulties inherent in this solution, 
which dispenses with the need for any tilting of the pro- 
pelling jets or of their exhaust gases but which will re- 
main a speculative solution until some jet constructor 
takes up the experiment, which is of course more likely 
to occur in this country than in Europe for obvious rea- 
sons. 

Two designs may be conceived, according to the basic 
diagrams in Fig. 33, with associated or separated flows. 
In these diagrams, the main sections of the two flows 
are shown in their true relative magnitudes, the only 
varying sections apart from the intake, which is not 
shown, being the throat and outlet of the propulsion 
nozzles, which must differ somewhat. 

As in the United States, direct lift by jets has already 
been achieved in Great Britain, France and the USSR— 
in an experimental form only: it is true, but which none- 
theless shows the way for the future. This type of per- 
formance requires the vertical thrust to exceed the all- 
up weight by at least 15 per cent. 
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One immediately realizes that an equivalent conven- 
tional jet engine would use such thrusts for ultra-accel- 
erated—and hence very short—take-offs. Thus, from 
a comparison with the concept of the jet engine with 
my polyvalent secondary flow, one can estimate that 
conventional propulsive jet engines with the same fron- 
tal area would have a total thrust that would provide 
an initial acceleration of 1.5g and an extremely short 
take-off run—though not necessarily as short a land- 
ing run—-and would subject passengers of a commercial 
aircraft to a degree of strain that would probably turn 
out to be prohibitive. 

Fig. 34 compares this form of STO to a vertical riser 
conforming to my solution, for the case of a civil trans- 
port airplane capable in both cases of cruising at JJ = 
2.5 over stage lengths of the order of 3,000 kilometers, 
carrying a payload attaining 10 per cent. Of an all-up 
weight of 90,000 Ibs. such an aircraft would be powered 
by six jet engines with polyvalent secondary flow, and 
of 50 in. maximum diameter. Without laboring the 
point unduly, let it be said that transport by supersonic 
aircraft with zero-length take-off and landing charac- 
teristics is a solution offering numerous advantages for 
the future, and that its advent is perhaps certain. 

Unavoidably, jet-achieved VTOL, with provision for 
conventional emergency landing, involves a considerable 
penalty in power, but when it is combined with the ca- 
pacity for flying at supersonic speed—as much for mili- 
tary aircraft as for civil transports—does it not justify 
a search for appropriate powerplants? 

In presenting this subject among others to you, I 
hope I have been able to show you that although Jules 
Verne’s anticipations on the voyage from the earth to 
the moon are now giving rise to considerable and won- 
derful efforts in science and technique in the country of 
the Wright brothers, other scientists and technicians in 
Verne’s own country, without overlooking space flight 
research, are continuing to work for the great and very 
useful progress that aviation can still accomplish by 
following the way so brilliantly opened up by the Wright 
brothers. 
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Discussion of the Lecture 


Hugh L. Dryden, Deputy Administrator, NASA: 
The lecture just delivered by our illustrious French 
colleague, whose personal friendship I enjoy and value 
highly, is a fitting memorial to Wilbur and Orville 
Wright in its presentation of the current valuable 
contributions to aeronautical knowledge of an out- 
standing scientific leader and his collaborators. On 
your behalf and mine, I express appreciation to Pro- 
fessor Roy for his lecture and ask him to extend our 
compliments also to the ONERA group on their 
accomplishments under his guidance. 

Professor Roy speaks briefly of the new field which 
he terms ‘‘Cosmonautics’’ and we in the U.S. call 
“Space Science and Technology” or, in some quarters, 
“Astronautics.” Had we had the courage to adopt 
his terminology, the new National Aeronautics and 
Space Administration might have retained the old 
NACA initials. Though we like many others dream 
of travel to the stars and throughout the cosmos, 
we felt that for some time to come our travels would be 
confined to our solar system and so avoided the grander 
implications of “‘cosmonautics.” We in NASA agree 
fully with the stress placed on the continuing im- 
portance of aeronautics research. 

I have been greatly impressed by the broad attack 
of ONERA, its fundamental research in materials, 
dynamics of structures, and aerodynamics, its novel 
applications of physical principles and ingenious 
instrumentation, and its bold experimentation. Its 
work has supplemented and complemented the work of 
NASA and similar research organizations in other 
countries with unique contributions, demonstrating 
anew that creativity and originality are not confined 
to any one group or country. 

Being originally an aerodynamicist, I was especially 
interested in the work on transition and flow separation, 
particularly in the study of flow separation on three- 
dimensional bodies. The work on flow visualization 
and hypersonic theory is noteworthy. 

The extensive work of ONERA on controlled and 
forced circulation by boundary-layer suction and 
blowing has long been known and admired. We are 
intrigued by Professor Roy’s description and film of 
the Deltaviex experimental flying model, an excellent 
illustration of the bold experimentation to which I re- 
ferred previously. 

Finally, I would call your attention to the many 
evidences throughout the paper of the personal contri- 
butions of Professor Roy himself, ranging from hy- 
personic flow theory to the analysis of propulsion 
systems for STOL-VTOL aircraft and his concept of 
the turbojet with “polyvalent secondary flow” as the 
most suitable type of power plant for these applications, 
particularly for supersonic transports capable of vertical 
take-off. This personal scientific productivity of 
the leader is a key element in the notable record of 
ONERA contributions. 


A. H. Flax, Vice-President—Technical Operations, 
Cornell Aeronautical Laboratory, Inc.: No one can 
fail to have been impressed, not only by the technical 
breadth and depth of Professor Roy’s lecture, but also 
by the genuine enthusiasm with which he delivered it. 
In this respect, he is a worthy successor to Professor 
Panlevé, his countryman, to whom he alluded in his 
historical introduction. If I may add a note to this 
introduction, it is recorded' that on taking his first 
flight Panlevé was in high spirits and was waving and 
gesturing gaily to the assembled crowds before the 
take-off. In the course of so doing, he inadvertently 
pulled a rope Wilbur Wright used to stop the engine, 
thus adding a light touch to the proceedings. As 
Professor Roy has told us, however, they ultimately got 
off and made a record-breaking flight. 

It is difficult to select, from so overwhelming a 
display of scientific and technical accomplishment as 
we have had put before us here, some particular points 
for detailed comment—one can hardly do justice to 
even a small number of the subjects mentioned. I 
shall limit myself to a few points bearing on fields of 
recent interest to myself and my colleagues. 

The materials field steadily assumes increasing 
importance in all kinds of engineering endeavor. The 
need to work at higher temperatures, in nuclear and 
cosmic ray environments, and under high mean stress 
coupled with severe fatigue loadings has engendered 
seemingly impossible requirements for materials. The 
empirical routes of materials development, dating from 
the Bronze Age, although they have served well in the 
past, seem no longer sufficient for current technology. 
Application of the rational techniques of modern 
physics and mechanics to materials problems certainly 
seems to have been slower than the importance of 
the problem warrants. The study of X-ray spectra 
produced by electron beams for analysis of materials, 
which is discussed in the paper, can be a valuable step, 
not only for introducing new methods of microanalysis, 
but also for more intimately acquainting materials 
engineers with fundamental physical processes in 
solids. Hopefully, this could lead to larger benefits 
in the understanding of material properties. A 
development with parallel implications in this country 
has been the recent application of electron paramagnetic 
resonance measurements to studies of polymerization 
and depolymerization of plastics.? It is clear that there 
is much more of this sort of thing to be done, and that 
it will be done in many countries. 

The very high frequency fatigue machine reported 
on in the paper is an interesting development. Yet one 
cannot avoid suspecting that there must be high- 
frequency effects on fatigue properties, at least for 
some materials, so that such a machine cannot be used 
uncritically. On the other hand, the discovery and 
interpretation of such effects cannot but provide 
valuable information on the basic microscopic prop- 
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erties of materials, and is all to the good. The 
comprehensive survey of high-frequency phenomena 
in metals by Zener*® suggests that this may become a 
most important tool in the hands of metal technologists. 

The studies of crackling of cobalt during its phase 
changes have been paralleled in this country and 
Germany by investigations of sounds emitted by metals 
under stress.“ ° In the latter studies, attempts were 
made to analyze the results on a quantitative basis, 
with very limited success thus far. Yet here again is 
evidence of whole worlds of phenomena, either partly 
or wholly outside the range of human senses, but 
easily observable with modern instruments, entirely 
unexplored as to their engineering significance. 

It is interesting that flight flutter testing, long viewed 
with fear and suspicion, following Von Schippe’s ill- 
fated experiments, has now been generally accepted in 
France as it has been in this country. The French 
equipment for such work, displayed in the Figures, is 
indeed impressive even in the light of our own very 
considerable accomplishments in this field. 

The energy method for the prediction of flutter speeds 
when a large number of degrees of freedom are involved, 
if I understand it correctly, cannot be more than a first 
improvement to the primitive energy method, since the 
modes are considered only in pairs for determination of 
relative amplitude and phase. Although such a method 
may have important specific applications it does not 
appear to be possible to violate the fundamental axiom 
of flutter analysis—that is, there is no royal road to the 
answer. When making flutter analyses by any method 
whatever,one always feels sure he is missing something, 
and the more approximate the method the worse this 
feeling becomes. The only relief for this unhappy 
feeling in difficult cases is recourse to wind-tunnel or 
flight tests. 

In wind-tunnel techniques, I am willing to make 
something of an international trade with Professor 
Roy. I will buy his magnetic balance, and will sell him 
in return some of our automatic digital read-out 
equipment for pressure measurement, which eliminates 
the manometer photographs from the picture entirely. 

The progress in three-dimensional boundary-layer 
theory and particularly its application to leading-edge 
separation on delta wings was especially interesting to 
me, since the subject has engaged our attention at 
Cornell Aero. Lab. for over 10 years now. Our first 
attempts to describe forces and moments on the basis of 
cross-flow theory were empirical and tentative, but 
remarkably successful for so small an effort. Subse- 
quent efforts® to refine the concepts quickly led to the 
need for a choice between amassing vast quantities of 
empirical data and working up a detailed picture of 
the viscous flow around the leading edge. Under the 
press of other programs we actually did neither, but 
have favored the latter course. It is therefore grati- 
fying that ONERA, and particularly Professor Roy, 
have undertaken this investigation and have made 
substantial progress. It is impossible to leave this 
subject without remarking on the flow visualization 
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technique employed in these tests. Here one feels 
that the French talent for the esthetic component in 
science and technology has come to the fore. I have 
seen many methods for flow visualization, each having 
its own special utility and some producing rather 
striking pictures. But rarely have I seen such beautiful 
pictures which at the same time provided valuable new 
information on the details of the flow. 

Viewing the lecture as a whole, it is difficult to escape 
the feeling that there is far less information exchange 
and cooperation between ourselves and the French 
than there should be. In some respects, the connec- 
tions have been good, particularly where AGARD has 
been involved, but in other things we are working on 
the same problems, but might as well be on different 
planets. The barriers of geography, language, national 
pride, and military security are formidable, but not 
insurmountable. No one who has heard this lecture 
cali go away with the feeling that we have a monopoly 
on scientific and technical ingenuity this side of the 
Iron Curtain. It should be clear to all that, in the 
difficult and exciting days which are ahead, all the 
talents which we and our friends can command will be 
needed. This paper has not only exposed what has 
been done in the past but has also pointed up how 
valuable the French contribution might be in the 


future. 
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John J. Green, Defence Research Member, Canadian 
Joint Staff, and Defence Research Attaché, Canadian 
Embassy: I was most pleased and honored to be 
invited to join with Dr. Dryden and Mr. Flax in dis- 
cussing this Twenty-Second Wright Brothers Lecture. 
Professor Roy and I have been associated in the work 
of AGARD ever since its creation in 1951, and I like 
to think that some of the work of ONERA which 
Professor Roy has discussed today had its inspiration 
in the panels and committees of AGARD. If this is 
indeed the case, then it is most appropriate that 
Professor Roy should be giving the Wright Brothers 
Lecture this year because much of the inspiration 
generated by AGARD stems from the great support 
given to it by the United States. 

This is the first time that a Frenchman has been 
chosen to deliver this most important lecture. This 
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too is appropriate, for when the great accomplishments 
of the Wright brothers became known, France was 
soon in the forefront of development and achievement, 
and the early history of aviation is bright with the 
names of famous French pioneers. These men too 
were no doubt inspired by the triumph of the Wright 
brothers and especially by the demonstration flights 
and record-breaking performances of Wilbur Wright 
at Le Mans in 1908, to which Professor Roy has referred 
in his lecture. 

The history of French aeronautics reveals much more 
than the flair for aircraft design possessed by French 
engineers. There is an unorthodoxy—an individualism 
if you will—which suggests the search for elegant 
solutions free of compromise. The lecture we have 
heard today suggests very strongly to me that this 
same imaginative, unconventional approach is manifest 
in French aeronautical science. Throughout Pro- 
fessor Roy’s talk there are many examples of this 
search for new techniques to replace traditional meth- 
ods. At the risk of being considered indiscriminate, 
I will mention one or two examples which appealed 
to me particularly. The elegant approach employed 
in flight research into aircraft structural vibration and 
damping is most commendable, as are the various 
aerodynamic techniques for the investigation of 
problems both in the wind tunnel and in flight. I 
was intrigued, too, by the unique magnetic suspension 
for wind-tunnel models, which perhaps should be 
named the “hobo’’ method—it has no visible means 
of support. 

Professor Roy referred to the “unceasing develop- 
ment of wind tunnels of all types’’ at ONERA. I 
should perhaps caution him here that he might find 
himself in the position in which G. W. Lewis of the 
NACA found himself a number of years ago. The 
story goes that Dr. Lewis was escorting a distinguished 
British visitor around Langley Field, showing him 
one wind tunnel after another. Since the visitor made 


at 2-year intervals in future. 
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no comments and maintained what was becoming an 
awkward silence, Dr. Lewis at length ventured to 
ask him what he thought of it all, to which, so the 
story goes, the visitor replied, ““A vulgar display of 
wind tunnels.”’ 


It is significant to me that Professor Roy's lecture 
concludes with a discussion of some of the French 
work on STOL/VTOL aircraft, for, to my mind, 
practical progress in this aspect of flight is the greatest 
challenge today to the aeronautical engineer. There 
is no shortage of power, but to apply this efficiently, 
sagaciously, and perhaps quietly in order to shorten 
the take-off and landing run of aircraft would have a 
more profound effect on aviation than almost any 
other development has had in a history brief in years 
but studded with brilliant achievement. All countries 
with aeronautical capability seem to be taking up the 
challenge, and I am glad to say that my own country, 
Canada, is making some contributions. On December 
10, 1958, a Canadian designed aircraft, the De Havil- 
land ‘‘Caribou,’’ the fully loaded weight of which is 
about 25,000 Ibs., was demonstrated at Fort Myer— 
landing and taking off from the same place that Orville 
Wright had flown from in September, 1908, to which 
Professor Roy referred in his opening remarks. Your 
newspapers said it “stopped on a dime.’ This may 
be a slight exaggeration. Flying at about 22,000 Ibs. 
gross weight, it appeared to me to take-off and land 
in a distance of about 250 ft. 


It would seem that the short take-off and landing 
aircraft has arrived and that it now remains for us to 
press on towards even more ambitious performances. 
There is room here for great ingenuity and much 
research effort, and neither of these are lacking on this 
continent or in Europe. The work at ONERA and 
elsewhere in France is excellent in this regard. 


It is a pleasure to thank Professor Roy and to 
compliment him on his distinguished address. 
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SUMMARY 


Assuming the temperature of the upper surface of a rectangular 
sandwich plate to be 77°F. higher than that of the lower one, 
equations are developed for bending and torsional moments 
and transverse shear forces. Special attention is paid to the re- 
quired edge reactions. 

Except for cases where all edges are clamped or all edges are 
simply supported, transverse shear forces occur, so that the shear 
deformation of the core has to be considered. Plates with two 
opposite edges simply supported and the other two free or clamped 
are dealt with. Also, the case where three edges are simply 
supported and the fourth edge is free is computed. Graphs of 
some numerical results are presented. 


SYMBOLS 

a = length of plate in the X-direction 

b = length of plate in the Y-direction 

m = integer number 

Pm = given by Eq. (42) 

t = total thickness of sandwich plate 

to = thickness of core 

ty = thickness of faces 

w = deflection, positive if downwards 

= coordinates 

As, Bu; Cu, Din, Em: Fe = Constants 

B, = upward force to be exerted by cor- 
ner bolt 

D = flexural rigidity of sandwich plate 

PO. és = shear stiffness of sandwich plate 
with respect to shear force 
Qn. t. 

E = elastic modulus of faces 

G = modulus of rigidity of core 

I = moment of inertia 

Ma. ti = bending moment per unit width 

Muy = torsional moment per unit width 

Miu = moment per unit width in Sec- 
tions I and IT 

MM, u = average moment per unit width 

) a = transverse shear force per unit 
width 

R = concentrated reactions at corners 
of plate 

’ is = excess temperature of upper with 


respect to lower surface if ther- 
mal gradient is 7/t; equal to 
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Thermal Stresses and Deflections in 
Rectangular Sandwich Plates 


t 
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Bell Aircraft Corporation 


[t/(te + ty)|Ta if gradient is 


arbitrary 

i = difference between average tem- 
peratures in the two faces 

T2 = tensile force in attaching bolt 

Vay = statically equivalent reaction 
per unit length of edge 

x, ¥ = coordinate axes 

= functions of y only 

a = coefficient of thermal expansion 

Qm = mzb/2a 

Bm = mrb/a 

Ym = pmb/2 

5m = pmb 

Am = mr/a 

v = Poisson's ratio 

Subscripts 
SF refers to the X- and Y-directions 
n,t refers to the directions normal 


and tangential to edge 


INTRODUCTION 


obs TEMPERATURE of the upper surface is assumed 
to be T°F. higher than that of the lower one, but 
it is assumed not to vary in the X- and Y-directions. 
The thermal gradient in the faces is assumed as 77/1. 
The gradient in the core is immaterial, since the normal 
stresses a, and o, in the core are neglected, as are those 
perpendicular to the middle surface. In cases where 
the transverse thermal gradient in the faces is different 
from 7/t, T can be considered to be equal to [t/(t, + 
where is the difference between the average 
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The plate is assumed 
The proper flexural 


temperatures in the two faces. 
to be free to expand in its plane. 
rigidities of the faces are neglected. 

The sign convention is similar to that in reference 1. 
The cases where all edges are clamped or simply sup- 
ported lead to similar results as for homogeneous 
plates,’ where deflections from transverse shear are 
neglected, since here the transverse shear forces appear 
to be zero. However, these results are considered in 
more detail, and it is shown that for all edges simply 
supported the concentrated corner reactions R are 
infinite. A method is given for computing the actual, 
finite force in the attaching corner bolt. In the other 
three cases, with two opposite edges simply supported 
and the other two edges free, or clamped, and that with 
three edges simply supported and the fourth one free, 
the shear deflections are taken into account, and 
formulas are derived for all important magnitudes, es- 
pecially the edge reactions and the concentrated corner 
reactions. 


PLATE CLAMPED AT Four EDGES 


A similar case for homogeneous plates is treated in 
reference 1, p. 54. The equally distributed clamping 
moments at all sides are 


M, = aTD(1 + »)/t (1) 


so that the curvature produced by the nonuniform 


heating is eliminated over the entire plate. This leads 
to bending moments 
M, = M, = M, (2) 


over the entire plate (see Fig. 1), and twisting mo- 
ments 1/,, and transverse shear forces Q, and Q, equal 
to zero. 

From reference 2, for a clamped sandwich plate 
three boundary conditions have to be satisfied. As- 
suming the sandwich plate to be supported over the 
entire thickness, they are 


(Ow/On) — (Q,/De,) = 0, Q:/De, = 9 (3) 


Since Q, = Q: = 0, for the homogeneous plate these 
conditions are also satisfied, and, since over the entire 
plate Q, = Q, = 0, a clamped sandwich plate will be- 
have exactly as a homogeneous one. The bending 
moments are given by Eqs. (1) and (2), where D is now 
the flexural rigidity of the sandwich plate (shown in 
cross section in Fig. 2), so that 


D = + t,)?/2(1 — (4) 


w= 0, 


PLATE SIMPLY SUPPORTED AT FouR EDGES 


The same case for a homogeneous plate is treated 
in reference 1, pp. 176-180. It is found that for the 
deflection w the following differential equation holds: 


= (0°w/dx?) + (0°w/dy?) = aT(1 + v)/t (5) 
Hence, the transverse shear forces are 


Q, = —D(0/dx)V*w = 0 
Q, = —D(0/dy)V*w ot (6) 


II 
Il 
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For a simply supported sandwich plate, the three 
boundary conditions are 


M, = 0, QO, = 0 (7) 


which are also satisfied for the homogeneous plate. 
Hence, for the sandwich plate, 0, = Q, = 0 also, so 
that it behaves as a homogeneous plate. Therefore, 
from Eq. (j) on p. 178 of reference 1, 


w= 0, 


sin (mmx/a) 


m=1,3,5,... m* 
cosh (mmy/a 
E osh (mm} (8) 
cosh a, 


and, from Eqs. (k) on p. 179, 
4aT(1 — v*)D 
_ 


M, 
at 
sin (mmx/a) cosh (mmy/a) 
m=1,3,5,... m cosh a, 
' (9) 


sin (max/a) cosh (mry/a) 


m=1,3,5,... m cosh a», J 


where D is now given by Eq. (4). The bending mo- 


ments are maximum at the edges, 


(M,), = (My),-9, = @T(1 — »*)D/t = 
Et,(t, t,)*aT /2t (10) 


For practical applications, it is important to find the 
required strength of the attachments that connect the 
edges of the plate to its supports. From Eq. (8), the 
twisting moments are 

nn Ow  4aT(1 — v*)D 
M,, = (1 — x 
= cos(mmx/a) sinh (mmy/a) 


m=1,3,5,... 


(11) 


m cosh a», 


M,, is maximum at the corners. For example, 


= [4aT(1 — v*)D/nt] X 


= (tanh a,,/m) (12) 

m=1,3,5,... 
so that ./,, becomes infinite in the corners. Assuming 
that the plate is supported such that the attachments 
(bolts) can only transmit forces perpendicular to its 
middle plane, concentrated reactions KX are exerted at 


7 
ied 
ut 
ns. a 
yse | 
Te 
nt 
+ 
ARS 


Fic. 3. 


the corners, for example, 
= 2(May), 9, (13) 


and these reactions are directed upwards. Hence, the 
forces R are infinite. From Eqs. (6) and (11), for 
example, for the edge x = 0, the statically equivalent 
reactions, continuously distributed along the edge and 
considered positive if acting upwards, are 


(Vz), = [Q, — (0M,,/Oy)]r0 = 
— [4aT(1 — v?)D/at} > x 


m=1,3,5,... 


[cosh (mmy/a)/cosh a,,| (14) 


Their minimum absolute value occurs at y = (0). For 
b = aand using Eq. (4), 


y-0 = + X 
[1/cosh (mm/2)| = 


m=1,3,5,... 


O.84[Et,(t, at| (15) 
At the corners, say at y = 0/2, 


—[2Et,/(t. + t,)2aT/at] X 


( Vz) y=a/2 


(cosh a,»,/cosh a,) (16) 
m=1,3,5,... 
which is infinite. 

Actually, in practical cases, R and V, will not lead 
to infinite forces in the attachments. For a square 
plate (6 = a) let the attachment be effected by bolts 
as shown in the neighborhood of the corner x = 0, 
y = b/2 = a/2 in Fig. 3. Then the force in the 
corner, which can be concentrated in the corner bolt 
and is compressive, has to be in equilibrium with the 
total downward load exerted by the tensile force in the 
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bolts from y = 0 to y = 0.45 6 = 0.45 a4 and from x = 
0.05 a to x = 0.50 a. Hence the upward force to be 
exerted by bolt 1 (Fig. 3) is 


0.454 50a 
where, since M,, = —My,,, 


0 


(My), 0, y=0 45a 


0.50a 0.50a 


(May), 0.500 (My), s00 
x=6.05a x=0.50a 


Hence, from Eqs. (11) and (17), 
4aT(1 — v*)D 
— 

mt 


| sinh 0.45mr 
m=1,3,5,... ™ cosh 


B, 


= cos 0.05m7 tanh 0.5m7 

m= 1,3,5,... m 

9.48a7T(1 — 


The tensile forces in the other bolts can be calculated 
assuming that they take the downward forces —(V,) 


r=0 


or (Vy), 9.02. For example, bolt 2 in Fig. 3 has to 
exert a tensile force 


V 


(Mry), y=0 (May), 0, y=0.40a (20) 


= 
II 


where /,, is given by Eq. (11). 

Also, /,, does not actually become infinite in the 
corners. The highest moments occur in sections I or II 
in Fig. 3. In section I the total moment is 0.025 
+/2aB,, or per unit width an average 


= 0.025+/2 aB,/0.1+/2 a = (1/4)B, = 
0.38[Et/(t. + ty?aT/t} (21) 
which is smaller than J/, or J, from Eq. (10). From 
Eq. (20), Tz = 0.22[Et(t. + t,)*aT/t]. Hence, the 
average bending moment in section II is, from Fig. 3, 


Muy = [0.05+/2aB, — 0.025+/2 a(2T>) ]/0.15+/2 a = 


Further sections get smaller average moments. The 
moments J/, and /, from Eq. (10) cause at the ends 
of section II bending moments 1/11, of half the value 
given in Eq. (10). Assuming a parabolic moment dis- 
tribution in section II, the maximum moment Wyqy po, 
in its center follows from the condition 


Mire + (2/3) (Mit mar — Mire) = 


or 


Maan = (1/2) (3 Mir Mi) = 2) x 
(1.29 — 0.25) [Et;(t. + (23) 
= 0.52[Et((t, + t)*aT/t] 
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being practically equal to the maximum bending mo- 
ments ./, or MW, from Eq. (10). 


PLATE SIMPLY SUPPORTED AT Two OPpposITE EDGES 
AND FREE AT THE OTHER Two EDGES 


In a similar way as for a homogeneous plate, for a 
plate clamped at all sides the equally distributed 
clamping moments would be given by Eq. (1). Fora 
sandwich plate D is determined by Eq. (4). To obtain 


the case where in Fig. 1 the edges x = 0, x = a are 
simply supported and the edges y = +0/2 are free, 
edge moments MW, = —aT7TD(1 + v)/t have to be 


superimposed along all edges, so that, for the case 
superimposed upon that of the plate clamped at all 
sides, for the edges x = 0, x = a the boundary conditions 
are 
M, = -—aTD(1+y)/t, w=0, Q, =0 (24) 
while, for the edges y = +0/2, 
M, = —aTD(1 + v)/t, =0, Q, =0 (25) 


The expressions, given in reference 2 for the bending 
and twisting moments, contain the shear stiffness Do, 
and Dg,. From reference 3, for a sandwich plate with 
an isotropic core, 


Do, = Da, = Gl(te + ty)?/te] 
Since t; < t,, this can be written as 


Do, = Do, = G[(t.? + 2t.t;) ‘te | = 
G(t, + 2t,) = Gt (26) 


Using this expression, from p. 27 of reference 2, 


M, 


— D} (d0°w/dx?) + v(d°w/dy?) — 
(1/Gt) [(0Q,/dx) + »(0Q,/dy) 
M, = + v(0°w/dx?) — 
(1/Gt) [(0Q,/dy) + v(2Q,/dx) |} 
M,, = — v)/2|D}2(0’w/dxdy) — 
(1/Gt) [(0Q,/d¥) + (2Q,/ax)}} 


Inserting Eqs. (27) in the equation of equilibrium 


(02M, /dx2) — 2(02M,, /dxdy) + (02M,/dy2) = (28) 


(27) 


one obtains 
= 0 (29) 


as for a homogeneous plate. Eq. (29) is in accordance 
with Eq. (79) of reference 4, since in the present case 
g=m, =m, = 0. 

In order to satisfy, besides Eq. (29), the first condi- 
tion (24), the deflection w is assumed as 


w= — [aT(1 + »)/t] [x(a — x)/2] + 


Yn Sin (30) 


where Am = (31) 


Eq. (30) satisfies the second condition (24). Atx = 0 
and x = a, from Eq. (24) 0°w/dy? = 0 and 0Q,/dy = 0, 
so that from Eqs. (27) and (30), at these edges, 


M, = —[aTD(1 + »v)/t] + (D/Gt) (0Q,/dx) (32) 


Since Q, is antisymmetrical with respect to the line 
x = a/2, let 


(33) 
so that 


=1,3,5 


Hence, at x = 0 and x = a, 00,/0x = 0, so that from 
Eq. (32) the first condition (24) is indeed satisfied. 
With 


m=1,3,5,... 


Vym SiN (34) 


the third condition (24) is also satisfied. 

Since the first term of Eq. (30) vanishes after applying 
the operation V‘ to it, insertion of Eq. (30) in the differ- 
ential equation (29) leads in the usual way (see refer- 
ence 1, p. 127) to the solution 


Yn = Am cosh + sinh Any + 
Cy sinh + DadAmy Cosh Any (35) 


From symmetry about the X-axis C = D = 0, so that 
w= —[a7(1 + v)/t] [x(a — x)/2] + 
(Am + BryrrAmy sinh sin 


m=1,3,5,... 


(36) 
To find expressions for Y,,, and Y,,, it is observed that 
QO, = (0M,/d0y) — (OM,,/0x) (37) 

Further, since the transverse load g = 0, 
(00,/ax) + (2Q,/ay) = 0 (38) 


Using Eqs. (27) and (38), Eq. (37) becomes, after re- 
arranging, 
Q, — — v)D/2Gt] [(0°Q,/dx?) + (0°Q,/dy?)] = 

— D(0/dy) [(0?w/Ox?) + (0?w/dy?)] (39) 


whence, after insertion of Eqs. (34) and (36), 


[4Gt/(1 — v) sinh Any (40) 


Letting Y,,, = e?” and taking account of the fact that 
Q,, and therefore Y,,,, is antisymmetrical about the 
X-axis, the solution of the homogeneous part of Eq. 
(40) becomes 


Yim = En sinh pry (41) 
where 
bm = + [2Gt/(1 — (42) 
The particular integral is found by assuming 
= Fp sinh \ny (43) 


which, after insertion into Eq. (40), gives 


Fn = —2m®DBy (44) 


be 
2 
0Q,/Ox = — 
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a 


so that 
Yym = Em sinh pny + Fm sinh 45 | 
Pmy y (45) > xX 
With Eqs. (44) and (45), Eq. (34) becomes 
QZ, = (E,, sinh — 
2m®DB,, sinh \,,V) sin (46) 
Inserting Eqs. (33) and (34) into Eq. (38) leads to b 
2. = > (1/Am) (d Vym/dy) COS Ap Xx 
m =1,3,5,... 
(47) 
= [Em(Pm/Am) cosh Pmy 

2m® DB» COSH COS ApX 
é The constants A,,, B,,, and E,, in Eqs. (36), (46), and ’ 
= (47) have to be determined from the boundary condi- Y 
| tions (25). From the third condition (25) and from Fic. 4. 

Eq. (46), 

= Expressing £,, in terms of B,, from Eq. (48), A, can be 
j Em = (2Am®D sinh a,/sinh ¥m)Bn (48) expressed in terms of B,,, which gives 

where An = {[(1 + »)/(1 — »)] — coth (50) 
EB am = Xm(b/2); Ym = Pm(b/2) (49) Finally, using Eqs. (27) and (38), the first condition 


From the second condition (25) and from Eq. (27), at (25) requires that, for y = 6/2, 
y = 6/2, — Df (0°w/dy2) + v(0°w/dx?) — 

2(0'w/dxdy) — (1/Gt) + (20,/2x)] = 0 = 
so that, after using Eqs. (36), (46), and (47), Using Eqs. (36) and (46) and expressing a7 D(1 + v)¢ 


as 
Wm? [(Am + Bm) sinh am + Byam cosh a,| = 
(1/GtXm) [Em(Pm? + sinh + v) /t = 
sinh Am | (4 faTD(1 + 1] > (ain (51) 


m=z=1,3,5,... 


results in the condition 


Am2((1 — v)Am + {2 + [2(1 — By) cosh am + Am2(1 — v)QnBm sinh ay — 

4 [1 — v)/Gt]Enpm cosh ym = [4(1 — v)/arm] [aT(1 + v)/t] 
a Since from Eqs. (48) and (50) E,, and A, can be expressed in terms of B,,, with the use of Eq. (42) this yields 

4 8a(1 — v*)T 

E + v + Al sinh 2am 2(1 V) Qm + coth Ym(cosh 2am 1) 


(52) 
Taking account of the original moments 1/, = 1/, = M,, from Eqs. (1), and using Eqs. (36), (38), (46), and (47), 
3 the resulting bending and twisting moments are 
M, = [aTD(1 + v)/t] — D} (0°w/dx?) + v(d0°w/dy?) — [(1 — v)/Gt] (0Q,/dx)} 
=D > — »)Am — cosh Amy + (1 — ¥)BmXm sinh Any} + 
[(1 — v)/Gt] (EnPm cosh — 2AmtDB,, cosh Amv)) sinA,x (53) 


M, 


laTD(1 + v)/t] D} (0°w/Oy?) + v(O?w/Ox?) + [(1 — v)/Gt] (00,/dx)} 
— YAn + Bn] cosh Any + (1 — »)BrXmy sinh Any} — 
m2=1,3,5,... 
[((1 — v)/Gt] (EnPm COSH — Cosh Amy)) Sit Ame (54) 
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N 


[(1 — — (1/Gt) [(0Q,/dy) + (0Q,/ax)}} 
> {Am?[(Am + Bm) sinh + ByAmy Cosh — 
13.5,... 
(1/2AmGt) [Em(Pm? + Am?) sinh pny — 4m®DBy» sinh CoS Ame (55) 


From Eqs. (25), at the free edges M,, = 0, so that at all corners the concentrated reactions R = +2 M,, are zero. 
Using Eqs. (47) and (55) the tensile forces per unit length along the edges x = 0 and x = a that have to be trans- 
ferred by the attaching bolts are, for example, for x = 0, 


(1 — »)D > Am®[(Am + 2Bm) cosh + sinh + 


2 {1 — v)D/Gt}rn? + 1} cosh Any — 
n(Pm/Am) {[(L — »)D/Gt)m? + 2} cosh pny 


or, using Eqs. (48) and (50), 
—(Vz)z=0 = D Brdm* ({5 V) Qm coth a, + v)Ddn2/Gt]} cosh + 


m=1,3,5,... 


(1 — v)Amy sinh Any — (2Pm Sinh am/Am sinh Ym) | [(1 — v)DXym2/Gt] + 2} cosh pny) (56) 


PLATE SIMPLY SUPPORTED AT Two OPPOSITE EDGES AND CLAMPED AT THE OTHER Two EDGES 


If, in Fig. 1, the edges x = 0, x = aare simply supported and the edges y = +6/2 are clamped, after assuming first 
all edges to be clamped, so that ./, and /, become equal to /, from Eq. (1), moments 1/7, = —aTD(1 + »v)/t are 
superimposed along the edges x = 0,x = a. Then, for the superimposed case, the boundary conditions along these 


edges are given by Eq. (24), and along the edges y = +0/2, from Eq. (3), 
w= 0, (Ow/dy) — (Q,/Gt) = 0, Q, = 0 (57) 


From the boundary conditions (24) along x = 0, x = a, and symmetry of w about the X-axis, w, Q,, and Q, are given 
again by Eqs. (36), (46), and (47). Then, from the last condition (57), 


Em = (2Xm'D cosh COSH Ym) Bm (58) 
From the second condition (57) one obtains 
An, = {[(2\n?D/PmGt) tanh Ym — Coth am — (2Xmn2D/Gt) — 1} Bm (59) 
and finally from the first condition (57), with the substitution (51), 


ta(l T ] 
(60) 


Amat (2m*®D/PmGt) [tanh Ym coth am — (Pm/Xm)| — 1} cosh am — ap) 
With these constants 1/,, M,, M,,, and (V,),~0 follow from Eqs. (53)—(56). Using Eqs. (58) and (59), from Eq. (56) 


= 


= D >> (1 v) [(2Am Gtpm) tanh Ym — | coth Am + (3 v)} cosh AmY + 
m=1,3,5,... 


(1 — v)Amy sinh — (2 cosh a,/cosh Ym) — v)DAm?/Gt] + 2} cosh pay) (61) 


The tensile forces per unit length along the edges y = b/2 that have to be transferred by the attaching bolt are, 
using Eqs. (46), (55), (58), and (59), 


(V,),-02 = [Q, — (OM,,/Ox) = —2D DY Brdm? Sinh Om Ay X (62) 


m=1,3,5,... 


The required downward concentrated forces RX at the corners, for example, for x = 0, y = 6/2, are 


(R) = = 4D Bm(Am*/Pm) cosh am tanh Ym (63) 


m=1,3,5,... 


From Eqs. (54), (58), and (59), the clamping moments along the edge y = 5/2 are 


(M,),-02 = [aTD( — »)/t] -D (An2Bm/sinh am) {2(1 — v) (Am®D/PmGt) tanh 7», cosh? ay, — 
m=1,3,5,... 


(1 — v)am + [1 + » — 2(1 — v) (A,?D/Gt)] sinh a, cosh am} sin Amx (64) 
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Fic. 5. Moment loading per unit width due to a difference T in face temperatures (sandwich panel with all edges clamped). 


PLATE SIMPLY SUPPORTED AT THREE EDGES AND FREE AT THE FOURTH EDGE 


In Fig. 4 let the edge y = b be the free one. Then the boundary conditions for the superimposed case are again 
given by Eq. (24) for x = 0,x = a. For the edge y = 0, 


M, = —aTD(1+»)/t, w=0, = 0 (65) 
and for the edge y = 3, 
M, = —aTD(1 + »)/t, Myr = 0, Q, = 0 (66) 


From Eq. (24) w can again be represented by Eq. (30), where Y,, is given by Eq. (35). Similarly, Q, and Q, can again 
be assumed as Eqs. (33) and (34). Since presently C,, and D,, are not zero, Eq. (39) now leads to 


(d?Vym/dy?) — Am? + [2Gt/(1 — v)D]} Vim = [4Gt/(1 — v) sinh Any + Dm cosh (67) 
instead of to Eq. (40), so that the solution of the homogeneous equation becomes 
Yom = Em sinh pny + Fn cosh py (68) 


where P,, is given by Eq. (42). The solution of the complete equation (67) gives 


| 
“Me 


SiN 
m=13,5,... 
(69) 


= [Em sinh pry + cosh pay — 2An®D(B, sinh Any + Dycosh Any) | sin Anx 


m= 1,3,5, 


QQ, = 


By inserting Eqs. (33) and (34) into Eq. (38), 


(1/Am) (d¥ym/dy) COS 
m=1,3,5,... 
(70) 
(1/Am) cosh Pmy + FuPm sinh 2m'*D(Bm cosh + Dn sinh Amy) | cos 
m=1,3,5,... 


The constants A», Bn, Cn, and D, in Eq. (35) and E,, and F,, in Eqs. (69) and (70) follow from the six boundary 
conditions in Eqs. (65) and (66). They are 
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4a(l — v*)T 
Du = 
(1 — v)Bn + [3 +v»+2(1 — v)A,,2(D/Gt)] sinh B,, cosh B,, — [201 — v)D/Gt] Anpm tanh 6, cosh? pm, 
where Bn = Amd; Pmb (72) 
An = 4aT(1 + v)/Am* at, B,, = 0 | 
= Dn} + v), (1 v) | Bos tanh Bm} [4a7(1 + v) at | tanh Buus (73) 


Fn = 2DmXm*D cosh B»/cosh 6», 


M,, M,, and M,, can be calculated from the second members of Eqs. (53), (54), and (55) where w is given by Eqs. 
(30) and (35) and Q, and Q, by Eqs. (70) and (69), respectively. For the tensile forces along the edges that have to 
be transferred by the attaching bolts, one obtains 


= G5 — » — (1 — tanh B, + [211 — v)DA,,2/Gt]} sinh Any + 


m=1,3,5,.. 


(1 — v)Amy CoSh — 2(Pm/Am) [C1 — v)Ddn2/Gt] + 2} (cosh B,,/cosh 6,,) sinh Pmy) 
[4a(1 — v*?)TD/at] (tanh 6,, sinh A,,v — cosh (74) 
m=1,3.5,... 
—(V,),-. = D An? — tanh B, + [2(1 — v)DA,,2/Gt] [(cosh B,,/cosh 6,) — 1]} D, — 


m=1,3,5,... 


[4a(1 — v*)T/),,%at] tanh Bn) sin Amx (75) 
The concentrated downward force FX at the corner x = 0, y = O is 


(R), 0, = 2( Mary) r=0, y=0 


Gt | cosh it) _ 8a(1 — v)TD tanh 
(1 — v)DX?,, | cosh 6, 


») 
= 2(1 — v)D ) (= tanh Bn — 
v 


(76 
at m=1,3.5, Xm 


Since at the free edge /,, = 0, at the corners x = 0, y= 
band x = a, y = b the forces R are zero. 


SOME NUMERICAL RESULTS 


Alexander Krivetsky of Bell Aircraft computed some 
numerical results and presented them in the chart and 
graph given in Figs. 5 and 6. Fig. 5 gives the bending 
moments Vf, = M, = M,, from Eq. (1) for a plate 
clamped at all edges, as a function of the difference in 


6*a:24" 


temperature 7 and the ratio ¢,/t. If the temperature 
varies linearly through the thickness ¢, the bending wi - 
stresses in the outer fibers can be calculated directly as =|260 
| 
o, = o, = (M,/I) (t/2) = 240 | 
‘ -- oll | 
aTD(1 + = aTE/2(1 — ») (77) 
If the thermal gradient in the faces differs from 7'//? 5 of / 
where 7 for this case is defined in the Symbols, the 3 | | 
stresses in the outer fibers are 3 
=| 160) FOR ALL EDGES CLAMPED 
S| 40 + | 10 
where 7) is the temperature difference between the a | | | 
outer fibers. | 
Fig. 6 shows the bending moments at the clamped 100) | | Lim 
edge from Eq. (64) for a square plate of 24 by 24 in., = 
as sketched in the insert, simply supported at two a sails clu, 
“1G. 6. Moments and stresses along clamped edge (plate 
opposite edges and clamped at the others, if f - 0.9 — simply supported at two opposite edges and clamped at the 
t; = 0.02in. and T = 200°F. It shows the distribution others). 


of the clamping moment 7, between x = 0 and x = 12 
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in. It also shows the distribution of the stresses o, 
from M, at the outer fibers, assuming a thermal gra- 
dient 7/t in the faces. At the corner x = 0, y = 12 
in., where o, is zero, obviously 7, = (1/2)aET7, so that 
with a = 0.000028, E = 107 psi and T = 200°F., o, = 
12,800 psi in accordance with Fig. 6. 

For a square plate, clamped at all edges, from Eq. 
(77) o, = 12,800/(1 — 0.3) = 18,300 psi along the 
entire length of the clamped edges (as well as every- 
where else), as indicated by the dashed line in Fig. 6. 
Hence, for the case with two opposite edges clamped, 
and simply supported at the other ones, as considered 
in Fig. 6, the maximum clamping moments are larger. * 
This is understandable, since canceling of the clamping 
moments at the edges x = 0 and x = a will cause excess 
clamping moments along the edges y = +0)/2 in order 
to keep the slopes Ow/Oy along these edges equal to 
zero. This effect will become smaller if one moves 
farther from the edges x = 0 and x = a, which explains 
why the maximum stresses o, occur close to the edges 
x = O and x = a and not in the center of the clamped 
edges. 

It should be noted also that, for a thermal gradient 
7 /t in the faces, in general, the stresses a, in the outer 
fibers are not completely given by the formula 


o, = (M,/D) (t/2) = [Et/2(1 — (M,/D) (79) 


M, depends on the curvature of the plate as far as it is 
caused by the deformation from bending only. From 
Eq. (27) this effect is given by that of the total curva- 
tures 


* This is analogous to what happens in the case of a square 
plate with edges a, under equally distributed transverse loading 
q, where, with all edges clamped, the maximum clamping mo- 
ment is —0.0513 ga?, whereas, with two opposite edges clamped 
and the other two simply supported, the maximum clamping 
moment is —0.070 qa?. : 


— [(0?w/dy?) + v(0*w/dx*) | 


diminished by the curvatures caused by the change of 
the shear deformation with x and y, 


—(1/Gt) [(02,/oy) + ¥(0Q,/0x) 
or, from Eq. (38), 
— »)/Gt](0Q,/0x) 


However, the curvature of the faces, of which the proper 
flexural rigidity is neglected throughout, is equal to the 
total curvature. Therefore, to the bending stresses 
from Eq. (79) excess bending stresses 


o, = [E/(1 — (t,/2) — »)/Gt] (0Q,/dx) = 
(Gt,/Gt) (0Q,/Ox) (SO) 


have to be added, where G; is the modulus of rigidity of 
the faces. However, in the present case, since at the 
clamped edge, from Eq. (57), Q,, and therefore 0Q,/Ox, 
are zero, a, from Eq. (80) is zero. 

If the thermal gradient in the faces differs from 7° 1, 
the stresses o, from Eq. (79) have to be diminished by 
stresses 


o, = Ea(T — T)/2(1 — »v) (81) 


where 7) has the same meaning as in Eq. (78) and 7 is 
defined in the Symbols. 
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The Tailored-Interface Hypersonic Shock 
Tunnel’ 


CHARLES E. WITTLIFF,* MERLE R. WILSON, ** anp ABRAHAM HERTZBERG*** 


Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


The ‘‘tailored-interface’’ hypersonic shock tunnel provides a 
means for producing the high Mach Number, high stagnation 
temperature flow conditions encountered in hypersonic flight. 
Various gasdynamic phenomena associated with shock tunnels 
are discussed, and experimental evidence of the successful ap- 
plication of this technique is presented. As an indication of its 
research application, the results of heat-transfer experiments 
on a hemisphere-cylinder model are presented and compared 
with theory. 


INTRODUCTION 


I RECENT YEARS considerable attention has been fo- 
cused on problems associated with hypersonic flight. 
This interest, for the most part, can be attributed to 
the high-temperature phenomena behind strong shock 
waves. At flight Mach Numbers above 6, the tempera- 
ture behind a normal shock exceeds 1,500°K. and in- 
creases to about 10,000°K. at Mach 22. In this tem- 
perature range, the thermal energy of the air is sufficient 
to excite the higher rotational and vibrational states of 
the molecules, cause dissociation and ionization, and 
produce chemical reaction of the constituents. Since 
these real gas effects may alter the flow field about a 
body, it is important to have an understanding of this 
interaction. 

In this relatively unexplored flow region, one of the 
primary needs has been for laboratory facilities cap- 
able of reproducing the high-temperature and high 
Mach Number conditions encountered in hypersonic 
flight. The shock tube offers a means by which the 
high temperature conditions can easily be created and 
maintained in a nearly steady state for brief periods of 
time.' In addition, the equilibrium state of the gas, 
which is homogeneously heated by the shock wave with- 
out contamination, can be precisely determined. These 
and other features of the shock tube have led to its 
widespread use in the fields of aerodynamic,?~‘* physi- 
and chemical research.°—!* However, the maxi- 
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mum flow Mach Number which can be obtained in air 
in the conventional shock tube is approximately three.'* 
This Mach Number limitation has, in the field of hyper- 
sonic research, restricted the application of the shock 
tube to studies of flow problems which are not Mach 
Number sensitive.” Thus, the conventional shock 
tube has been successfully used to simulate the very high 
stagnation-point heat-transfer rates encountered in the 
hypersonic re-entry of a ballistic missile. 

However, there are many hypersonic problems in 
which Mach Number simulation is necessary. In or- 
der to meet this requirement, development of the hyper- 
sonic shock tunnel began at the Cornell Aeronautical 
Laboratory in 1950 with the addition of a simple diver- 
gent nozzle to the downstream end of a conventional 
shock tube.'* This nozzle expanded the supersonic 
flow behind the shock wave to hypersonic speeds, so 
that flow Mach Numbers up to 10 with full stagnation 
temperature simulation were obtained for about | milli- 
second. However, in order to obtain this flow duration, 
a relatively long shock tube was necessary. Because 
of viscous attenuation in this slender tube, the strength 
of the shock wave was found to be much less at the en- 
trance of the nozzle than it was initially. Thus, total 
enthalpy was lost, and flow conditions in the test sec- 
tion became unsteady. Detailed results of these meas- 
urements have been reported in reference 19. 

A study of more sophisticated shock tunnel designs 
showed that an increase in testing time per foot of shock 
tube of 8 to 60 times that of the conventional shock tube 
could be obtained.””: *! Such a large increase in testing 
time per foot would permit use of a shorter shock tube 
giving the same, or longer, test period relatively free of 
attenuation effects. The ‘‘tailored-interface’’ method, 
which resulted from this study, has recently been in- 
corporated in the C.A.L. hypersonic shock tunnel. 

A description of this technique and experimental 
verification of its usefulness are presented in this paper. 
In addition, some of the gasdynamic phenomena which 
are especially important in the design and operation of 
the tailored-interface shock tunnel are discussed. To 
illustrate the research application of this shock tunnel, 
measurements of heat-transfer rates on a hemisphere- 
cylinder body are presented, and the results are com- 
pared with theoretical predictions. 


THE TAILORED-INTERFACE METHOD 


In the “‘tailored-interface’’ modification, the incident 
shock wave is reflected at the downstream end of the 
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the same as in the driven gas. This information plus 
the specific heat ratio of the driver gas are sufficient 
to determine, first, the temperature of the driver gas 
ahead of the reflected shock, and, finally, the initial 
temperature and pressure of the driver gas.’ ** Thus, 
when the initial pressure and temperature of the driven 
gas, the desired Mach Number, and the driver gas com- 
position are specified, the required driver gas pressure 
and temperature to achieve tailoring can be determined. 
These quantities have been calculated for helium and 
hydrogen driver gases, initial air pressures from 0.01 to 
10 atmospheres, and flight Mach Numbers from 6 to 22. 
The results, presented in Figs. 2 and 3, are plotted, for 
convenience, against both shock wave Mach Number 
and the flight Mach Number corresponding to an equiv- 
alent stagnation enthalpy at an ambient temperature 
of 218°K.; i.e., in the stratosphere. 

These Figures show that helium at room temperature 
(300°K.) will tailor at an equivalent flight Mach Num- 
ber of 7.5, and room-temperature hydrogen will tailor at 
Mach 10.6. Tailoring at intermediate Mach Numbers 
can easily be achieved using a mixture of these two 
gases. To obtain stagnation enthalpies corresponding 
to flight Mach Numbers above 10.6 using the tailored- 
interface technique, a heated helium or hydrogen driver 
gas is necessary. This heating can be accomplished 
by various methods; for example, by an electrical heat- 
ing unit or a heat exchanger, by an adaptation of the 
double driver or ‘‘buffer’’ gas technique,*~*” by com- 
bustion,' or by an electric arc discharge.! Using helium 
heated to about 3,000°K. or hydrogen at 1,000°K., the 
tailored-interface shock tunnel should permit tempera- 
ture simulation at flight Mach Numbers up to 22. 


Testing Time 


The maximum available testing time is determined 
by the interaction of the reflected shock wave and the 
unsteady expansion wave propagating downstream 
from the primary diaphragm station, provided the 
driver is long enough so that the expansion waves prop- 
agating upstream into the driver do not overtake the 
downstream propagating waves after reflecting from 
the end of the driver (Fig. 1). Calculation of the test- 
ing time, based on the wave diagram in Fig. 1, indicates 
values approximately 8 to 12 times greater, depending 
on the Mach Number, than can be achieved with a 
conventional nonreflected shock tunnel of the same 
length. 

By placing a supersonic convergent-divergent nozzle 
at the primary diaphragm station, the supersonic por- 
tion of the unsteady expansion fan (Fig. 1) can be re- 
placed by a steady flow expansion and the testing time 
thereby increased to the entire time it takes the re- 
flected shock wave to travel from the hypersonic noz- 
zle to the diaphragm station and return (Fig. 4). This 
“steady”’ configuration provides a test period from 14 
to 60 times greater than the nonreflected method, de- 
pending on the Mach Number. However, the steady 
configuration requires that the area ratio of the steady 
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(steady configuration ) 


flow nozzle at the diaphragm station be changed for 
each Mach Number. The required driver gas pres- 
sures and temperatures for tailoring are approximately 
the same as those for the unsteady configuration (Figs. 
2 and 3). 

The variation of testing time per joot of shock tube 
length with equivalent flight Mach Number is presented 
in Fig. 5 for these two configurations. It is important 
to note that, although the testing time decreases appre- 
ciably with increasing Mach Number, the number of 
model lengths of flow remains relatively constant be- 
cause of the increased velocity. For example, in the 
C.A.L. shock tunnel using the unsteady configuration, 
over the Mach Number range from 6 to 18, the number 
of model lengths of flow decreases from 74 to 38, for a 
6-in. model, while the testing time decreases from 6 to 1 
milliseconds. 


Test Section Conditions 


The test section Reynolds Number per foot for flight 
Mach Numbers from 6 to 22 is presented in Fig. 6. 
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These values were calculated using real gas properties 
for air*~*! in thermodynamic equilibrium.* The test 
section conditions shown apply to any tailored-inter- 
face shock tunnel when stagnation enthalpy and Mach 
Number are simulated at an ambient temperature of 
218°K. Any shock tunnel will have some pressure 
limitation which will restrict its range of full simulation 
of Mach Number and altitude. Such a restriction is 
shown in Fig. 6 by the curve indicating a pressure of 
2,000 atmospheres behind the reflected shock wave, 
which is the pressure limit of the present C.A.L. shock 
tunnel. 


SPECIAL GASDYNAMIC CONSIDERATIONS 


There are a number of gasdynamic effects that have an 
important bearing on the design and performance of 
shock tunnels in general, and that warrant special con- 
sideration. Included in these are interaction of the re- 
flected shock wave and the boundary layer on the walls 
of the shock tube, mixing at the interface between the 
driver and driven gases, the Taylor instabliity of the in- 
terface, loss of stagnation enthalpy by radiation, con- 
vective heat transfer at the throat of the hypersonic 
nozzle, thermodynamic nonequilibrium in the nozzle, 
and adequacy of the available test time. Although 
these phenomena do not introduce any problems in the 
present operation and performance of the C.A.L. tail- 
ored-interface shock tunnel, some of them may present 


* The authors wish to thank C. Beckett for making available 
the latest NBS calculations of the equilibrium properties of 
argon-free air at densities up to 316 times sea level density.*! 
These tables were necessary in order to calculate the p; = 10 
atm. curve in Fig. 6. 
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difficulties if the range of operation is extended to higher 
Mach Numbers. 


Interaction of Reflected Shock and Boundary Layer 


The interaction of the reflected shock wave with the 
boundary layer behind the incident shock has been stud- 
ied analytically and experimentally by H. Mark,** who 
showed that a large scale shock-boundary-layer effect 
exists when low energy flow in the boundary layer does 
not have sufficient total pressure to pass through the 
reflected shock without separating. However, it has 
been observed that, at high densities in which the 
boundary layer behind the incident shock undergoes 
early transition from laminar to turbulent, this in- 
teraction is not very pronounced. The turbulent 
boundary layer feeds energy into the low energy portion, 
increasing the total pressure and inhibiting separation. 
In the C.A.L. tailored-interface shock tunnel, which op- 
erates at high densities, no effects of this shock-bound- 
ary-layer interaction have been found in the test air. 


Interface Mixing 


In all of the tailoring calculations it has been assumed 
that the interface between the driver and driven gases 
represents a discontinuity in temperature and density. 
In reality, the interface is a constant pressure mixing re- 
gion, the extent of which is believed to depend on the 
opening process of the diaphragm, the speed of the 
shock, and the geometry of the shock tube. Mixing 
between the hot air and cold driver gas would be ex- 
pected to reduce the amount of air available for testing. 
The experimental results obtained thus far, however, 
do not indicate any noticeable effect of this sort. 

It should be noted that interface mixing would be 
enhanced by any strong Taylor instability** relating to 
the deceleration of two gases of different density by the 
reflected shock wave. For the ‘‘tailored-interface,’’ the 
air density is always greater than that of the driver gas, 
and this represents a stable system. 


Radiation Losses 


If the almost stagnant air behind the reflected shock 
wave remains in the shock tube at a high temperature 
for longer than a few milliseconds, it is possible for it 
to lose a significant portion of its energy by radiation to 
the walls of the shock tube. For example, if the testing 
time were 20 milliseconds, which could easily be ob- 
tained in the steady configuration, the radiation loss at 
2,000 atm. would be about 6 per cent at a stagnation 
temperature of 6,000°K. and 15 per cent at 8,000°K. 
This loss would increase to about 40 per cent at 100 
milliseconds for 8,000°K., and would be even greater 
for lower stagnation pressures. The above values are 
based on an emissivity of unity selected on the basis of 
the available data.!* '¥ Calculations indicate, how- 
ever, that radiation losses never exceed 5 per cent of the 
stagnation enthalpy over the operating range of the 
C.A.L. tailored-interface shock tunnel. Any facility 
which heats a fixed mass of air in a closed chamber will 
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experience similar radiation losses. Hence, it is desir- 
able to obtain the data in the shortest possible testing 
time in order to minimize the uncertainty in the stagna- 
tion enthalpy. 


Nozzle Throat Heat Transfer 


Another heat-transfer problem that is important in 
the design of a tailored-interface shock tunnel is the 
convective heat transfer at the throat of the nozzle. 
The results of an analysis that considered both the com- 
pressible laminar*! and incompressible turbulent* 
boundary-layer cases indicate that the nozzle heating 
problem can be very serious for testing times greater 
than several milliseconds for both cases. For example, 
at a stagnation temperature of 6,000°K., an uncooled 
stainless steel nozzle would begin to melt after 11 milli- 
seconds. The testing time in the present tailored-in- 
terface shock tunnel is not long enough for the nozzle 
surface to reach the melting point, and, in fact, after 
some 200 tests no damage due to heating or erosion has 
been observed. 


Thermodynamic Nonequilibrium 


Since the temperature and density decrease rapidly 
through the hypersonic nozzle, it is possible that ther- 
modynamic nonequilibrium conditions may occur in the 
nozzle because of relaxation effects. Nonequilibrium 
conditions will exist if the local flow time becomes 
short compared to the local relaxation times for such 
phenomena as molecular vibration, dissociation, ioniza- 
tion, and chemical reaction between the different gas 
species. Criteria for predicting such behavior have 
been given by Penner* and Bray.” 

Although Penner’s criterion is for determining 
whether the exhaust gases from a rocket engine will be 
near thermodynamic equilibrium or nearly frozen after 
adiabatic expansion through a nozzle, it has been used 
to predict nonequilibrium effects in a hypersonic noz- 
zle.*®* More recently, Bray has proposed a criterion 
that appears to be more applicable to hypersonic wind- 
tunnel nozzles having large expansion ratios. Applied 
to the case of oxygen dissociation, for example, Bray’s 
criterion states than nonequilibrium, or ‘‘freezing,’’ will 
occur when the equilibrium rate of change of the degree 
of dissociation is of the order of the atom recombination 
rate. 

Because of the complexity of the chemical kinetics of 
a real gas at high temperatures, Bray utilizes Lighthill's 
model for an ideal dissociating gas* to obtain relatively 
simple mathematical expressions for the changes of 
state. A calculation has been made for air using Bray’s 
approach and the recent chemical-kinetic data for oxy- 
gen reported by Byron.” A pressure and temperature 
behind the reflected shock wave of 1,000 atm. and 
4,500°K., respectively, were selected, because at these 
conditions the relaxation effects of oxygen dissociation 
should predominate over oxygen and nitrogen vibration 
and nitrogen dissociation. For example, Logan® 
shows that, at a temperature of 4,500°K. and a density 


of one-tenth sea level density, oxygen dissociation con- 
tains 40 per cent of the total internal energy, nitrogen 
vibration about 10 per cent, nitrogen dissociation about 
2 per cent, and oxygen vibration only a negligible 
amount. The remaining internal energy is contained 
in the active degrees of freedom, which have very fast 
relaxation times. Even if the relaxation time for ni- 
trogen vibration were much slower than that for oxygen 
dissociation, the latter would predominate in the over- 
all effect on the flow because of its greater percentage of 
the internal energy. 

In this calculation two different assumptions were 
made regarding the effectiveness of nitrogen molecules 
in the recombination of the oxygen atoms; first, that 
the nitrogen was completely ineffective and served only 
as a diluent, and, second, that the nitrogen was com- 
pletely effective. These two extreme cases bracket the 
actual conditions. These calculations for the C.A.L. 
shock tunnel nozzle indicate that freezing will occur at 
nozzle area ratios of 45 and 100 for the respective as- 
sumptions regarding nitrogen. Although these area ra- 
tios are quite small compared to the overall area ratio of 
several thousand, indicating freezing early in the nozzle, 
the temperature, density, and degree of dissociation de- 
crease so rapidly that very little energy of dissociation 
is locked in the frozen state. In the calculated exam- 
ple, at an area ratio of 45, the ratios of the temperature 
and density to their stagnation values are 0.36 and 
0.0052, respectively, and the percentage of dissociation 
is 0.08. The amount of ‘‘frozen’”’ energy in dissociation 
is about 0.8 per cent of the total enthalpy, and the as- 
sumption of thermodynamic equilibrium for the test 
section calculations appears justified for these condi- 
tions. However, for much lower pressures a signifi- 
cant amount of the total enthalpy may be frozen out, 
and the assumption of thermodynamic equilibrium 
would not then be valid. 


Adequacy of Available Test Time 


In view of the brief testing times available in shock 
tubes and shock tunnels, it is appropriate to inquire 
whether or not the flow duration is adequate for aerody- 
namic investigations. Flow phenomena about a body 
may be broadly classified into three groups—the ex- 
ternal flow field, the boundary layer, and the wake. 
The establishment of quasi-steady states for each of 
these phenomena will be by wave processes, diffusion, 
and mixing, respectively. 

The external flow field generally should achieve 
quasi-steadiness within a few body lengths of flow, since 
the wave processes propagate at a sonic speed or 
greater. This has been well substantiated by schlieren 
photographs. 

The unsteady growth of a boundary layer can also be 
related to body lengths of flow, and can be expected to 
become quasi-steady about as quickly as the external 
flow field. Of course, changes in a boundary layer take 
place by molecular or turbulent diffusion, which is a 
much slower process than wave propagation. How- 


he 
ho 
ct 
eS 
as 
he : 
eS 
n- | 
nt 
yn. 
4 
ed 
es 4 
j 
he 
ne 
1g 
g- 
T, 
tO 
1e 
1e 
s, 
e 
it 
it 
n a 
a 
0 
of 
y 
: 


Fic. 7. Cornell Aeronautical Laboratory 11-in. by 15-in. Hyper- 
sonie Shock Tunnel. 


ever, except on the leeside, the boundary layer will be 
relatively thin, so that the adjustment of the boundary 
layer to flow changes is very rapid. 

The time required to establish the quasi-steady situa- 
tion in the wake is not well understood. Wave proc- 
esses are involved in establishing the overall shape, and 
the edge is a layer of high shear, as is a boundary layer. 
These two aspects of the wake suggest early flow estab- 
lishment. If the circulation within the wake is estab- 
lished by a diffusion process, a very long time would 
be required, since the wake size is of the same order as 
the body diameter. However, the pumping action re- 
sponsible for the circulation may well involve wave proc- 
esses within the wake, in which case establishment 
times would be short. This is supported by recent 
shock tube experiments*' *? which indicate that the 
wake of a hemisphere-cylinder becomes quasi-steady in 
less than 100 microseconds, or about 15 body lengths of 
flow. 

In view of these considerations, the testing time in 
the C.A.L. tailored-interface shock tunnel should be 
adequate since, as previously mentioned, it provides 
flow durations in excess of 30 body lengths. 


APPARATUS 


Description of Tunnel 


The C.A.L. tailored-interface shock tunnel consists of 
a high-pressure driver section 14 ft. long with a 3.5-in. 
inside diameter and a low-pressure driven section, also 
14 ft. long, with a 3-in. inside diameter. These sec- 
tions are designed for a working pressure of 2,000 atm. 
Located at the downstream end of the low-pressure tube 
is a convergent-divergent nozzle having a contraction 
ratio of 18.9:1, which is sufficient to ensure almost 
complete reflection of the incident shock wave. 

The divergent portion of the nozzle has a two-stage 
bilateral expansion. This arrangement permits an 
initial expansion in the horizontal plane ending in a 1.5 
by 12.25 in. cross section, and a terminal expansion in 
the vertical plane to a test section size of 15 by 11 in. 
The opening of the terminal nozzle is variable in order 
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to change the test section Mach Number, and is 
smaller than the exit of the primary nozzle in order to 
provide for a boundary-layer bleed. 

A thin secondary diaphragm of Mylar, separating the 
shock tube from the nozzle and test section, permits the 
nozzle, test section, and dump tank to be evacuated to 
very low pressures (on the order of 1 micron Hg abso- 
lute) prior to each test. This procedure is necessary 
to minimize the nozzle starting time.*® <A starting 
time from 0.5 to 1.0 millisecond is obtained in this 
nozzle. The high pressure behind the shock wave 
easily ruptures this diaphragm so that its presence 
does not affect the flow establishment. However, 
it was found that minute Mylar fragments eroded blunt 
models and resulted in flow disturbances when they re- 
bounded from the model.** 

In order to avoid this particle damage, a simple de- 
flection nozzle was inserted between the two expansion 
nozzzles. This deflector section has a two-dimensional 
wedge which spans the width of the flow and is mounted 
at a negative angle of attack. The flow passing below 
the wedge is turned upward 10° through a Prandtl- 
Meyer expansion. The terminal nozzle and test sec- 
tion are then mounted behind the wedge also at a 10° 
inclination. The entrance of this nozzle is thus shielded 
from diaphragm particles which are effectively centri- 
fuged out of the test flow, and no further erosion of the 
models has been observed. Air passing above the 
wedge, plus the unused air below, is dumped into a large 
tank. A drawing of the nozzle and test section region 
is shown in Fig. 7. 


Data Measurements and Instrumentation 


The data measured in the shock tube during each 
run are the initial pressure in the low-pressure (driven) 
tube, the driver pressure at which the primary dia- 
phragm bursts, the pressure history in the driven tube 
behind the incident and reflected shock waves, and the 
speed of the incident and reflected shocks. These data 
are used to calculate, taking into account the real gas 
effects, the flow conditions behind the reflected shock 
which are the stagnation conditions for the nozzle flow. 

Local heat-transfer rates on a model in the test sec- 
tion are determined from measurements of transient 
surface temperatures by means of a thin film resistance 
thermometer.” The theory for the heat conduction in 
the body is then used to relate these measured tempera- 
tures to the heat-transfer rate.” 

Rayleigh pitot pressure measurements in the test 
section are made using both a Statham pressure trans- 
ducer (Model PM222b * 25-200)** and a Kistler SLM 
gage (PZ-14). The test section static pressure is meas- 
ured at the sidewall using an Altec-Lansing Model 
21BR-200-2 microphone as a pressure transducer. All 
of these pressure transducers are suitably shock 
mounted for vibration isolation. 

Schlieren photographs of the flow about models are 
obtained using a multiple spark light source“ and a ro- 
tating drum camera. With this equipment, five 
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schlieren photographs can be taken during each run 
at preset time intervals varying from 50 to 1,000 mi- 
croseconds. 


EXPERIMENTAL RESULTS 


Shock Tunnel Performance 


The first experiments that were made in the tailored- 
interface shock tunnel were aimed at verifying the pre- 
dicted tailoring conditions. During these initial tests, 
in which helium was used as the driver gas, pressure 
histories were recorded at three stations along the low- 
pressure tube, and schlieren pictures were taken of the 
flow around a blunted cone using the multiple spark 
light source. The pressure measurements showed a 
constant pressure behind the reflected shock for about 
4 milliseconds, an indication of successful tailoring. 
Further evidence that the test section flow was steady 
during this period was presented by the schlieren photo- 
graphs. 

Fig. 8(a) shows the pressure records taken during one 
of these runs. The incident shock Mach Number for 


the run shown was 3.78, which is slightly less than the 
theoretical value of 4.1. Flashing of the multiple light 
source is indicated by the pips appearing on the traces. 
The 4 millisecond period of constant pressure is in good 
agreement with the predicted testing time (see Fig. 5). 
The available testing time in the test section will be 
somewhat less than this because of the nozzle starting 
phenomenon. Tests were continued over a range of 
shock Mach Numbers from 3.4 to 4.2 and the results 
show that room temperature helium successfully tai- 
lors over this range. Fig. 8(b) shows another pressure 
history for a run where the shock Mach Number was 
almost identical to the theoretical value. These results 
indicate that it is not necessary to maintain the theore- 
tical conditions precisely in order to achieve tailoring. 

Another significant point is that the attenuation of 
the shock wave is truly negligible over the last half*of 
the shock tube where the wave speed is measured. This 
can, for the most part, be attributed to the low value of 
the length-to-diameter ratio (L/D = 56) of the driven 
section.!® 

Following the experiments with helium, several tests 
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TABLE 1 
Summary of Mach Number Measurements 
Nominal Nominal Nominal Pitot Static 
Driver Test Stagnation Stagnation Static Pressure Pressure 
Gas Gas Pressure (psia) Temp. (°K.) Temp. (°K.) Mach No. Mach No. 
Helium Air 1,000—10, 000 2,000 50-75 13.4 + 0.5 12.3 = 0.5 
Hydrogen Air 5,000 4,000 220° 11 + 1 = 


were made using hydrogen at room temperature as the 
driver gas. For this case, the theoretical calculations 
predict tailoring at a shock Mach Number of 6.14, 
corresponding to a flight Mach Number of 10.6. The 
pressure histories at the 10-ft. and 13-ft. stations are 
shown in Fig. S(c). For this test the shock Mach 
Number was 6.20. The arrival of the incident and re- 
flected shock waves at the 7-ft. station appears as a pip 
on the upper trace. Although the pressure behind the 
reflected shock ultimately achieves a constant value, it 
does not attain this value as quickly as in the helium 
driver tests [see Fig. 8 (a), (b)].. A comparison of the 
records from many hydrogen driver gas tests indicates 
that this type of pressure history is present in all the 
runs. Since these tests indicate that it is not a function 
of either shock Mach Number or pressure level, this ef- 
fect does not appear to be associated with either the 
tailoring process or the interaction of the reflected 
shock wave and the boundary layer. In addition, wave 
speed measurements indicate very little attenuation. 
Combustion at the hydrogen-air interface is a possible 
explanation of this phenomenon. To check this possi- 
bility, several tests were made using nitrogen as the 
driven gas and a hydrogen driver, since this combina- 
tion would tailor at nearly the same conditions as hy- 
drogen and air. These tests exhibit a pressure history 
behind the reflected shock similar to that of the helium- 
air tests [Fig. 8 (a), (b)], adding additional weight to the 
interface-combustion explanation. 

Having established the tailored-interface modifica- 
tion of the shock tunnel as a practical technique, a ser- 
ies of experiments was made to measure the test section 
Mach Number. Using the pressure transducers previ- 
ously described, the Rayleigh pitot and the wall static 
pressure were measured in the test section. Mach 
Numbers have been computed from the ratios of these 
pressures to the pressure behind the reflected shock 
wave, assuming that the flow is isentropic. The re- 
sulting Mach Numbers are shown in Table 1. The 
plus and minus variations represent the scatter from 
run to run as well as the variation with stagnation 
pressure. For the helium tests, the two methods of de- 
termining Mach Number give values that consistently 
differ by about one Mach Number. It is not known at 


this time whether this is due to instrumentation error or 
a stagnation pressure loss from either viscous effects in 
the primary nozzle or wave phenomena in the deflec- 
tion nozzle. The lower Mach Number measured dur- 
ing the hydrogen tests results from increased real gas 
effects at the higher stagnation temperature. 

The effective nozzle area ratios, inferred from the 
pitot pressure Mach Numbers, are about 3,000 and 
2,500 for the helium and hydrogen tests, respectively, 
compared to the geometrical area ratio of 5,000. This 
reduction, which is attributed to boundary-layer 
growth, is comparable to that of other hypersonic 
tunnels. 


Aerodynamic Experiments 


The first aerodynamic tests made in the tailored-in- 
terface shock tunnel were a series of experiments with 
both helium and hydrogen drivers to measure heat- 
transfer rates on a hemisphere-cylinder model. The 
temperature histories at several angular positions are 
shown in Fig. 9(a), and the variation of the stagnation- 
point heat-transfer rate with time is shown in Fig. 9(b) 
for a typical run. The stagnation heat-transfer rate is 
relatively constant from 0.8 to 4.0 milliseconds at an 
average value of 52 B.t.u./ft.*-sec. In calculating the 
heat-transfer rate from the temperature histories, the 
temperature records are read at 100 microsecond inter- 
vals. This accounts for the discrete points shown in 
Fig. 9(b). Although the maximum variation in the 
heat transfer is approximately +18 per cent of the 
average value, SO per cent of the data points, in this 
time interval, lie within +10 per cent of the average. 
This scatter includes the experimental errors as well as 
the effects of flow fluctuations and electrical noise. 

A comparison of the experimental stagnation-point 
heat-transfer rates with theoretical predictions of Fay 
and Riddell* has yielded additional information re- 
garding the test section Mach Number. Typical re- 
sults are summarized in Table 2. 

The experimental stagnation-point heat-transfer rates 
for hydrogen/air tests agree very well with the values 
predicted theoretically for a Mach Number of 11.2. 
This Mach Number value agrees with that determined 
from the pitot pressure measurements with the hydro- 


TABLE 2 
Comparison of Experimental and Theoretical Stagnation-Point Heat-Transfer Rates 


Stag. Stag. 

Driver Test Shock Pres. Temp. 
Gas Gas Mach. No. (psia) (°K.) 
Air 3.85 2,860 1,880 

Air 3.93 2,970 1,930 

Air 6.55 3,000 4,140 

Air 6.45 3,640 4,050 


Exper. Heat Theoretical Heat Transfer 
Transfer 
(B.t.u./ft.*-sec. ) M,. = 13.4 = 12.0 = 11.2 
52.0 39.8 51.3 
53.6 41.2 53.3 
128.7 125.7 
157.5 154 2 
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gen driver. In the case of the helium driver gas tests, 
however, the experimental data agree with the theory 
at a Mach Number of 12, which is consistent with the 
static pressure measurements rather than the pitot 
pressure measurements. It should be noted that, when 
the differences in the test conditions between runs are 
accounted for in the theoretical heat-transfer calcula- 
tions, the differences in the experimental data from run 
to run are satisfactorily resolved. 

A comparison of the experimental heat-transfer dis- 
tribution over the nose of the model with Lees’ theoret- 
ical prediction” for a laminar boundary layer is shown 
in Fig. 10. The data include tests in air using both he- 
lium and hydrogen drivers as previously tabulated plus 
one test in nitrogen using a hydrogen driver. The theo- 
retical curve is for y = l.4and J, = 10. The agree- 
ment with the theory is quite good for all the tests, 
which include a range of stagnation temperatures from 
1,700°K. to 4,100°K. in air and a value of 6,300°K. in 
nitrogen. The variation of the data during a typical 
run is indicated by the vertical bars alongside the data 
points. 

Receut research in the C.A.L. tailored-interface shock 
tunnel has included a study of the hypersonic shock- 
boundary-layer interaction on a flat plate and also of 
the heat transfer on a slender yawed cone.” It should 
be noted that it is in research of this type, where high 
Mach Number and high stagnation temperature condi- 
tions are required, that the shock tunnel is most valu- 
able. 


CONCLUDING REMARKS 


The tailored-interface shock tunnel offers a means 
whereby the high Mach Number and high stagnation 
temperature flows necessary to simulate hypersonic 
flight conditions can be easily produced. An impor- 
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tant feature of this type of shock tunnel is that it pro- 
vides a longer testing time than a conventional shock 
tunnel and, hence, a means of alleviating flow unsteadi- 
ness in the test section caused by shock-wave attenua- 
tion. In addition, there appear to be no unusually 
difficult problems related to its design and operation. 
In its present state of development as a hypersonic 
test facility, the shock tunnel is particularly well suited 
to heat-transfer studies. Recent advances in force and 
pressure measuring techniques, however, promise to ex- 
tend its application to a wider range of investigations. 
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The Smallest Height of Roughness Capable of 


Affecting Boundary-Laver Transition 


A. M. O. SMITH* ann DARWIN W. CLUTTER** 
Douglas Aircraft Company, Ine. 


SUMMARY 


An investigation was made to determine the smallest size of 
isolated roughness that will affect transition in a laminar-bound- 
ary layer. Critical heights for three types of roughness were 
found in a low-speed wind tunnel. The types were (1) two- 
dimensional spanwise wires, (2) three-dimensional discs, and (3) 
a sandpaper type. In addition to type of roughness, test vari- 
ables included the location of roughness, pressure distribution, 
degree of tunnel turbulence, and length of natural laminar flow. 

The most satisfactory correlation parameter was found to be 
the roughness Reynolds Number, based on the height of rough- 
ness and flow properties at this height. The value of this critical 
Reynolds Number was found to be substantially independent of 
all test variables except the shape of roughness. This parameter 
also correlates well other published data on critical roughness in 
low-speed flow. The value of the roughness Reynolds Number 
necessary to move transition forward to the roughness itself was 
also determined for the three types of roughness and was found 
to be approximately constant for a given type of roughness. 

An investigation of the limited amount of available data on 
critical roughness in supersonic flow indicates that the effects of 
roughness may still be correlated by the roughness Reynolds 
Number. The value of this Reynolds Number depends primarily 
on the Mach Number at the top of the roughness. When this 
Mach Number is greater than 1.0, the roughness Reynolds Num- 
ber based on conditions behind a shock is probably the charac- 
teristic parameter. 


SYMBOLS 
A = constant 
= speed of sound 

c = chord of airfoil or test surface 

d = diameter or width of roughness element 

H = boundary-layer shape parameter 6*/6 

k = roughness height 

= Reynolds Number, unless otherwise noted 

R, = chord Reynolds Number U’,,c/v 

R, = roughness Reynolds Number u;,k/v 

R5* = Reynolds Number based on displacement thickness 
U6*/v 

Rg = Reynolds Number based on momentum thickness U’@/y 

i’ = local velocity outside boundary layer 

u = velocity within boundary layer in x direction 

ve = friction velocity V ty p 

x = distance along surface measured from forward stagna- 
tion point 

y = distance perpendicular to surface 

a = angle of attack 

8 = shape parameter of Hartree boundary-layer profile, 
U = 

6 = boundary-layer thickness 

6* = boundary-layer displacement thickness 

6 = boundary-layer momentum thickness 
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X = length of a wave in the surface 

“a = dynamic viscosity 

v = kinematic viscosity 

p = mass density 

Subscripts 

crit = evaluated at the critical roughness 

H.W. = evaluated at the hot wire 

k = evaluated at the roughness element 

i = evaluated at transition 

w = evaluated at the wall 

0 = natural transition 

oo = conditions in undisturbed flow 

6 = evaluated at edge of boundary layer 

1 = evaluated ahead of normal shock 

2 = evaluated aft of normal shock 
INTRODUCTION 


De PUBLISHED DATA on the effects of isolated 
roughness elements on transition fall primarily 
into two categories: (a) the effect of a roughness ele- 
ment on the transition point as a continuous function 
of all roughness heights, and (b) the more restricted 
problem of determining the largest height of roughness 
that does not affect the transition point. Category (b) 
is the principal subject of this paper. 

Although the published data on critical roughress 
heights cover a large range of critical roughness R: y- 
nolds Nuinbers (Rx crit = crit = 60 to 1,100), 
the data do not cover the range of relatively high veloc- 
ity and small critical height. The test surfaces used in 
most of the investigations had rather low transition 
Reynolds Numbers. Moreover, the various investiga- 
tors have presented their results in a variety of ways by 
means of a variety of correlation parameters, such as 
k/5*, Rsx, R,., etc. Most of the data were obtained in 
low-speed flow. 

The primary part [Part (I)] of this paper concerns 
itself with low-speed flow and is essentially a digest 
of reference 1. Part (II) extends the results to high- 
speed flows. 


PART (I) CRITICAL ROUGHNESS IN 
LOW-SPEED FLOW 


The principal purpose of the present investigation was 
to deter.nine the smallest or critical roughness height 
over a large range of Reynolds Numbers, as affected by 
(1) body pressure distribution, (2) roughness shape, 
(3) roughness height, (4) roughness location, (5) tandem 


¢ This research was sponsored by the U.S. Navy Bureau of 
Aeronautics, under Contract NOa(s) 54-773c. 
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arrangement of roughness elements, and (6) wind- 
tunnel turbulence. 

Of course, an immediate purpose was to establish, if 
possible, a general method for correlating the critical 
roughness heights with parameters of the flow (such as 
Re or R,). Fulfillment of this objective was achieved 
partially by a general examination of already available 
test data and partially by creation of new data, as de- 
scribed later. Because of the variety of presentation 
made in the past, considerable effort was expended to 
bring most of the data, including our own, together into 
a consistent presentation. 


Factors AFFECTING TRANSITION DUE TO ROUGHNESS 


Transition from laminar to turbulent flow is brought 
about by various types of disturbances in the laminar 
boundary layer. On flat or convex surfaces the magni- 
tude or growth of disturbances is dependent on one or 
more of the following: 


(a) Pressure gradient 

(b) Free-stream velocity and viscosity 

(c) Free-stream turbulence 

(d) Internal flow disturbances, such as blowing into 
the boundary layer, thermal gradients, and vi- 
brations 

(e) Secondary-flow effects, as, for example, the ef- 
fects of sweep 

(f) Surface roughness 


In this investigation we are primarily interested in 
disturbances introduced by the last-named factor— 
that is, roughness—but whether or not these disturb- 
ances affect transition will depend on both the static 
and the dynamic properties of the boundary layer, 
which are dependent on the other factors listed above. 

Neglecting free-stream turbulence and pressure gradi- 
ent, the transition point x7 due to a given shape of 
roughness is a function of the following variables: 


k = height of roughness 

Xz = position of roughness 

Ux = velocity in the undisturbed boundary 
layer corresponding to the top of the 
roughness 

U = velocity outside the boundary layer 

(Ou/Oy), = velocity gradient at top of roughness 

v = kinematic viscosity 


The location of transition x7 can be written functionally 
as 


xp = f(R, xx, Ux, U, Ou/Oy,, v) (1) 


Obviously these variables are not independent of each 
other, for, in a particular flow, 


= f(k, xx) 
(Ou/Oy), = f(k, Xx) 
Therefore, Eq. (1) can be reduced to 
xr = f(k, xx, U, v) (1a) 


and also 
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These variables may be grouped into various nondimen- 
sional parameters. One such parameter is the rough- 
ness Reynolds Number, 


R, = (2) 


If k is not large compared to the boundary-layer 
thickness, we know from the approximately linear 
nature of the bottom of the boundary-layer profile that 


u,/U = A k 0. (3) 


where A is a constant. 
The roughness Reynolds Number may now be written 


Ry = (uk/v) = = 
A (A) 


Eq. (4) is a functional relation of three dimensionless 
parameters. 

Eq. (4) is valid for small k only; if k is so great as to 
be outside the boundary layer, 


u,/U=1 
and = (k/0)Ro, (5) 


That is, R, varies linearly with k/6, for large k, whereas 
for small k it varies quadratically as k/6,. 

Of course, there are other nondimensional parameters 
that could be obtained from the variables above; a 
relation involving R;, was chosen because a Reynolds 
Number is inferred as the critical parameter if it is as- 
suined either that a wake is created or that vortices are 
shed from the roughness element at a certain critical 
Reynolds Number and that one or both of these dis- 
turbances either cause immediate transition or, when 
amplified, cause transition farther downstream. The 
criterion that R, ¢;;; = constant was first proposed by 
Schiller in 1932 (see reference 2, page 311). 


PREVIOUS INVESTIGATIONS 


Most previous investigators of the effect of roughness 
on transition have presented their results in the form of 
parameters similar to those of Eq. (4) above (see refer- 
ences 3 through 14). The parameters are similar in 
that they include a roughness Reynolds Number based 
either on the roughness height or on the boundary- 
layer thickness at the roughness and a roughness- 
height parameter such as k/6, or k/6,*. The investiga- 
tions indicate that the shape of the roughness has a 
major effect on the values of the critical-roughness pa- 
rameters, but, on the other hand, they indicate that 
surface curvature, history of the boundary layer, and 
free-stream turbulence have only minor effects on the 
values, even though these latter factors are major ones 
in establishing natural transition on a smooth surface. 

The various investigations show disagreements on 
several points: (1) the effect of roughness shape on the 
critical roughness height, (2) the ratio of critical height 
to the height of the boundary layer at the roughness, 
and (3) whether or not the length of natural laminar 
flow greatly affects the value of critical roughness. 
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Although the previous investigations covered a wide 
range of types, heights, and positions of roughness, they 
did not give information on the case of a small roughness 
element close to the leading edge. This is perhaps the 
most important case in determining the tolerable rough- 
ness on a body. Part of the purpose of the investiga- 
tion reported here is to fill this gap and reconcile the 
discrepancies in the existing data. 


TEST FACILITIES AND TECHNIQUES 


Wind Tunnel 


All tests in the present investigation were made in the 
Douglas El Segundo 30- by 45-in. wind tunnel. Air 
speeds as high as 185 ft. per sec. were used during the 
tests. The settling chamber at the entrance of the con- 
traction cone of the tunnel accommodates six removable 
screens, 40 mesh by 0.0065-in. wire, spaced 10 in. apart. 
For the majority of tests, all six screens were used, but, 
for the purpose of learning the effect of turbulence, a 
few tests were run with one screen and some with no 
screens. With all screens in place the turbulence is 
quite low; transition measurements on the model, 
which are described below, gave an amplification factor 
of e*-7, which compares favorably with values obtained 
in other low-turbulence tunnels and in flight.’ 


The 108-In.-Chord Test Surface 


The model was designed as a nearly constant-veloc- 
ity, two-dimensional surface spanning the tunnel from 
top to bottom. The model profile began as a sharp 
30° wedge, fairing into a flat-sided surface 25 in. aft of 
the leading edge, being so designed as to give a constant 
velocity over the curved portiou.'® The model is 
shown in Figs. 1 and 2. The flat section was 2.254 
in. thick and extended aft from 25 in. to 96 in. from the 
leading edge. The last 12 in. of the chord was a simple, 
radius-nose, straight-sided flap, which could be deflected 
between —13° and +13°. 

Surface joints could not be made to fair as smoothly 
as was desired on the test surface; therefore the entire 
test surface was covered with a single sheet of “‘Mylar’’ 
polyester film (manufactured by the Dupont Com- 
pany). The ‘“Mylar’’ film was 0.005 in. thick and had 
an average root mean square roughness of 2.5 and a 
maximum of 4 microinches. The film was folded 
around the nose and carried aft to the leading edge of 
the flap, being secured to the model surface with strips 
of pressure-sensitive tape (adhesive on both sides). 
This ‘“‘Mylar’’ film proved satisfactory for fairing out 
the discontinuities of the metal surface, and the film 
stood up well for the duration of the tests. 

Since the test section of the wind tunnel is only 4 ft. 
in length, the leading edge of the model extended 4.5 in. 
into the contraction cone of the tunnel and extended 
aft 55.5 in. into the diffuser. A profile of the tunnel 
wall is shown in Fig. 1, along with the model. 


Types of Roughness 


The roughness tested was primarily of three types: 
spanwise round wires, 1/16-in. diameter circular discs, 


and sandpaper. A few tests were made with the Mylar 
sanded at right angles to the direction of the airflow, 
but these will not be discussed in this paper; for further 
details, see reference 1. The wires may be thought of as 
representing spanwise skin joints; they are, of course, 
two-dimensional in nature. The discs may be consid- 
ered as representing spots of dirt or incompletely 
countersunk flush rivets; they are a three-dimensional 
type of roughness. The sandpaper and scratched 
Mylar cannot very well be classified dimensionally. 
One or the other or both might represent roughness 
such as may arise from raindrop and dirt erosion, dust 
on the surface, mud spots, or the like. 

Stretch-straightened lengths of brass or copper wires 
were stretched tightly spanwise across the test surface 
at the selected location. They were secured to the sur- 
face by taping the ends and by nitrate dope applied to 
the full span of the wire. The wires used in the tests 
varied from 0.002 in. to 0.043 in. in diameter. Tests 
were not confined to a single wire since, in addition, 
several runs were made with wires in tandem—that is, 
one located a few inches behind the other. 

The second type of roughness tested consisted of 
spanwise rows of brasscylindrical discs of 1 / 16-in. diame- 
ter. The discs ranged in height from 0.0029 in. to 
0.030 in. The spanwise spacing of the discs was 1.25 
in. throughout, except for a few runs in which only a 
single disc cemented to the center was tested. 
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Fic. 1. Details of the two-dimensional test surface. 


Fic. 2. Typical china-clay transition pattern on the two- 
dimensional test surface. The irregularity of the transition locus 
is due partly to transverse contamination from the tunnel walls 
and partly to irregularities in the roughness element. The 
rearmost point was always used in reduction of data. 
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Fic. 3. Variation of velocity and boundary-layer parameters 
with chordwise position for a = 0°. 


Two materials were used to represent the sandpaper 
type of roughness—No. 100 grit Crystolon, an abrasive, 
and ground steel shot, a steel grit used in sandblasting. 
To apply the grit, the Mylar surface of the model was 
masked off locally to leave a 1/4-in.-wide spanwise gap. 
Nitrate dope was applied over the exposed strip and the 
grit was blown onto the wet dope. Before the dope had 
hardened, the masking tape was removed so that no 
rough edge was left in the dope by the tape. In this 
way satisfactory | /4-in.-wide strips of roughness were 
made. 


Measurement of Roughness Height 

Heights of wire, disc, and sandpaper roughness were 
always measured on the model by means of a roughness 
gage consisting primarily of an Ames dial gage, reading 
to 0.0001 in. The measurement sought and used was 
maximum height of roughness along the entire span. 
For the wires, every inch of each wire was checked. For 
the rows of discs, the height of each individual disc was 
measured. For the sandpaper type of roughness, the 
entire area of roughness was searched to find the highest 
point. When measured this way the actual roughness 
heights were found to be appreciably greater than the 
nominal heights. 

The maximum error in determining the roughness 
height was 7 per cent for the wires and discs, with the 
exception of the 0.002-in. wire, where the error was 
about 15 per cent. For the sandpaper type of rough- 
ness, the error was less than 12 per cent. 


Methods of Determining Transition 

Two methods were principally used to determine 
transition—the china-clay and the stethoscope tech- 
niques. These two techniques were checked with the 
hot-wire anemometer and surface pitot tube technique. 
Usually the various techniques indicated values of Rr 
that agreed with each other within 4 per cent. 


RESULTS 


Transition Measurements With and Without Roughness 


The tests were made with the model at angles of at- 
tack of +1°, 0°, —1°, and —2°. The velocity distri- 
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bution for the model at an angle of attack of 0° is 
shown in Fig. 3. The boundary-layer parameters 
needed in the data reduction are also shown. These 
parameters were determined by the method presented 
in reference 17. 

Fig. 4(a) shows the variation of the natural transition 
Reynolds Number with chord Reynolds Number for 
four angles of attack. The highest transition Reynolds 
Number Ry = 7.4 X 10° was obtained at an angle of 
attack of —2°. Fig. 4(b) shows the corresponding 
variation of R»y,, the Reynolds Number based on the 
momentum thickness at transition, w.th chord Reyn- 
olds Nuimber. 

Fig. 5 illustrates the effects on transition produced by 
the various types of roughness having the same height. 
It shows that, as the air speed is increased beyond the 
point where the roughness first begins to cause prema- 
ture transition, the rate of forward movement of transi- 
tion with speed increases rapidly until transition occurs 
at the element itself. The rate of this movement for- 
ward is greater for the three-dimensional disturbances 
than for the two-dimensional. <A single disc has no 
effect on natural transition until it produces a wedge of 
turbulence whose apex is at the disc itself. Fig. 5 also 
exhibits the fact that three-dimensional roughnesses 
are less apt than two-dimensional to promote transition. 


Transition Measurements With Several Values of Tunnel 

Turbulence 

The effect of varying the degree of tunnel turbulence 
on the natural transition Reynolds Number is shown in 
Fig. 6. The Figure shows the decrease in transition 
Reynolds Nuinber as the number of screens in the tun- 
nel was reduced from six to one or none. The tunnel 
turbulence was not measured in this series of tests, but 
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Fic. 4. Variation of the natural transition Reynolds Number 
and the Reynolds Number based on momentum thickness at 
transition with chord Reynolds Number for various angles of 
attack. 
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previous measurements and calculations indicated the 


following: for no screens, 


100V (u’? + v’? + w’) /3/U & 0.26 to 0.30 


for Six screens, 


100V (u’? + v2 + w’2)/3/U 0.03 


where wu’, v’, and w’ are root mean square values of the 
u, v, and w fluctuations, respectively. 

Fig. 7 shows the effect of increased tunnel turbulence 
on the heights of roughness required to cause premature 
transition. There is a very slight increase in critical 


roughness with the increased tunnel turbulence. 


Definition of Critical Roughness 
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Fic. 6. Variation of natural transition Reynolds Number 


Observation of plots of (x/c)r versus R, shows that it 
is frequently difficult to determine exactly when a given 
roughness becomes effective in reducing the length of 
natural laminar flow. Therefore, a new definition has 
been used in the data reduction, and this definition is 
presented in this Section. 

Plots can be made of the ratio of the transition Reyn- 
olds Numbers with and without roughness Rr/Rr, 
versus the height parameter k/6,. Typical curves of 
such plots are shown in Fig. 8. As can be seen from 
these plots, it is difficult to determine exactly the value 
of k/@, at which the roughness becomes effective in 
reducing the natural transition Reynolds Number— 
that is, when the ratio of transition Reynolds Numbers 
Rr / Rr, becomes just less than unity. Therefore, the defi- 
nition of critical roughness for a particular configuration 
was taken as the value of k for which Rr has decreased to 
0.95 Rr, Using this definition, the critical value of 
k/0, can be determined from the plots of k/@, versus 
R,/ Rr, for a particular roughness height and location, 
as is shown in Fig. 8. This critical value of k/@, can 
then be used to determine the critical value of Ro,, 


R,,, ete. The same type of plot—that is, Rr/Rr, 
NO ROUGHNESS 
SANDPAPER TYPE 
ee ROW OF DISCS 
SINGLE DISC 
(x) 
Hz),= (2), 
Re 
Fic. 5. Variation of transition point with chord Reynolds Num- 


ber for various types of roughness having the same height 


with chord Reynolds Number with six, one, or no turbulence 
screens installed in the tunnel. 


versus k/@#,—was used to determine the roughness 
values when transition reached the roughness element 
itself. 

The value of the critical roughness Reynolds Number 
can be determined from the relation 


Ry, = = (R/Ox) (Ux 


where all variables are determined at critical conditions. 
The value of u,/ U can be determined from the known 
model velocity distribution (8 from L’ CxP/2-)) 
and the tables of u/U’ of reference 18. 


Method of Presentation of Critical-Roughness Parameters 


The results of this investigation will be presented on 
plots similar to that shown in Fig. 9, which shows the 
relation between Ry, Ro,, ux/U, and k/6, for several 
boundary-layer profiles. Fig. 9 shows the relation be- 
tween Ry, k/6;, and R, for R, = 200 and various pres- 
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Fic. 7. Variation of transition point with chord Reynolds 
Number for single spanwise wires with six, one, or no turbulence 
screens installed in the tunnel, ata = —2°. 
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LINE DETERMINING CRITICAL ROUGHNESS 
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Fic.8. Variation of ratio of transition Reynolds Number, with 
and without roughness Rr/R7z,, with k/6 for single spanwise 
wires at a = 0°. 
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Fic. 10. Critical values of k/@, and Ro ~!/? for single spanwise 
wires. Open symbols are for a = 0° and —2°. Solid symbols 
are fora = —1° and +1°. Half-solid symbols denote values at 
increased tunnel turbulence. Lines of Ry = constant represent 
Blasius boundary layer (8 = 0). 
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Fic. 12. Critical values of k/6,; and Ra ~'/? for 1/4-in. strips 
of sandpaper type of roughness. Values from reference 14 are 
also shown. Lines of R; = constant represent Blasius boundary 
layer (8 = 0). 
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Fic. 9. The relation between Ry, Ra, uz/U, and k/6 for several 
boundary-layer profiles. 
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Fic. 11. Critical values of k/@, and Ré,~'/? for spanwise rows 
of discs. Lines of Ry = constant represent Blasius boundary 
layer (6 = 0). 
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Fic. 13. Critical values of k/@, and Re,~'/? for spanwise 
wires of this and previous investigations. Lines of Ry, = con- 
stant represent Blasius boundary layer (8 = 0). 
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sure distributions (8 = 0, 0.1, 0.5, and 1.0 from the 
Falkner-Skan family of velocity profiles). Over the 
range of 6’s shown, the shape of the boundary-layer 
profile has only a minor effect on the k/# — i. -- 
R, relation. For the present investigation, 8 varied 
from 0.025 to 0.64 in the test region. Natural transi- 
tion sets a lower limit for the values of critical Rs, 
and this lower limit may be estimated from stability 
theory and experimental results." The lower limit is 
shown in Fig. 9, for the curves of R, = 200 and the 
various 6’s and a small-perturbation amplification fac- 
tor of e*; that is, /Bidt = 9 (reference 15). It should 
be noted that this lower limit of Ry, is a function of 8 
only. Also shown in Fig. 9 is the relation between 
Ro, k/@, and R, for R,’s equal to other constant 
values for 8 = 0; that is, flat-plate flow. The relation 
between u;,/U and k/6, for 8 = 0 is also presented. 
The presentation of the data on a three-parameter 
plot with R,, Re,~'?, and k/6, as the variables is con- 
sistent with the relation for R, developed above, 


R, A(k/0x)?Ro, (k/6,, small) (4) 


Critical-Roughness Data 


Fig. 10 shows the values of critical roughness of wires 
of various sizes and locations at various angles of at- 
tack. The critical roughness points tend to follow the 
lines of R, = 100 with a scatter of approximately 50 
per cent, except at the most forward location of rough- 
ness—that is, at x, = 0.25 in.—where the values of R, 
were as high as 260. The values of critical roughness 
appear to be independent of angle of attack, with the 
possible exception that, for those wires located close 
to the leading edge (x, = 0.25 in.), there is a trend for 
R, and k/6; to decrease as a varies from 0 to —2°. 

The increase of R, as k approaches the leading edge is 
consistent with the observations presented by Preston 
in reference 19. Preston points out that there is a lower 
limit to the Reynolds number Ry, at which fully de- 
veloped turbulent flow can occur. This limit can not 
only be predicted by theory but is also observed experi- 
mentally. The lower experimental limit is about Ry = 
320. To cause transition at a point where Rg < 320, 
the device producing transition must add a drag to the 
boundary layer to raise Ry to a value of 320. Clearly, 
any roughness has a drag and, consequently, increases 
the momentuin thickness of the boundary layer. Of 
course, in addition to the criterion that R» > 320, dis- 
turbances must be present and must be amplified, to 
establish turbulent flow. The amplification of the dis- 
turbance depends on both Ro and the pressure gradient. 
In the present setup, as one approaches the leading edge. 
Rs becomes small and the pressure gradient becomes 
more favorable. Both of these factors have a stabiliz- 
ing effect on any disturbance in the boundary layer and, 
consequently, the roughness must produce larger and 
larger disturbances as it approaches the leading edge if 
it is to produce‘early transition. 

In any case, the variation of R,. .,i with pressure dis- 


tribution is small. The effect of tunnel turbulence on 
the values of critical roughness of the wire type is indi- 
cated in Fig. 10. There appears to be a slight increase 
in the values of Ry .-, and (k/6,)-ri with increased tun- 
nel turbulence. 

Fig. 11 gives the critical values of roughness for the 
spanwise rows of discs and for single discs. The Figure 
shows more scatter and higher critical values of R, and 
k/6, than does the similar plot for the two-dimensional 
roughness. The scatter may be due to the effect of 
shape, which subject is discussed further below. The 
single discs show higher critical values of R, than do 
rows of discs. This fact is not surprising, for the dis- 
turbance produced by a disc spreads in a wedge-shaped 
wake behind the disc and the interaction of several of 
these disturbance wedges produces earlier transition 
than that caused by a single disc. It is worth noting 
again that a single disc had no effect on the natural 
transition point until transition was produced at the dise 
itself. 

The critical values for the sandpaper type of rough- 
ness are presented in Fig. 12. In the Figure the two 
symbols with a line joining them are for the same rough- 
ness configuration for which it was impossible to deter- 
mine exactly the critical values—thus, the range. The 
critical values shown were calculated from conditions 
at the leading edge of the 1/4-in. strip of roughness. 
They generally fall in the range of R, from 250 to 300. 
Fig. 12 shows that the span of the strip of the sandpaper 
type of roughness has little effect on the values of critical 
roughness. 


COMPARISON AND DISCUSSION OF THE RESULTS OF THIS 
AND PREVIOUS INVESTIGATIONS ON CRITICAL 
ROUGHNESS 


Correlation of Data With Various Parameters 


A study of various parameters for correlating the 
critical roughness values, including those parameters 
used by the previous investigators, resulted in the 
choice of R, ¢-i: as the most suitable parameter. In 
some instances, the parameters used previously give 
better correlation for particular conditions, such as a 
particular size, location, or shape of roughness, but 
have one or more of the following disadvantages: (a) 
they are valid over only a small range of Reynolds 
Nuinbers, (b) they are restricted to particular regions 
on a body, or (c) they depend on the length of natural 
laminar flow. 

One of the parameters used previously will be com- 
pared with R, for correlating the results of this investi- 
gation. Several investigators select k/5* or k/6 as one 
parameter for correlating the effect of roughness on 
transition. If the parameter is really fundamental, 
k/6,* or k/6, should be a constant for critical roughness 
of a given shape. If k/6,* or k/é, is constant, it fol- 
lows that &,,;, must approach zero at the leading edge 
of a flat plate; otherwise k/6,* would become infinite. 
Experience indicates that the permissible roughness in 
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this example is something greater than zero. Hence 
k/5,* is not fundamentally a sound parameter. 

In the present investigation @, is chosen in preference 
to 6,* or 6g as the parameter to represent the boundary- 
layer thickness, because it is usually the most significant 
of the three common thickness parameters. The varia- 
tion in the shape parameter /7 = 6*/@ is relatively small 
over the test region (//, varies from 2.26 to 2.56). 
Thus, a plot of critical roughness in k/@;, m,* 2 co- 
ordinates would look like one in k/6,*, Rs*,~'!? coordi- 
nates. Dryden,’ using the results of several investiga- 
tions on two-dimensional roughness, shows that a plot 
of Rr/Rr, versus k/6,* gives a single curve for the region 
where the roughness moves transition from its natural 
location x7, to the location of roughness itself x,. A 
similar plot of the results of this investigation shows no 
such correlation. Fig. 8 is typical of the Rr/Rr* 
versus k/6, plots of this investigation. The plot shows 
a wide variation of Rr/Rr, with k/6,, the value of k /6,. 
increasing as the roughness moves forward on the test 
surface. One explanation of this apparent discrepancy 
is offered. Dryden does not give the location of the 
roughness elements for all the data he presents, but a 
check of three of the four sets of data he used showed 
that the locations of the roughnesses vary from 20 to 
67 in. from the leading edge. The present investiga- 
tion dealt mostly with locations 16 in. or less from the 
leading edge. The curves of R7/Rr, versus k/6, tended 
to move closer together as the roughness was placed 
farther aft on the test surface. Moving the roughness 
even farther back would undoubtedly have restricted 
the variation of k/@, still more. This is physically 
plausible, for as x, increases, the distance between the 
roughness location and natural transition (between 
x, and x7,) decreases, and the range of 6,, and thus of 
k/0,, for a given k correspondingly decreases for the 
movement of transition from x7, to x,. Inspection of 
Figs. 9 and 10 shows that, as x, increases, the critical 
parameters k/@, and Ro," * both decrease, and even- 
tually Re, would approach the limiting value given by 
stability theory for a ‘‘clean’’ surface. 

Other investigations—that of Klebanoff et al.‘ on 
three-dimensional roughness and that of Brinich”’ on a 
wire type of roughness in supersonic flow—whose re- 
sults are presented after the method of Dryden, show 
that moving the roughness position forward increases 
the values of k/6,*. The roughness positions for both 
of these investigations are closer to the leading edge 
than those presented by Dryden. 


Comparison of Critical Roughness Data 


Before comparing the values of critical roughness 
from previous investigations and those of our investiga- 
tion, the definition of critical roughness will first be 
discussed. The previous investigators used one or more 
of the following techniques to define critical roughness: 
a hot wire placed just forward of the clean-airfoil transi- 
tion point, the china-clay technique, or observation of 
the drag rise on an airfoil. With the data of this inves- 
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tigation, several calculations have been made for which 
critical roughness is defined as the height that reduces 
the transition Reynolds Number to 99 per cent of its 
natural value, that is, Rr/Rr, = 0.99 instead of Rr/ Rp, 
= 0.95, as in the definition given earlier. A comparison 
of the critical values of roughness as determined by the 
two definitions shows that the maximum difference in 
k/6, determined by the two definitions is less than § 
per cent for wire type roughness and even less for three- 
dimensional elements, since for them the transition 
point moves forward more abruptly once the roughness 
has become critical. 

The critical-roughness values observed for spanwise 
wires in this investigation and in the investigations of 
Cowled and Fage are presented for Comparison in Fig. 
13. The several sets of data show similar variation of 
critical k/6, with R»,~'/*; for example, the sharp rise in 
(k/O).,;, and R, .,;, for the wires nearer the leading edge. 
The amount of scatter of the various results is about the 
same. The values of R, ,,,, lie between 40 and 260, 
with a mean of about 125. 

The critical values for wires in tandem are shown in 
Fig. 13 along with the single-wire data. The values 
shown are based on the wire that gave maximum value 
of Ry crit The results help to verify Cowled’s discovery 
that transition was dependent only on the most critical 
roughness when two wires were in tandem. 

A study of the various investigations of three- 
dimensional roughness has shown greater values of 
R; ¢rit than those for the two-dimensional roughness, 
and also a great deal more scatter. The study has 
shown that some of the scatter is due to the shape of the 
roughness. Fig. 14 shows the variation of R, .,i: with 
k/d for several investigations. The diameter of the 
roughness elements is d. There is a tendency for the 
values of R, .;i, to increase with an increase in k/d. 
The one exception is the work of Gregory and Walker, 
where single roughness elements rather than spanwise 
rows of elements were tested. In their tests on single 
roughness elements, R;, .;:: increases with an increase in 
k or a decrease in x, for a constant k/d ratio. It is 
also worth noting that as the value of k/d decreases, 
the three-dimensional type of roughness approaches a 
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Fic. 14. Critical values of Rx and shape parameter k/d for 
various three-dimensional types of roughness of this and previous 
investigations. A limiting line of minimum Rx ¢;i¢ is shown. 
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two-dimensional type and, thus, the lower values of 
crit follow. 

No attempt has been made to determine the effect 
of spanwise spacing of the three-dimensional elements 
on the values of R, «ri. However, as was noted above, 
single discs gave higher values of R,. .,;: than a row of 
such discs. Klebanoff et al.‘ present one example 
showing that a decrease in the spacing between elements 
gives a slightly lower value of k/6,* for a given Rr. It 
is reasonable to expect that decreasing the spacing be- 
twveen elements would lower the values of be- 
cause the row becomes more nearly two-dimensional in 
nature. 

The critical roughness values observed for the sand- 
paper type of roughness in the present investigation and 
in an investigation by the NACA" are shown in Fig. 
12. Both investigations were made on 1/4-in.-wide 
strips of roughness. The NACA investigation shows 
higher values of critical roughness than those obtained 
in the present investigation. One possible explanation 
of this difference is that the density of the grains of sand 
was greater in this investigation, and thus the strips of 
roughness acted more like a two-dimensional roughness 
than those used by the NACA. 

One type of surface irregularity of major importance 
that was not investigated in the present wind-tunnel 
tests is surface waviness. Cowled!! presents the re- 
sults of such investigations. His report shows that 
waves give much higher critical values of R, than the 
other types of roughness, values falling in the range of 
Ry crit from 1,000 to 4,200. There is no consistent cor- 
relation shown on the k/6, versus Re,~'/* type of plot. 
An approximate single-curve correlation is obtained by 
plotting Cowled’s data on a k/d versus Re,~'’* plot, as 
is shown in Fig. 15. It should be noted that A, the 
wavelength, is constant for Cowled’s data. This plot 
suggests that, for the wave type roughness, thickness 
ratio rather than height alone is of importance in de- 
termining whether or nor the wave will precipitate pre- 
mature transition. 


Effect of Pressure Gradient on Critical Roughness 


The various investigations show that pressure-dis- 
tribution effects are of minor importance in determining 
critical roughness. Critical values that have been ob- 
tained from the various investigations cover a wide 
range of pressure distributions. In the regions of 
laminar flow, the pressure distribution of the present 
investigation covers a range of Hartree’s 6 from —0.2 to 
+0.64. The regions where roughness was located 
covers a range of 6 from 0.025 to 0.64. Cowled’s 
roughness lay in a range of 8 from 0.4 to 0.6. Fage 
tested in a range from 0.02 to 0.1. As is shown in 


Figs. 10 and 13, the values of critical roughness are all 
in reasonable agreement. 

The range in pressure distributions that has been 
covered by the various investigations should be suf- 
ficient for any practical application of the roughness 
data, such as the design of low-drag airfoils. The 
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maximum value of 8 that is ever likely to occur on air- 
foils is about 0.55; e.g., an NACA 66,2-420 airfoil 
at an angle of attack of 16° gave a 8 of 0.57 on the lower 
surface in a wind-tunnel test. It is impossible to ob- 
tain values of 8 this high and simultaneously have lam- 
inar flow on both sides of an airfoil. Hence, the maxi- 
mum §’s of interest are even less. 


Effect of Tunnel Turbulence on Critical Roughness 


Two facts seem to indicate that the degree of tunnel 
turbulence has a minor effect on critical values of R,. 
First, satisfactory correlation of the data from various 
investigations has been obtained for roughness of a 
given shape. These investigations were carried out in 
a variety of tunnels, in which it may be assumed that 
the level of tunnel turbulence varied widely. Second, 
tests made in this investigation show that increasing 
the turbulence level tenfold (from 0.03 per cent turbu- 
lent to 0.30 per cent) increases the value of Ry oie ap- 
proximately 50 per cent (see Fig. 10), that is, reducing 
the number of turbulence screens from six to one to 
none increases the average R;,, .,;, for the wire configura- 
tion tested from approximately 130 to 165 to 195, re- 
spectively. 

The increase in the critical-roughness Reynolds 
Nuinber with increased tunnel turbulence is consistent 
with small-perturbation theory. First, the length of 
natural flow is less for increased turbulence, and thus a 
perturbation introduced by the roughness element has 
less time to grow before it reaches natural transition. 
Second, premature transition may be caused not only 
by the perturbation introduced by the roughness but 
also by the sum of this perturbation and those present 
due to free-stream turbulence. If the level of free- 
stream turbulence is increased, the perturbation con- 
tributed by the roughness elements may have to be in- 
creased correspondingly to cause premature transition. 


Composite Plot of All Critical-Roughness Data and Its 
Practical Utilization 


The results of the various investigations for all types 
of roughness are presented in Fig. 16. The critical 
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Fic. 16. Comparison of critical roughness values for various 
types of roughness of this and other investigations. Lines of 
Rx = constant represent Blasius boundary layer (6 = 0). 


values extend over a large range in R,, but the values 
may be grouped according to type of roughness. The 
most critical types are the two-dimensional elements, 
which have R; as low as 40 and which extend as high as 
R, = 400, with a mean of 125. The second most criti- 
cal elements are the three-dimensional elements, which 
fall in the range of R, .,i: from 100 to 800, with a mean 
of about 350. The least critical are surface waves, 
which have R; ,;i, extending from 1,000 to 4,000. 

Frequently, the predominant interest lies in the re- 
gion where the roughness has no effect on transition. 
The limit for this region of tolerable roughness can be 
obtained from Fig. 16. The lowest observed roughness 
is R, = 40, but to allow for a safety factor to account 
for shapes worse than those investigated, etc., the ‘‘safe 
region”’ of roughness is selected as that below R; = 25. 
Using R, = 25, the tolerable roughness was calculated 
for an airfoil and a body of revolution. Values for 7 
= 0.4 are presented in Table 1. The values indicated 
are for the most critical positions on the bodies, a point 
slightly forward of 1 per cent chord. 

At least one series of flight tests that has been made 
to determine the effect of isolated surface roughness on 
transition corroborates the values in Table 1. John- 
son,”! reporting the results of two flight tests made at 
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TABLE 1 
Maximum Allowable Roughness for M = 0.4 at Most Critical 
Point on the Body 


Altitude Altitude 


Altitude 


Configuration 10,000 ft. 35,000 ft. 60,000 ft. 
NACA €5-010 Airfoil, 
Ci =0 0.00092 in. 0.0017 in. 0.0042 in. 


Elliptical Body of Rev- 
olution of Fineness 


Ratio 9:1 0.00064 in. 0.0012 in. 0.0029 in. 


an altitude of 1,000 ft. and a Mach Number of 0.62, 
states that “‘... near the leading edge at about 2'/, 
per cent chord, a pimple of about 1 thou. in height 
causes transition while further aft at 20 per cent chord 
a pimple 3 thou. high is needed.” 

It is interesting to compare the critical roughness for 
transition in laminar flow with the roughness neces- 
sary to cause a drag rise in turbulent flow. In turbu- 
lent flow a surface is considered to be hydraulically 
smooth if kvx/v < 5 (see Schlichting®*). We follow 
Schlichting’s notation except to replace his k, by k; 
ve is the friction velocity defined as V¥a/ p. Since 
= p(Ou/Oy),., we have, remembering that u, 
R(OUu/OY) w, 


The quantity u,k/v is our usual roughness Reynolds 
Number, and, thus, the long established criterion for a 
hydraulically smooth pipe—namely, kvx/v < 5—is ap- 
proximately the same as the lower limit of critical rough- 
ness Reynolds Number (arbitrarily chosen as 25 above), 
a most interesting coincidence. Goddard?’ also noted 
this similarity. In Part (II) we shall return to this 
similarity in discussing the effect of compressibility on 
critical roughness. 


Roughness Used as a Tripping Device 


In wind-tunnel and towing-tank tests it is often de- 
sirable to fix transition at a given location. Roughness 
is the most common means of accomplishing this result. 
Here the vaiue of R; necessary to bring transition right 
up to the roughness is of importance, for the roughness 
should be no larger than that necessary to trip the 
boundary layer because it adds to the form drag and 
may have other adverse effects. The roughness neces- 
sary to move transition to the roughness itself is pre- 
sented in Fig. 17. The data of this and the other inves- 
tigations show that the values of R, necessary to trip 
the boundary layer are, approximately, 


Roughness Average R, 
Spanwise wire 300 
Row of 1/16-in.-dia. discs 600 
1/4-in. wide strip of sandpaper type 400 


R, increases above these values as the tripping device 
approaches very close to the leading edge, because of 
the criterion, noted earlier, that Rs > 320 for a fully 
turbulent flow. Preston"! notes that, when the natural 
Ro is so small that one must depend entirely on the drag 


It 
on t 
Rey 
tod 
cons 
deal 
deve 
sion 
sion 
beca 
dim 
dim 
rate 
two- 
sion 
thre 
of tk 
dime 
trem 
the ¢ 
thro 
Now 
para 
dian 
havi 
obvi 
ata 
Becz 
side 
gene 


238 

of t 

Ry 

i bow 

incc 

and 

lenc 

pert 

| mer 

| turt 

| ane 

oe flat 

° this 

4 

tet? 

corr 

bap | The 

6/6, 

2, 

e 

oy 

Ba, 8° Sond a 

ag 

AE 

W 

the 

amp 
ary | 

ficat 

tially 


itical 
ide 


in. 


THE SMALLEST HEIGHT OF ROUGHNESS 239 


of the tripping device to produce Reg > 320, then the 
R, of a wire-type tripping device required to trip the 
boundary layer is from 600 to 800. 

Hama’ presents ‘‘an efficient tripping device’ that 
incorporates the favorable properties of both the two- 
and three-dimensional elements for producing turbu- 
lence immediately. The two-dimensional element is 
favorable in that, for the same size, it produces a larger 
perturbation in velocity. The three-dimensional ele- 
ment is favorable in that it produces vortices that will go 
turbulent quicker than will two-dimensional disturb- 
ances. Hama’s device consists of a spanwise row of 
flat triangular patches lying on the surface with their 
apexes pointed upstream. Hama estimates the height 
Reynolds Number necessary to produce transition at 
this device to be 45. Values of k/0, and Re, taken from 
a figure of reference 13, in which turbulence is shown 
behind the trip, are plotted in Fig. 17, along with the 
corresponding values for the other types of roughness. 
The values plotted are not necessarily the lowest values 
of k/6, and R, producible by Hama’s device. 


Concluding Remarks on Critical Roughness in Low-Speed 

Flow 

It has been shown that the effect of critical roughness 
on transition is primarily a function of the roughness 
Reynolds Number, and, in turn, its value has been found 
to depend greatly on shape. These findings are not in- 
consistent with that part of boundary-layer theory that 
deals with the growth of small perturbations and their 
development into turbulence. Briefly, the two-dimen- 
sional roughness is more critical than the three-dimen- 
sional because it introduces larger disturbances. But, 
because the two-dimensional element produces two- 
dimensional disturbances, which must evolve into three- 
dimensional vortices before producing turbulence, the 
rate of movement of transition forward is less for the 
two-dimensional elements. In the preceding discus- 
sion it has been shown that as the shape ratic k/d of 
three-dimensional elements becomes small, the effect 
of these elements on transition approaches that of two- 
dimensional elements. Now what about the other ex- 
treme, when k/d becomes very large—that is, what is 
the effect of a fine wire normal to the surface, extending 
through the boundary layer far into the free stream? 
Now, instead of a height parameter like R,, the major 
parameter is undoubtedly Rz = ud/v, the local wire- 
diameter Reynolds Number. From the well-known be- 
havior of wakes and vortex streets behind cylinders it is 
obvious that turbulence will be generated when Ry 
at any point along the wire exceeds a certain value. 
Because u is zero on the surface and increases to LU’ out- 
side the boundary layer, turbulence should first be 
generated outside the boundary layer. 

When a wire is laid spanwise along the model, within 
the boundary layer, the disturbances it creates are 
amplified by the velocity distribution within the bound- 
ary layer, perhaps by the Tollmien-Schlichting ampli- 
fication mechanism. Certainly the disturbance ini- 
tially caused by the wire becomes stronger downstream. 
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ness when transition occurs at the roughness itself. Lines of Rx 
= constant represent Blasius boundary layer (8 = 0). 


But the disturbance created by a long thin wire per- 
pendicular to the model acts in a quite different fashion. 
If the wire does create some turbulence, this turbulence 
is generated substantially outside the boundary layer 
and directly behind the wire. Unless the turbulence is 
developed immediately, it will not develop at all, be- 
cause disturbances introduced by the wire will be 
damped in the free stream, more or less in accordance 
with well-known laws of turbulence decay. 

This idealized example, together with that of a con- 
ventional two-dimensional wire type of roughness, 
helps to bound the types of effect caused by different 
kinds of roughness. On the one hand, for the perpen- 
dicular wire, if turbulence is present at a downstream 
point, it is obvious that stronger and stronger indica- 
tions will be found forward of this point right to the 
wire. On the other hand, for the wire on the surface 
of the model and immersed in the boundary layer, 
turbulence may be present at a downstream point but 
not slightly forward of this point. 

The downstream development of disturbances caused 
by other shapes of roughness will be somewhere between 
those of the examples above. The disturbed region 
caused by a three-dimensional roughness high in the 
beundary layer may be more similar to that behind the 
perpendicular wire. The disturbed region caused by a 
three-dimensional roughness that lies very low in the 
boundary layer may be more similar to that behind the 
conventional two-dimensional wire. In order to under- 
stand fully the effect of roughness on transition, a de- 
tailed research on the flow around and behind roughness 
elements is necessary. Some studies of the micro- 
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Fic. 18. Variation of critical roughness Reynolds Number with 
the local roughness Mach Number in supersonic flows. 


mechanism of the disturbed flow behind an element 
have been made (see references 6, 7, 10, 12, and 24), 
but many questions still remain unanswered, particu- 
larly the dynamic aspects dealing with wave motion. 

In retrospect, when a wide range of conditions is con- 
sidered, it is clear that no better correlation parameter 
exists than R, for determining the degree of roughness 
necessary to affect transition. While there is a great 
deal of scatter in experimental values of R,; ¢,i, the 
value is not affected by the several variables considered 
in this investigation. The large spread of values is 
much less when considered in terms of R,'/?, a quantity 
more nearly proportional to the height of the roughness. 
In brief, R, is a very sensitive parameter. 

As was noted in the Introduction, Schiller was first to 
propose the criterion that R, .,;. = constant. Further- 
more, purely from considerations of the wake behind 
cylinders in a uniform stream, Schiller proposed that 
Rx crit = 50 was about the value to be expected for a 
wire type of roughness. 


PART (II) EFFECTS OF MACH NUMBER 
ON CRITICAL ROUGHNESS 


Useful data on critical roughness at high Mach Nuim- 
bers are extremely limited. From wind-tunnel tests, 
it is known that the height of a roughness trip or the 
roughness Reynolds Number defined as kU’/vs neces- 
sary to cause transition increases markedly with Mach 
Nuinber.* This increase is not surprising for the fol- 
lowing reasons: first, the boundary-layer thickness in- 
creases with Mach Number for a given U’/v; second, 
the wall temperature rises as J/° (unless the wall is 
cooled), and thus the viscosity increases; and, third, 
the Reynolds Number above inay be written as 


Ri, = RU/vs = (1/0.663)(U/a)(k/d) 


where v; = 0.663 ad (by basic kinetic gas theory), a = 
speed of sound, and X = mean free path of molecules. 
If R,, were to remain constant for all Mach Numbers, 
k/X could become quite small. This is not physically 
reasonable. 

An approach to the effect of Mach Number on critical 
roughness is suggested by the work of Liepmann and 
Goddard.”**6 They conducted tests to determine the 
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effect of surface roughness on skin friction in turbulent 
flow and found that the ‘‘critical roughness’ below 
which the surface is hydraulically smooth is a function 
of kvx/v,. only, just as in the incompressible case. 
Furthermore, the critical roughness below which the 
surface is hydraulically smooth is 


Rerit =~ 10», Ux 


and is independent of Mach Number. This is approxi- 
mately equal to the thickness of the laminar sublayer 
for a smooth surface for both compressible and incom- 
pressible flows. As noted above in the Section ‘‘Com- 
posite Plot of All Critical-Roughness Data and Its 
Practical Utilization,”’ 


kvx/v = V /p 


for the low-speed flow. This fact and the work of 
Liepmann and Goddard on turbulent boundary layers 
suggest the local Reynolds Number, 


R.. = ‘yy. 


as the parameter to correlate the effect of Mach Num- 
ber on critical roughness. Then, if there is any effect 
of Mach Number, it will be as a function of local Mach 
Nuinber, 


M,. We: ‘ay 


The local values of R, and 1, must be based on proper- 
ties in the neighborhood of the roughness. The top of 
the roughness is both a convenient and a logical refer- 
ence point; therefore, hereafter the subscript ‘‘k’’ will 
denote properties at the top of the roughness. 


CRITICAL ROUGHNESS AT SUPERSONIC SPEEDS 


The only available sources of useful supersonic data 
on critical roughness that enable the determination of 
R, = u,k/v;, are references 20, 27, and 28. The results 
of reduction of the data on wires are shown in Fig. 18. 
The calculations were made using the boundary-layer 
temperature and velocity profiles calculated by the 
Chapman and Rubesin method.?* Sutherland’s re- 
lation for the variation o; viscosity with temperature 
was used. The data of reference 27 were obtained for 
wire-type roughness on a cone whose surface was either 
insulated or cooled. In the case of cooling, the rough- 
ness element was always in front of the region of maxi- 
mum cooling, which caused relatively large tempera- 
ture gradients at the roughness location. These tem- 
perature gradients were approximated from the data 
presented in reference 27 and used in the determination 
of the velocity and temperature profiles by the method 
of reference 29. The NACA data” were also obtained 
for wires but on an insulated surface. Fig. 18 shows the 
variation of R, with where = u;,/a,. It should 
again be emphasized that the properties u,;, »,, and a, 
are values calculated at a height k in the boundary layer 
when roughness is not present. 

At roughness Mach Numbers equal to or greater than 
1, one would expect a shock wave to occur in front of or 
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on the roughness element—e.g., see Fig. 18 of reference 
27 and the various photographs presented in reference 
32. If the shock wave produced by the roughness is of 
the form shown in Fig. 19—that is, the shock is normal 
at the top of the roughness, changing to an oblique 
shock above the roughness—then the flow conditions 
at the top of the roughness are thcse behind a normal 
shock. The characteristic roughness Reynolds Num- 
ber is now R,, = uy,k/vy,, where the subscript 2 refers 
to conditions behind a normal shock at /,. A calcu- 
lation of R,, gives the broken line in Fig. 18 for 1/; > 1. 
It should be pointed out that, in the supersonic region, 
R, is plotted against J/,,, not \/,,. The latter would, 
of course, always be equal to or less than 1. Thus, as- 
suming a normal shock when greater than 1, Fig. 18 
shows that in the subsonic region R,, ,,;; increases nearly 
as \/,,°, dips at W,, > 1, then again increases as .1/,,”. 

It is of interest to note that a Reynolds Number 
R,,*, based on velocity and temperature downstream of 
a normal shock and pressure upstream of the normal 
shock, has nearly a constant value for the limited 
amount of data on critical roughness at supersonic 
M,,.. This is shown by the “‘star’’ points in Fig. 18, 
where R,.,* is defined by 

It is known from experimental work on laminar flow 
over small steps or roughness elements (for example, 
see reference 32) that the pressure measured on the sur- 
face about the element does not necessarily change 
greatly even though shocks are observed at the step or 
roughness element. This would be expected both be- 
cause the shock almost immediately evolves into a 
simple Mach line and because the extent of the shock 
is too sinall to cause a general pressure rise. 

A still more nearly constant value of R; throughout 
both the subsonic and supersonic regime is obtained if 
R,, is computed using (1) stagnation temperature, (2) 
static pressure upstream of any shock, and (3) a veloc- 
ity “,, determined by normal shock relations. Use of 
upstream static pressure is plausible, as has already been 
noted. Furthermore, use of stagnation temperature is 
plausible because the fluid in contact with the roughness 
will be close to this temperature. Unfortunately, no 
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Fic. 19. The approximate physical process by which a rough- 
ness produces turbulence in a boundary layer. 
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uce transition at the roughness, as a function of the local rough- 
ness Mach Number. 


plausible physical justification for using a velocity «4, 
downstream of a normal shock has been found. These 
remarks and facts are included in the hope they will 
aid other investigators in studying the effect of Mach 
Number. 

Fig. 1S shows that the critical values of R, for MM, 
less than 0.5 fall below those values of R, observed by 
the authors for low-speed flows. Possibly the explana- 
tion is as follows. Both Brinich and van Driest give 
the nominal heights of their roughness elements. The 
calculation of their R,’s was based on these heights. 
Our experience shows that the maximum height is 
undoubtedly higher than the nominal height. A 
greater height would cause the values of both R, and 
AM, to increase; and the increase would follow closely 
the line of data points in Fig. 18, since the variation of 
local yu, v and M in the boundary layer closely follows 
this line for the free-stream conditions of the points 
shown. 

A similar variation of R, with .1/, is observed for 
roughness elements necessary to cause transition at the 
element. References 20, 28, and 30 present observa- 
tions on such roughnesses at supersonic speeds. The 
variations of R,, or R,, with 1/7, for these observations 
are shown in Fig. 20. When the roughness element 
protrudes completely through the boundary layer, rela- 
tively large values of R,, were observed——for example, 
note Brinich’s data in Fig. 20. 


DISCUSSION 


The concept that the roughness senses the purely 
local flow, which of course implies conditions down- 
stream of a shock for supersonic flow, is in agreement 
with observations made with hot wires. In a study of 
heat transfer from fine wires in supersonic flows, Laufer 
and McClellan*! found that the Nusselt Number was 
a unique function of the Reynolds Number when both 
of these parameters were based on properties behind a 
normal shock. In other words, heat-transfer proper- 
ties of a small hot wire are independent of Mach Num- 
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ber when based on fluid properties downstream of a 
normal shock. This statement was well verified by the 
experimental data presented in reference 31. 

If it is separation from the roughness that produces 
vortices or disturbances that cause transition, as was 
suggested above for low-speed flow, perhaps an idea 
of the effect of Mach Number on critical roughness 
Reynolds Number can be obtained by studying the 
effect of Mach Number on laminar separation on bodies 
at very low Reynolds Numbers or by studying separa- 
tion caused by buimps or steps on a surface in super- 
sonic flows. The authors know of no experiments 
which give the effect of Mach Number on separation 
from bodies at low Reynolds Numbers, but separation 
caused by bumps or steps on a surface have been 
studied by Chapman, et al.** and Greber, et al.** Ref- 
erence 32 shows that the Reynolds Number Ryy = 
AX U/v, based on the length of the separated region 
AX in front of a bump (see Sketch), which is just large 
enough to cause turbulent reattachment, increases with 
the free-stream Mach Nuinber, and that this increase is 
similar to the variation of critical R,, with 1/,, found 
for the roughness wires above. A calculaticn of Ryx, 
= U,, AX/v for points where 1; > 1 shows that 
increases nearly as but dips at J/; close to 1. 
Although the Reynolds Number Ryx, is much larger 
numerically than the roughness Reynolds Number R,, 
the variation of the two Reynolds Numbers with local 
Mach Nuinber is similar. It should be noted that the 
height of the steps of reference 32 was considerably 
greater than the boundary-layer thickness ahead of the 
separated region, and also that if a shock was caused 
by the step, it formed immediately in front of or on the 
step but behind the separated region. Since the sepa- 
rated region (of length AX) was in front of the shock, 
the only reason for discussing the variation of Ryy, = 
AX U2/v2 with M,, is that the length of the separated 
region gives an idea of the magnitude of the disturb- 
ance. 

Greber, et al.,** have made experimental and theoreti- 
cal studies of the interaction of oblique shock waves 
with laminar boundary layers, and these studies may 
help to explain the dip in R,,, at 17, slightly greater than 
1. They have been able to predict the pressure rise 
that causes separation and also the length of the sepa- 
rated region. Furthermore, they point out that the 
conditions near separation are independent of the 


JOURNAL OF THE AERO/SPACE SCIENCES—APRIL, 1959 


agency that induces separation—e.g., a shock wave or 
a step in the surface—and they also point out that the 
separation phenomena are essentially self-similar. 
That is, all significant parameters at separation are 
functions of the local Mach Number and Reynolds 
Number. The criterion for separation which they 
present can be written as 


Cy, = 1)? 


when C,, = pressure rise at separation and (Cy), = 


(7/q)s, Where subscript s denotes values evaluated at the 
height of the separated region, or, in the case where 
separation is caused by a surface protuberance, the 
value of 7 /q at the height of the protuberance. Clearly, 
this expresssion for C,, peaks at J = 1. 

In the case of a normal shock at a surface protuber- 
ance of height k, it can be shown that 


R, = AM,4/(M,? — 


where A is a constant. 

Thus, we see that in supersonic flows—i.e., 1/, > 1— 
the pressure rise necessary to cause separation gives a 
peak in R, at J, = 1. This, of course, does not ex- 
plain the increase of R; with Mach Number for sub- 
sonic or supersonic flows. It is known that the stabil- 
ity of a laminar separated region increases markedly 
with Mach Number.*° 

Before summarizing what is known of the effects of 
Mach Number on critical roughness, an illustration will 
be given in which the variation of R, with J, pre- 
sented in Fig. 1S helps explain an observed phenomenon. 
The NACA has observed that, under conditions of ex- 
treme surface cooling, the transition Reynolds Number 
decreased with decreasing surface temperature.**: ® 
This is in disagreement with the prediction of stability 
theory and with what has been observed for moderate 
surface cooling. The authors of reference 34 question 
whether or not transition reversal can be attributed 
primarily to surface roughness, their reasoning being 
that cooling produces only a small change in 6. 
Calculations based on data presented in references 34 
and 35 and an assumed & show that in one case in which 
cooling produced only a 50 per cent increase in k/6, it 
produced a tenfold increase in R,. It appears that if 
the criterion for critical roughness is R, as presented in 
this Section, transition reversal due to extreme surface 
cooling can be attributed to surface roughness. 

Summarizing, it can be said that the limited amount 
of data on critical roughness at high Mach Number is 
not inconsistent with the mechanism, suggested above, 
by which the roughness causes transition in low-speed 
flow—that is, critical roughness is that which first 
introduces local separation or vortices. Calculations 
based on this principle and the variation of R, with /; 
presented in Fig. 18 (the broken line) give tolerable 
roughness that would be encountered in flight. Such 
tolerable roughnesses are presented in Fig. 21, both for 
an insulated flat plate and for a flat plate with a wall-to- 
free-stream temperature ratio of 1. The curve of Fig. 
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18 was extended linearly to make the calculations at the 
higher Mach Numbers. The Figure shows that toler- 
able roughness first decreases with Mach Number and 
then increases because of aerodynamic heating, particu- 
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Roughness Average R,, 

Spanwise wire 300 

Row of 1/16-in.-dia. dises 600 

1/4-in. wide strip of sandpaper type 400 


larly because of the effect of temperature cn viscosity. 
Also the Figure shows the effect of surface temperature 
on the tolerable roughness at 1 in. from the leading 
edge; as the Mach Number is increased in the sub- 
sonic region, the curves of tolerable roughness for the 
two surface temperatures (the insulated wall and 7, = 
T;) diverge because of aerodynamic heating. At 
higher Mach Numbers further increases in temperature 
and a corresponding increase in viscosity are caused by 
the shocks produced by the roughness. The increase 
in viscosity along with theincreasing Mach Number 
increases the tolerable roughness until eventually the 
roughness extends through the boundary layer into the 
free stream. When the roughness extends far enough 
R, becomes independent of wall conditions, as is indi- 
cated in Fig. 21. Also, the tolerable roughness be- 
comes independent of its position on the plate if it ex- 
tends far enough out of the boundary layer. Further 
experimental data on critical roughness in supersonic 
flows is needed to validate the values of tolerable 
roughness presented in Fig. 21 and also to validate the 
mechanism on which the calculation of these tolerable 
roughness values was based. 


CONCLUSIONS 


This investigation to find the smallest height of 
roughness capable of affecting boundary-layer transition 
leads to the following conclusions. 


Low-Speed Flow 


(1) For a given shape of roughness, the value of 
R; crit has been found to be relatively constant and to be 
substantially independent of both height and location of 
roughness, the pressure distribution, and the degree of 
free-stream turbulence. 

(2) The shape of roughness has a major effect on the 
values of R, The value of R; increases as the 
type of roughness changes from a three-dimensional 
type to a two-dimensional type. The present investi- 
gation established values of R, .,:: for three types of 
roughness as follows: 


Range of 

Roughness 
Spanwise wires 40-260 
Row of 1/16-in.-dia. discs 100-550 


1/4-in. wide strip of sandpaper type 178-330 


(3) The results of the various investigations covering 
all types of roughness show that the smallest roughness 
that has any effect on boundary-layer transition is that 
corresponding to a roughness Reynolds Number of 
about 40. 

(4) The values of R, necessary to move transition to 
the roughness itself are, approximately, 


High-Speed Flow 


(5) The small amount of available, usable test Cata 
shows that R, may be the characteristic parameter for 
correlating the effects of roughness on transition, just 
as in the subsonic case. The value of R, varies with 
the local Mach Number J/, at the top of the roughness. 
For \/, greater than 1, the Reynolds Number R,,, 
based on conditions behind a shock, is perhaps the best 
characteristic parameter. The value of R, increases 
with ./, and, if a normal shock is assumed, the value 
peaks at 1/, = 1, decreases for \/;, slightly greater than 
1, then again increases for larger .\/;, approximately as 
M,?, according to the sparse ainount of data. The 
value of R,. .7;: for a wire at supersonic ./, is lower than 
the corresponding value for a three-dimensional rough- 
ness, just as in the subsonic flows. Further studies of 
critical roughness in supersonic flow are needed. 

The Physical Process—It is believed that the domi- 
nant process by which isolated, protruding types of sur- 
face roughness first affect transition is now essentially 
established, at least in the gross or first-order sense. 
Premature transition is caused when the roughness Rey- 
nolds Number becomes sufficient to develop a dead- 
water region and/or vortices behind the roughness. 
Three parameters seem to determine when this region 
will just begin to form—a Reynolds Number based on 
strictly local properties of the flow, a Mach Number 
based on the same local properties, and, finally, the 
shape of the roughness. All other factors, such as free- 
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Fic. 21. The variation with Mach Number of tolerable rough- 
ness on a flat plate at 1 in. and 12 in. from the leading edge. 
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—— FLOW DIRECTION 
-SONIC POINT 


Fic. A-1. Approximate location of sonic point on several typical 
bodies. 


stream Mach Number or pressure distribution, enter 
into consideration only to the extent that they deter- 
mine the local fluid properties at the reference point on 
the roughness. To first order, none of the observed 
facts are inconsistent with this model. Flow pictures 
of separation behind rods at very small Reynolds Nuin- 
bers, Laufer’s observations on hot wires, Goddard’s 
observations on roughness in compressible turbulent 
flow, and Chapman’s and Greber’s considerations on 
flow separation all reinforce the concept. The approxi- 
mate physical model is illustrated by Fig. 19. 


APPENDIX—AN IMPROVED PHYSICAL MODEL FOR THE 
COMPRESSIBLE FLOW PAST AN ISOLATED BIT OF 
ROUGHNESS 


In the main body of this article, in the section of the 
effects of Mach Number, two methods of computing a 
roughness Reynolds Number were introduced. One 
used the height k and the density, velocity, teimpera- 
ture, and viscosity computed for conditions behind a 
normal shock, corresponding to the appropriate stream 
filament at the height k&. The second used the velocity 
and temperature downstream of a normal shock but the 
pressure upstream. 

Some criticisms of the first method are: (1) no 
normal shock is observed to exist near a roughness ele- 
ment low in the boundary layer, and (2) there is no 
observable pressure rise near a roughness element; in 
fact, the main effect of a small bit of roughness is to 
create a Mach line. 

With respect to the second method, criticism (1) 
continues to apply, and, in addition: (3) the assump- 
tions do not fit together into any known physical model 
or flow situation. 

A model that is much nearer physical reality has since 
been conceived; it uses sonic velocity as the reference 
velocity and recognizes the fact that the roughness is a 
small bluff body. When M;, exceeds some high sub- 
sonic value, the flow becomes transonic and remains 
transonic as M,, is indefinitely increased. For a cone- 
cylinder body as well as a two-dimensional shape having 
the same profile, sonic velocity occurs at the intersec- 
tion of the nose with the afterbody. Hence, at the 
maximum breadth of these bodies a characteristic 
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velocity is the sonic velocity. This breadth can be 
thought of as corresponding to k. In general, in the 
transonic flow regime, for any type of nose the maxi- 
mum velocity is approximately sonic, and this velocity 
occurs close to the top of the nose. Liepmann and 
Roshko* discuss the problem of bodies having wedge 
and conical noses. Fig. 5 of Bogdonoff*” shows the 
sonic point for several axially symmetric and two- 
dimensional bodies. Lees* and Sibulkin® provide 
good source material for the theory of the sonic-point 
location and the flow in its vicinity. 

The essence of these considerations is that, because 
of the blunt nature of the roughness, the sonic point 
will always be near the top of the roughness, as illus- 
trated in Fig. A-1, and the maximum velocity never 
differs greatly from the sonic value, so that the sonic 
velocity appears to be a reasonably representative meas- 
ure. Then a suitable characteristic Reynolds Number 
is one formed as follows: (a) the characteristic height 
is k, (b) the characteristic velocity is equal to sonic 
velocity appropriate to the streamline at the height , 
(c) the pressure is equal to static pressure, and (d) the 
temperature is equal to the sonic temperature appro- 
priate to the streamline at the height k. The viscosity 
is determined from the temperature (d). The density 
is determined by means of the perfect gas law and the 
pressure (c) and temperature (d). Because the rough- 
ness is very small and because it is in the boundary 
layer, the use of static pressure is believed justified. 
No large region of locally high pressure can exist, be- 
cause pressures can be propagated by way of the bot- 
tom of the boundary layer. The formula for R, so 
defined, designating sonic conditions by ( )*, is 


= p(pr, (A-1) 


The symbol R,* should not be confused with R,,* in 
the main body of this article. 

Fig. A-2 presents the data of Fig. 18 reduced accord- 
ing to Eq. (A-1). The dashed line is taken from Fig. 
18 to indicate the previous normal-shock results. The 
critical Reynolds Number is indeed nearly constant. 
In addition, under the present hypothesis the high- 
speed data now appear more consistent with the low- 
speed results. Under the present system of data re- 
duction we refer to a velocity (a*) in the immediate 
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vicinity of the top of the roughness. The low-speed 
results used u,, a velocity at the correct height but far 
upstream of the roughness. Were we to use a local 
velocity in the immediate vicinity of the top of the 
roughness, considerations of the flow about circular 
cylinders adjacent to a wall indicate that the local 
velocity would be two to three times u,. Then R,!/? 
so evaluated would be some 50 per cent greater—that is, 
it would be in line with the high-speed data. 

The differences in concept of this Appendix, that of 
the main report, and the older u;k/v; formula are 
worthy of further explanation. The older wk v; 
formula referred to properties in the free stream far 
away in both the x and y directions. Roughly speak- 
ing, the u,k/», correlation refers to properties at the 
correct y elevation in the boundary layer but still far 
away in the x direction. The present correlation re- 
fers to the properties of the flow in the immediate 
vicinity of the roughness—that is, at about the correct 
value of both x and y. The only reason a correlation 
of this sort becomes possible at all in such a compli- 
cated, viscous, transonic flow is that velocities near the 
top of the roughness are known not to differ greatly 
from a*. 

The flow that is actually in contact with the rough- 
ness has come from the wall upstream, and of course 
the roughness has both physical and thermal contact 
with the wall. Therefore, Eq. (A-°) is not the only 
one possible under the present concept. Recovery 
temperature, wall temperature for cooled walls, and 
total temperature are all possible characteristic tem- 
peratures. However, until more data are available 
there is no way to discriminate among the various possi- 
bilities. 

An interesting deduction can be made from Eq. 
(A-1) for the special case of flow past an insulated plate 
of a gas whose Prandtl Number is unity. For these 
conditions 7* and a* are constant on all streamlines 
and R, as evaluated by Eq. (A-1) becomes independent 
of the shape of the boundary-layer velocity profile! 

The roughness Reynolds Number defined as R, = 
usk/vs is also independent of the shape of the velocity 
profile. Then, for the special case under consideration, 
assuming perfect gas relations and up ~ 7°."6, it is easy 
to compute the ratio of R, te R,*. When the Mach 
Number is great enough to develop sonic velocity at 
the roughness—that is, 1/7, > 0.7 approximately—we 
find 


1 — 1)/2]M,?]!-26 
[(y + 1)/2] 

At Mach 5 with y = 1.4 the ratio is about 38:1. 
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Some Exact Solutions to the 
Magnetohydrodynamic Equations for 
Incompressible Flow 


F. D. Hains 
Boeing Scientific Research Laboratories, Seattle, Wash. 
October 20, 1958 


— SOLUTIONS of the magnetohydrodynamic equations are 
considered for incompressible, nonviscous flow. The flow 
field and magnetic field are two-dimensional and functions of 
x and yonly. The steady flow field with magnetic interaction is 
assumed to be established so that no consideration is given here 
to the magnetic field when there is no flow or to the flow when the 
magnetic field is removed. In the subsequent analysis the con- 
ductivity o is constant. 
The governing equations are :! 


uuz + vu, + pr/p = —J3He/p (1) 
uv, + ov, + p,y/p = J3Hi/p (2) 

uz +o, = 0 (3) 

J; = + uH, — Ji = = 0 (4) 
My, + Hy = 0 (5) 

4nJ; = Hx, — (6) 


where the subscripts 1, 2, and 3 indicate components in the x, y, 
s directions, respectively. 
Introducing the stream function y, p can be eliminated from 
Eqs. (1) and (2) to yield 
Vil + Yrrr) Vil + Very) —H2J3y/p J3z/p (7) 
A solution is assumed in the following form: 
= (Gx + Cr)g(y) 
HM, = (Gx + C.)P"(y) 
H, = F(y) 
E; = (Cx + C:)M(y) 


(8) 
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(4) 
(5) 
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— 


Using Eqs. (3), (4), and (8), 


where C; and C2 are constants. 
Eq. (7) becomes 
— = (g'F — CP’g)F’ + 

2P'FCig’ + C,P'M + — 


4roCigF(M + g'F + CigP’) + FM’ (9) 
Eliminating J; from Eqs. (4) and (6) 
F = — CP’ — M)/g’ (10) 
and combining this with Eq. (5) 
+ [4aoCig — — 
(4aogg’’/g') CP’ + 420M’ — =0 (11) 


If one of the three quantities g, P, or ./ is specified, the re- 
maining two are found from Eqs. (9) and (10). This set of equa- 
tions will also give the trivial solutions J; = o£, which means 
that the streamlines coincide with the magnetic field lines. 

The stagnation-point flow is g = Ay + B where A and B are 


constants. Eq. (9) is simplified to 


[AF — + B)|F’ + 2P’FCQ,A + CP'M 
F(Ay + B) [M+ AF + C(Ay + B)P’] + FM’ 
(12) 
Integrating Eq. (11), 
M =(-P"/4ae) — C(AY + B)P’ + GAP + C; (13) 
so that Eq. (10) becomes 
F = -—C,P + C/A 
Substituting Eqs. (13) and (14) into Eq. (12) results in a 
third-order ordinary differential equation in P which can be 


solved numerically. 
A simple solution to this equation is obtained by taking 


(14) 


M=M= const. } 
F = Fy = const. } (15) 
Cc, =0 \ 
which automatically satisfy Eqs. (12) and (14). Eq. (13) yields 
P" = — C3) (y + C;) (16) 
so that the solution is 
u = 
v =0 
E=M™M (17) 
= My — C2) (y + cf 
He = Fy 
This is a uniform flow with parabolic magnetic-field lines. The 


pressure distribution can be obtained from the momentum equa- 
tions and the temperature distribution from the energy equation. 


REFERENCE 
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The Limiting Circulatory Lift of a Wing of 
Finite Aspect Ratio 


Barnes W/. McCormick 

Head, Aeronautical Engineering Department, School of Engineering, 
University of Wichita, Wichita, Kan. 

November 3, 1958 


XPERIMENTAL EVIDENCE INDICATES that the circulatory lift 
E of a jet-flapped wing of finite aspect approaches a maximum 
as the momentum coefficient of the jet is increased. Reference 1 
suggests that this phenomena is the result of strong tip vortices. 
In that same reference, an experimental curve is given which 
shows that the maximum circulatory lift of a finite jet-flapped 
wing varies nearly linearly with the aspect ratio. The analytical 
study presented here confirms these experimental results. It 
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4, 


Fic. 1. Wing section. 


will be shown that the maximum lift of a finite wing due to circu- 
lation, produced by any means whatever (suction BLC, blowing 
BLC jet-flaps, ordinary flaps) increases linearly with increasing 
aspect ratio. This behavior is attributed to the severe vertical 
displacement of the trailing vortex sheet at the high lift coeffi- 
cients. This vertically displaced vortex sheet not only decreases 
the angle of attack of each wing element, as calculated in the 
usual lifting surface theories, but also subtracts from the incom- 
ing free-stream velocity, which effect is usually ignored in finite 
wing calculations. For ease of calculation only the case of the 
elliptic wing is treated which is sufficient for demonstration pur- 


poses. 


THEORY 


The lift coefficient of a wing with a jet-flap can be expressed 
as the sum of two terms. 

Cr = Crp + Cy sin (a + 8) 

Cir is the dimensionless lift due to the circulation about the 
airfoil, produced by the angle of attack and the increment due 
to the jet. The remaining term is that dimensionless lift due to 
the momentum reaction of the jet in the vertical direction. 

The bound circulation, [, which may be abnormally high be- 
cause of the jet-flap effect, must satisfy Helmholz’s theorem of 
vortex continuity resulting in a strong trailing vortex sheet. 
As stated previously, it is usually assumed that this vortex sheet 
lies in the plane of the free-stream velocity, and its effect on the 
wing sections is calculated only as a decrease in the angle of 
attack of the wing sections. 

Consider now a section of the trailing vortex sheet and wing 
The vortex sheet is displaced downward due 
This angle of displacement in- 
creases as one moves downstream of the wing. Defining the 
downwash to be normal to the vortex sheet, the vortex sheet will 
be displaced from the horizontal at the wing by an angle, @, given 
by: 


shown in Fig. 1. 
to its own induced velocity. 


(1) 


By the law of Kutta-Joukowski, the resultant force in a given 
direction is equal to the product of the density, the resultant ve- 
locity normal to the direction, and the circulation I, or in vector 


= sin=! (w/V) 


notation 
eV XT (2) 


normal to the free-stream velocity; thus, 


F= 


By definition, lift is 
from Fig. 1 
L= p(V —wsindr 


— (w?/V)Ir (3) 


In the general case, the calculation of w and 6 would be lengthy ; 
however, consider the elliptic distribution of [ given by: 


r= — (y/b)? 


For this case, the relationship between w, and the wing section 
lift coefficients are well known. 


(4) 


(5) 


w = constant = 


w 
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Fic. 2. Comparison of theory and experiment. 


Consider, therefore, the mid-span section. 
Lo = — (w/V)?] 
By definition: 
Cip = Cy = Lo/(1/2)pV2Cy = 2(To/coV) [1 — (1/4) (To/bv)?] 
Or, since the aspect ratio is given by: AR = 4b/ me. 
Cr = [1 — (1/4) (To/dbV)?] (6) 


From this equation, it can be seen that C; reaches a maximum 
at some value of [)/bV. To find the value of Cz, Eq. (6) is 
differentiated with respect to T)/bV and equated to zero. 


Therefore: = 2/V3 (7) 
Thus: Chaz = AR 
~121R (8) 


From Fig. 2 it can be seen that Eq. (8) is in substantial agree- 
ment with tests performed by the NACA ona series of rectangular 
wings, as reported in reference 1. 
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A Creep Instability of Thin-Walled Tubes 
Under Internal Pressure 


lain Finnie 
Shell Development Company, Emeryville, Calif. 
November 10, 1958 


Ww" A THIN-WALLED TUBE of circular cross section, under 
constant internal pressure, deforms due to creep there is a 
continuous increase of radius and decrease of wall thickness. The 
stresses in the tube thus grow with time, which leads to an in- 
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creasing rate of deformation. The purpose of the present note 
is to point out that this effect leads to a maximum tube life which 
can be predicted from tensile creep data. 

The average tangential stress in the tube wall may be written 
as 


op = PR/h = (PR,/ho)et 


where P is the internai pressure, R and h the inside radius and wall 
thickness of the tube (original values Ro, fo), and ep and ep are 
the true (natural) tangential and radial strains. Alternatively, 
we may write the previous expression as 

= ore?“ (1) 
where o7, is the initial average tangential stress. 

It is found by the experiment! that in creep, as in other forms 
of plastic deformation, the volume of the deformed material does 
not change. For the tube under internal pressure there is no 
rotation of the principal strain axes during deformation, and the 
condition of volume constancy may be written as 


aE tea te=0 


In this equation, e, is the axial strain, and we have assumed that 
elastic strains are small enough, relative to creep strains, for elastic 
volume changes to be neglected. 

In an internal pressure test* under creep conditions, it is found 
that the axial strain is extremely small relative to the other prin- 
cipal strains.? Utilizing this fact, we have ey = —ep, and hence 
Eq. (1) becomes 

or = one (2) 

To relate stress and strain, we must now make some assump- 
tion about the creep behavior of the tubing. Many materials, 
when loaded by constant stress, show a creep rate which is ini- 
tially quite high but soon decreases to a steady value. This 
steady-state creep rate stays almost constant until rupture 
occurs.t When considering fairly large strains, as in the present 
case, it is thus often a good approximation to treat the creep rate 
as a function only of stress. The empirical relation which usu- 
ally gives a satisfactory fit to constant-stress creep data is 

é = Bo" 
where o and é are the stress and steady-state creep rate, and B 
and v are constants which may vary with temperature. 

For a tube under internal pressure, it may be shown that the 
tangential creep rate is given by 

ép = 


Substituting into Eq. (2) we obtain 
ép = ene”? 
where é7, is the initial tangential creep rate. Integration leads to 
= 1 — 2nér7,t (3) 


and it is clear that finite tube deformations are only obtained if 
t <1/2né7,. Thus, an upper limit to the life of the tube may be 
Ductile materials, in short-time 
Under 


obtained from creep data alone. 
rupture tests, show a rather large elongation to failure. 
these conditions, we would expect that ¢ = 1/2nér, would be a 
good approximation to the tube life. For specimens which fail 
at longer times, the elongations to fracture are generally smaller, 
and we would expect the upper limit (¢ = 1/2né7,) to consider- 
ably overestimate the rupture life. 

To complete the picture, we will now examine the conditions 
for rupture under increasing stress. We assume that constant- 
stress tensile rupture data may be fitted by the expression 


t = 


* Strictly speaking, we should specify that we are considering the case in 
which the axial load due to internal pressure is carried by the tube wall. 
This is almost always the situation in practice. 

+ The accelerating creep rate (tertiary creep) observed toward the end 
of constant-load tensile rupture tests is largely due to the increased stress 
produced by decreasing specimen diameter. However, it appears that 
limited tertiary creep also occurs in constant-stress tests.? 
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where ¢ is the rupture life and A and g are constants which may 
depend on temperature. If f is the rupture life which would 
correspond to the initial stress and ¢’ is the rupture life which 
corresponds to the varying stress, we write 


t’ 
J, (dt/t) = (dt/Ao™*) = 1 
0 0 


This assumes, following Robinson,‘ that in each increment of 
time a percentage of the rupture life is used up which depends 
only on the stress during that increment. When the accumulated 
total of the rupture life fractions equals unity, rupture is assumed 
tooccur. The preceding equation may also be written as 

(a/ao)* dt = to (4) 

0 
This equation applies only to varying uniaxial stress, and we 
should now consider the problem of predicting the life of a tube 
under internal pressure from tensile rupture data. The best ex- 
perimental evidence’ available at the present shows that for duc- 
tile metals the maximum tensile stress in a multiaxial stress sys- 
tem is the criterion of failure. Since the maximum tensile stress 
in a tube under internal pressure is the tangential stress, we may 
use the preceding equation with o = o7 and o = o7,. Making 
this substitution Eqs. (2)—(4) lead to 


(1 — dt = to 
/J0 


1 — 2(n — = (1 — if n (5) 


and hence 


In (1 — 2néq,t’) = —2nép to ifn = 


These equations relate the rupture lives at constant stress (¢)) and 
constant pressure (¢’) if the various constants (, p, and B) are 
known. In cases where the stresses are low, the time to rupture 
long, and the deformation to failure consequently small, the 
rupture life predicted by Eq. 5 may be considerably less than given 
by t’ = 1/2né7,. 

Unfortunately, there are few data against which we may com- 
pare the predicted maximum life of tubing under internal pressure. 
A single, rather limited comparison can be made from the test 
data of Fig. 1. This shows the tangential strain at the outside 
of a 2S aluminum tube which was loaded by internal pressure at 
250°C. (482°F.). The test conditions are such that we might 
expect the rupture life to be governed by the instability condi- 
tion. Diameter changes were measured during the test with a 
special high-temperature feeler gage. We can read from the 
figure that the rupture life ¢#’ = 72.3 hours and can estimate the 
initial (minimum) tangential strain rate as é7, = 8.0 XK 10~4 
hours!. Using the approximation ¢’ = 1/2né7,, we find n = 8.6. 
This is quite close to values of nm which may be estimated from 
data? in the literature for 2S aluminum at 500°F. and com- 
parable stresses and strain rates. 

Eq. 3 may be written as 

ep = —(1/2n) In (1 — 2nép,t) (6) 
Putting 1/2né7, = 72.3 hours and » = 8.6 in Eq. 6 leads to the 
broken curve shown in Fig. 1. Strictly speaking, the fact that 
this curve lies close to the measured values only proves that we can 
reconstruct the curve from its initial slope, é7,, and final abscissa, 
t’. However, further justification is given by the agreement of 
the predicted » with data in the literature. 
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° 


STRAIN 


TIME (HOURS) 


Solid curve shows tangential strain at outside diameter 
of a 2S aluminum tube, 1.600 in. outside diameter and 0.100 
in. wall thickness, under 290 psi internal pressure at 250°C. 
Broken curve shows strains calculated from Eq. (6). Thin solid 
line shows initial (minimum) creep rate. 
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A Simple Construction for the Determination 
of the Magnetohydrodynamic Wave Speed 
in a Compressible Conductor 


Nicholas Rott 
Professor of Aeronautical Engineering, Graduate School of 
Aeronautical Engineering, Cornell University, Ithaca, N.Y. 


November 19, 1958 


HE PROBLEM OF THE MAGNETOHYDRODYNAMIC plasma waves 
g ped solved by Herlofson! and van de Hulst,? and an excellent 
review of the subject was given by Friedrichs,* whose report, un- 
fortunately, belongs to the more or less unavailable ‘‘subliter- 
ature.”” The purpose of the present note is to point out a simple 
geometric construction, which gives not only the wave speeds as 
a function of the wave direction, but also the direction of the re- 
sultant particle motion (‘‘polarization’’) for the different modes. 
The usefulness of this construction, which, naturally, does not 
contain new information, is that it permits a quick survey of all 
possible plasma wave modes—a desirable result in view of the 
complications of the problem. 
Let a be the speed of sound in the perfect gas, in the absence 
of magnetic fields: 
a? = 7(p/p) (1) 


and let b be the Alfvén velocity in an incompressible conductor 
(of the same density p) in the presence of a uniform magnetic 
field H (using 1/K SQ units) : 

b? = pH?/p (2) 


Consider a plane wave in the compressible conductor whose 
normal forms an angle J with the direction of the magnetic field. 
The field components are H, = H cos 8 and H, = H sin # parallel 
and perpendicular to the wave normal; the corresponding com- 
ponents of the particle velocity behind the wave, in the plane 
formed by H and the wave normal, will be called Au and Av. 
[The motion in the third direction (z) will be considered later. | 
The symbol A will be used for the change of any quantity across 
the wave, and c is the wave speed which is to be determined. 
Continuity requires that 
pAu = cAp (3) 


To simplify the present analysis, a result obtained by previous 
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investigations will be quoted without proof: weak magneto- 
hydrodynamic waves are isentropic. This result, which can be 
used, as in gasdynamics, in lieu of the energy equation, states 
that 
Ap = a?Ap = p(a?/c) Au (4) 
The momentum equations are influenced by the presence of a 
body force, caused by the interaction of the magnetic field with a 
current sheet flowing ‘‘in’”’ the wave (in s-direction, perpendicular 
to the plane of 7 and the wave normal). The strength of the 
current sheet, in turn, is measured (according to Ampére’s law) by 
AH, so that the two momentum equations take the form 
Ap — pcAu + wH,AH, = 0 (5) 
pcAv — wH,AH, = 0 (6) 
Finally, the induction law (for infinite conductivity*) will be 
formulated, namely, that the magnetic flux across a particle-fixed 
closed line remains unchanged. Balancing the changes across a 
loop formed by lines parallel to the x- and s-directions, one finds 
cAH, = H,Au + H,Av (7) 
Elimination of Ap and AH, in Eqs. (5) and (6) by Eqs. (4) and 
(7) gives 
(a? — c? + by*?)Au + b2b,Av = 0 
—brb,Au + (c? — b,?)Av = 0 
where = bcos and b, = b sin J. Finally, 
Av/Au = b,b,/(c? — 6,2) = (c? — a? — b,?)/brb, (8) 
and ¢ fulfills the equation 
— c%(a? + b?) + a*b? cos? = 0 (9) 
The two solutions for c? are 
= (1/2)(a® + b?) (1/2)V at + — 272b? cos 29 (10) 
The angle y between the particle motion and the wave normal 
follows from Eq. (8). It is advantageous to construct tan2y 
as a product of the two expressions given in Eq. (8): 
tan? y = (Av/Au)? = (c? — a? — b,*)/(c? — 5,2) 
and to consider the angle 2y: 
cos 2y = (1 — tan*y)/(1 + tan’y) = 
(a? + by? — bz?) /(2c? — a? — 5?) 
which can be written by the use of Eq. (10) in the form 
cos = +(a? — b? cos 289) — (11) 
The well-known results expressed in Eqs. (10) and (11) lead 
to the following geometrical construction: let the direction of the 
magnetic field be along OH, and let ON lie in the direction of the 
wave normal, so that the angle HON = # (see Fig. 1). Consider 
B on ON and E on OH such that OB = BE = b?/2; this makes 
the angle NBE = 23. Let BA = a?/2 so that 


* This assumption can be dropped if the changes between two uniform 
states before and after a weak shock (or wavelet) are considered. 
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OA = (1/2)(a? +62) = (1/2)(c,? + 
according to Eq. (10); also, 


AE= (1 + bt — 2a*h? cos 20 

= — 

The circle around A with the radius AE cuts ON at the points 
and so that 


OCG, = a2, OCs = 


Moreover, EC; and EC, are the directions of particle motion, or 
the angle OQ. E = OGE = 

Proof: OAE can be shown to be 2y, by making ED perpen- 
dicular to ON, so that 


BD = BE eos 28 = (1/2)b? cos 20 
and AD = AB — BD = (1/2)(a? — b? cos 29) 


and finally, AD/AE = cos 2y, according to Eq. (11). 

The two directions of particle motion in the two modes are 
perpendicular to each other, as is already known. 

The construction (which by symmetry only has to be used in 
the range 0 < 3 < 7/2) has the drawback that the ‘‘velocity 
squares” are plotted, but it offers the advantage that it also 
gives the direction of the particle motion for the modes. It can 
be shown that if the roles of a? and 6? are reversed, the construc- 
tion remains valid if the roles of H and N are also reversed. 

Finally, a short remark can be made regarding a third possible 
mode, in which the velocity of the particles is directed in the z- 
direction, perpendicular to the two directions obtained so far. 
It can be shown that this mode remains uncoupled with the com- 
pressibility effects of the gas, so that its propagation speed equals 
the Alfvén velocity calculated with the field component parallel 
to the wave normal, viz., c; = b cos 3. It is found that c¢; is al- 
ways intermediate between c, and @, except for J) = O when c; 


coincides with cs. 
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Concise, Systematic Jacobian Thermodynamics 
for High-Temperature Real Gases 
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SYMBOLS 


4 = dimensionless Helmholtz function 
F = dimensionless Gibbs function 
G = dimensionless internal energy 
I = dimensionless enthalpy 
= Jacobian operator 
| = pressure 
R = the gas constant 
S/R = dimensionless entropy 
T = absolute temperature 
Z = P/pRT = compressibility factor 
= J(in p, In T) 


a 
b = J(Z,1np) = J(G, In T) 

c = J(Z,G) = (b? — mn)/a 

m = J(in T, Z) 

n = J(In p, G) 

bia = = Gp = —(0Z/dln T)p = —Zr 
c/a = Gp? + GrZp 

m/a = —(0Z/dl1n p)r = —Zp 

n/a = (0G/dln T)p = Gr 
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xand y = general variables 
In = natural logarithm symbol 
= density 


Mos pAY high-velocity flight phenomena are leading 
scientists into a large number of high-temperature prob- 
lems. One such problem is introduced by dissociation and ioni- 
zation of gases. By either decreasing the density or increasing 
the temperature of a gas, one can force a nearly ideal gas to change 
into a so-called ‘‘real’’ gas which takes into account the effects 
of dissociation, ete. Shaw,? by functional analysis (or Jaco- 
bians), has developed a unique orderly method for obtaining any 
thermodynamic partial derivative for a simple system — e.g., an 
ideal gas — in terms of reference partial derivatives. Presented 
here is a table of common thermodynamic quantities that are 
generalized in order to account for real-gas effects. 

A thermodynamic partial derivative can be represented by a 
combination of any six (or less) chosen thermodynamical partial 
derivatives plus certain state variables.* These six partials can 
be reduced to three by the three following steps: proper choice 
of independent variables (p, 7); Maxwell’s relation (G, = —Zr), 
and interrelation of the derivatives (ac — 6? + mn = 0). The 
terms in parentheses that follow each of these three steps are 
for the particular case which is to be presented here. That is, 
all of the thermodynamic quantities are taken to be functions of 
density and temperature. The dimensionless internal energy 
and the compressibility factor are used as the working dependent 
variables, thus establishing the three reference partial derivatives 
Gp, Gr, and Z,._ Although this choice of parameters is arbitrary, 
very accurate numerical values for internal energy and com- 
pressibility can be obtained by using spectroscopic data. Also, 
temperature and density have been found by experience to be 
convenient variables with which to work. 

The fundamental equation that relates a Jacobian table to a 
partial derivative is 
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The Jacobian formulas of importance for filling in the tables in 
terms of reference quantities are 


J(in P, x) = (1/Z)J(Z, x) + J(In T, x) + J(in p, x) (2) 
J(S/R, x) = J(G, x) + G J(In T, x) — Z J(in p,x) = (3) 


JU, x) = J(G, x) + J(Z, x) (4) 
J(A/RT, x) = Z J(\n p, x) — G J(In T, x) (5) 
J(F/RT, x) = J(A/RT, x) + J(Z, x) (6) 


where x is any thermodynamic variable. 

In order to obtain the value of J(x, y) from the tables, let the 
first variable (x) apply to the desired quantity in the left-hand 
column and let the second variable (y) represent the desired 
quantity in the row across the top of the tables. Using Table 1 
let us find the value of (OP /0p)r by means of reference quantities 
whose values are known. Making use of Eq. (1), we find that 


(OP/dp)r = (P/p)(O1n P/dln p)r = 
(P/p) J(in P, In T)/J(in p, In T) 
(P/p) [a — (m/Z)|/a = (P/p) [1 — (m/Za)| 
Or, by Table 2, we would find in a like manner that 
(OP/Op)r = (P/p) [1 + 

Table 1 is the more general table in that the choice of inde- 
pendent variables has not been made, whereas Table 2 makes p 
and TJ the independent variables —each value in Table 1 is divided 
by a = J(In p, In 7). Thus other tables similar to Table 2 can 
be made by dividing Table 1 through by other such quantities. 
Table 2 enables one to find the numerical value of any partial de- 
rivative (composed of the quantities heading the tables) in terms 
of just Gp, Z,, and G7, along with G and Z which are all known 
and are all functions of pand 7. 

The tables are a useful, convenient short-cut procedure for the 
handling of almost any thermodynamic partial derivative situ- 
ation. Other related manipulations and second-order partial 


(O0x/dy)2 = J(x, s)/J(y, 2) (1) derivatives! can be handled in much the same manner as pre- 
TABLE 1 
J(x, ¥) J(x,y) = —J(y, x) J(x,x) = 0 Example: J{[(A/RT), G] = Gb + Zn 
| ym | In P In p in T S/R Zz I A/RT F/RI 
| 
b — Ga — Za G 
b 
In p “see 0 a Gat+n —b n—b —Ga —b — Ga 
m 9 
In 7 —a 0 Za —b m m—b —Za m— Za 
—Ga 
| Gm Gm — Gb 
Z | 
G |} —n b Gb + Zn —Gb — Zn —Gb — Zm 
| 
| —m 0 Gu + 2b Gm + 2b 
- + 2b Z | —Gm + Gb Gm + Zb Gm + Zb 
I | ° —Gb Zn | —Gb — Zn - 
| —b + Ga + Za | : 
A/RT | Gm | Ga Za Gb + Zn Gb + Zn —Gm — 2b =a 0 —Gm — 2b 
} + Gb + Zn 
—Gm + G 2b 
Gm | b6+G — Gm + Zb 0 
+Za- — +Zn — Zh+c inte + Gbh+Zn+ 
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sented here. For further details on Jacobians, one should con- 
sult Shaw.? 
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On the Heating Effect in a Resonance Tube 


|. M. Hall and C. J. Berry* 

Aerodynamics Division, National Physical Laboratory, Teddington, 
England 

November 17, 1958 


- A RECENT NOTE! Sibulkin and Vrebalovich have drawn atten- 
tion to Sprenger’s work on heating in resonance tubes.? It 
is mentioned that the dissipative mechanism inside the tube 
might be shock waves, and recent work at the National Physical 
Laboratory has shown that this is the case. 

A 2-in. long, 3/16-in. square resonance tube, closed at one end, 
was excited by a supersonic jet issuing from an orifice of the same 
cross section as the tube. The pressure ratio across the orifice 
was about 4. When the mouth of the resonance tube was in one 


* The authors acknowledge the assistance of R. J. North who designed 
the Cranz-Shardin system. 


(a) 20 microseconds 


(b) 140 microseconds 


(c) 260 microseconds 
Ba 


(d) 425 microseconds 


Fic. 1. Shadow photographs of the resonance tube showing 
four stages in the cycle. The times shown are measured from 
an arbitrary zero. In (a) and (b) the shock wave, indicated by 
a white arrow, is moving toward the bottom of the tube while 
in (c) and (d) it is moving away from the bottom. 
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of the compression regions of the jet, the air column in the tube 
resonated, producing an intense note with a wavelength of ap- 
proximately four times the length of the tube. The most in- 
tense sound occurred when the mouth was in the first compres- 
sion region, and in this case the temperature at the closed end of 
the tube was about 140°C., the total temperature of the jet being 
about 20°C. 

Short-duration (about 1 microsec.) shadow photographs were 
taken using a Cranz-Schardin system which takes eight photo- 
graphs at intervals of the order of 50 microsec. Four photographs 
from such a series are shown in Fig. 1. The shock wave is clearly 
seen in its passage down the tube but is more difficult to distinguish 
on its return [(c) and (d) in Fig. 1]. A correlation between the 
flow inside and outside the tube is clear from these photographs. 
When the shock wave is moving into the tube, the external 
pattern is similar to the case of a jet impinging on a flat surface 
with spillage sideways as a wall-jet. In the second pair of photo- 
graphs, where the internal shock wave is moving out of the tube, 
the jet appears to penetrate into the tube with little or no spillage. 

The precise mechanism of the heating effect is not yet com- 
pletely understood although it is clearly related to the nonisen- 
tropic flow inside the tube. Further investigations are in progress 
which, it is hoped, will enable a complete explanation of the proc- 
ess to be obtained. 
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Linear Deformations of Conical Shells 


M. D. Mc Ilroy 
Bell Telephone Laboratories, Inc., Murray Hill, NJ. 
November 10, 1958 


eee NOTE DESCRIBES a complete set of solutions to the linear- 
ized problem of elastic deformations of thin shallow conical 
shells.* The relevant differential equations are' 


CV'F + tan ad*w/dr? = Ol 


DV‘w — r~' tan ad*F/dr? = 0 f (1) 


where F is a stress function, w is transverse displacement, a is 
elevation of cone generators above a base plane normal to the 
shell axis, and C and D are material constants. 

The system in Eq. (1) is transformed into a set of independent 
fourth-order ordinary differential equations of generalized hyper- 
geometric type,’ 

26(6 —1)} (2) 
6 = sd/dz, n=0,1,2.. 


by writing w — i(CD)~!/2F = W,(s) cos n0 (3) 


and performing straightforward algebraic transformations. The 
independent variable z is 


s=ir/a, a = cot aVD (4) 


Two independent “cone function” solutions of Eq. (2), found 
by the Frobenius method, and denoted Co,(z, ) and Coo(z, nm), 
are elementary functions. Two other solutions, found through 
Euler transform integral representations, and denoted Co;(z, ) 
and Co,(z, n), are transcendental functions expressible in terms 
of Bessel functions. In terms of the cone functions, solutions 
to Eq. (1) may be expressed by 

* For details of this work see a doctoral thesis accepted by the Massa- 
chusetts Institute of Technology, September, 1958, entitled On the solution 
of the differential equations of conical shells, and supervised by E. Reissner. 
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4 
= > [A; Re Co,(ir/a, n) + 
j=l 
4 B; Im Co,(ir/a, n)| cos (5) 
F//CD = > [—A; Im Co,(ir/a, n) + 
j=l 


B; Re Co,(ir/a, cos j 


where the A, and B, are arbitrary real constants. 

As solutions to problems concerning conical shells of sufficiently 
large extent that r/a > 1, the functions Co; and Co, embody 
boundary-layer (exponential edge effect) deformations. The 
functions Co; and Coz, however, have different interpretations 
accordingly as m = 0, m = 1, orn > 2. Associated with the 
constants of Eq. (5) are the following types of solution: 


A» interior state 
By, dislocation 
Az» dislocation 
By interior state 


n =0 A, rigid motion 

B, nugatory 

A, nugatory 

B, rigid motion 

n>2 Ao, Be interior states 


Here an interior state, as distinguished from a boundary-layer 
state, is any deformation of membrane and/or inextensional 
bending type. A nugatory solution is one for which all stresses 
and displacements vanish. On the assumption of rotational com- 
pleteness of the shell, dislocations (in the sense of Volterra) as 
well as nugatory solutions must be rejected. In their stead may 
be appended an arbitrary torsional deformation, axial rotation, 
lateral translation, and the singular solution 


w = A» Re Co.(ir/a, 1) cos 6 t 
(6) 

F/ CD = Im Co,(ir/a, 1) cos 6 + Ax(r/a)6 sin 
Asr— © ali meaningful solutions defined above approach in an 
appropriate sense approximate solutions given for the conical shell 
by Love,’ while as r ~ 0, the solutions approach those for a flat 
plate. 

The solution to any properly set boundary-value problem in- 
volving shallow (tan a < 1) conical shells can be expressed in 
terms of solutions outlined above. Also solutions to Seide’s 
approximate differential equations for the deep conical shell,‘ 
which are identical to the system in Eq. (1), are comprised here; 
except that the m of the shallow-shell theory must be replaced by 
n/cos a, Eq. (2) follows from Seide’s equations. The cone func- 
tions are adequate to any case for which /cos a is an integer, but 


ee a satisfactory extension of the cone functions to nonintegral order 
e is required to handle Seide’s equations in general. 
tee Formulas for the cone functions are: 
Co,(z, 0) = 1, Co(z, 1) = (1/2) g, 
"—2 — 2 — — 


Co(z,0) = —logs +s + ¥(2) + W(1), (8) 


Co.(z,1) = 327! — 2, 


+ 2) _, 


— 1)r(2n+ 
(n+ 1)r(n) 


"=? 1(2n — 2 — k)V(2n — 

T(n — — 1 — + 
[W(n — k) + — 1 —k) — — 2 —k) 
W(2n — k) — + 2) — W1 + + 2/(n? 1), 


n=2.3... 


Co,(z, 0) = f 1231/7) ds — 1, 
‘ (9) 
Co(z, m) = d"/ds" | 


Zz 
Co,(z, 0) = -2f K2(22!/?)dz, (10) 
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Laminar Jet Mixing of Electrically Conducting 
Fluid in a Transverse Magnetic Field* 


S. |. Pai 

Research Professor, Institute for Fluid Dynamics and 
Applied Mathematics, University of Maryland, 
College Park, Md. 


November 24, 1958 


ly THIS NOTE we discuss the effect on laminar jet mixing of an 
electrically conducting, incompressible, and viscous fluid 
under an applied transverse magnetic field. We assume a jet of 
this fluid issuing from a small two-dimensional slit and mixing 
with the surrounding medium of the same fluid. There is an ex- 
ternally applied transverse magnetic field Hy(x) = H,x~" where x 
is the distance along the axis of the jet and ” is an arbitrary con- 
stant. There is no external pressure gradient in the present 
problem. We shall consider only the case of small electrical 
conductivity o or small magnetic Reynolds Number so that the 
induced magnetic field may be neglected. Under this condition 
and the usual boundary-layer approximations, the fundamental 
equations of our problem are! 


u(Ou/Ox) + v(Ou/dOy) = —(Q/9)x-*" u 4+ v(07u/dy?) (1) 


(Ou/Ox) + (dv/dy) = 0 (2) 

where the external magnetic field is taken as 
H,(x) = (3) 
and Q = ue? Hi? (4) 


H)(x) is in the y-direction which is perpendicular 
to the jet axis. p is the density of the fluid and yu, is the mag- 
netic permeability. The MKS system is used. 

The boundary conditions of our problem are 


H, isa constant. 


y = 0: 


y= 


v = 0, (Ou/doy) = 0 
u=0 (5) 


We introduce the stream function y such that 


u = Oy/dy, v = (6) 
Eq. (2) is automatically satisfied. We replace the x — y co- 
ordinates by the following variables 

x, 0 = (y/x?/4) (7) 


In order to solve Eq. (1), we write the stream function y in the 
following form: 


y= 2 3 n= Qe" fi(n) 
(Qe")? +...) (8) 


* This work was supported by the Air Force Office of Scientific Research 
under Contract AF 49(638) 401. 
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where m = (4/3) — 2n. 


Substituting Eqs. (6) to (8) into Eq. (1) and coilecting terms 
of the same power of Q¢", we have 


fo’’’ + (fo fo’) = 0 (9) 


fil’ — (Bm — 2)fo'fi’ = fo’, ete. (10) 
With the boundary conditions 
=... =O;fo” =fi’ =...=0 (11) 
n= =...=0 


where prime refers to differentiation with respect to 7. 

Eq. (9) is the well-known equation of laminar jet mixing of 
ordinary fluid dynamics.? The equation for fi, fo, ete., are all 
linear differential equations, with variable coefficients, and they 
may be easily solved by numerical integration. 

For first approximation, if we take f = f, + Q¢”fi, the incre- 
ment of the axial velocity component u over the corresponding 
value in the ordinary fluid dynamics increases with Q. 


SIMILAR SOLUTION 

If nm = 2/3, m = 0. Such a solution may be called a similar 
solution because the shapes of velocity distribution u at various 
x-stations are similar. The solution (8) to (11) for this problem 
holds true for small value of Q because we use the method of series 
expansion. 

Some general conclusions for this similar solution may be drawn 
from the following considerations. In the similar solution 
f(¢, n) = f(m) which is satisfied by the equationt 


(ff’)’ Of’ (12) 
if we expand the function f(7) in power series of the parameter 


Q, we have the same results as given in Eqs. (8)-(11). 
Now we introduce the following transformation 


f=VOF(V0n) = (13) 
Eq. (12) becomes 

F'" + (FF’)' = F’ (14) 
where primes refer to differentiation with respect to s. Now it is 
interesting to notice that F is independent of Q explicitly—i.e., 
the same function F(z) holds for any value of Q within the ap- 

proximations of our analysis. 

The x-component of velocity in the jet is 

u = (1/3x!/3) (df/dn) = (Q/3x"/*) (dF/dz) (15) 


Heie uw increases with Q which checks with our result by the 
method of series expansion. 
The spread of the jet mixing region is proportional to 


b = (16) 
Hence the spread of the jet mixing region decreases as (1/V 0) 
decreases. 
The important parameter Q may be written as 
Q = (9U/L)RuRe (17) 


where Ro = we oUL = magnetic Reynolds Number (18) 


Ru = p-H,?/pU? = modified magnetic pressure number (19) 


If Hy is independent of x, Ry reduces to the conventional mag- 
netic pressure number. Here U and L are the characteristic 
velocity and length of the problem respectively. The parameter 
Q increases with the product of Ro and Ry. 


REFERENCES 
Rossow, V. J., On the Flow of Electrically Conducting Fluids Over a Flat 
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+ Prof. J. M. Burgers pointed out to the author that there is no analytic 
solution of Eq. (12) throughout the whole domain to satisfy the boundary 
condition (5). This point will be further discussed elsewhere. 
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le THE ANALYSIS of a beam column, one is sometimes con- 

fronted with the problem of effective end restraint. Although 
the rotational rigidity of the supports can often be determined 
with comparative ease and accuracy, the matter is usually re- 
solved by using an arbitrary ‘‘fixity factor’’ to reduce the effective 
column length, with no attempt made to also decrease the maxi- 
mum moment as a logical consequence. In either case, consider- 
able labor is required to compute the true moments, and the 
procedure is not suitable for routine analysis. 

The following discussion shows that it is possible to represent 
the effects of end restraint and lateral load as functions of a single 
nondimensional parameter. Thus, regardless of whether arbi- 
trary fixity factors or actual rotational rigidities are used, it be- 
comes feasible to represent many cases of practical interest by a 
series of plots which reduce the computational work consider- 
ably. The theory upon which this note is based can be found 
in a number of texts—e.g., Chapter 1 of reference 1. 

We consider the beam column shown in Fig. 1 with an arbi- 
trary lateral load, W, and end moments 1/4 and J/,% due to the 
elastic nature of the supports. Since the final rotation of the end 
tangents is a linear combination of the separate effects, we may 


write 


= ki Ma + + fa(W) 


Op = + + fa W){ 


where the &’s and f’s may be thought of as generalized stiffness 
coefficients. 

If, noting reference 1, we now refer to 1/4 and Mz, as the 
moments produced by the rotations against the elastic restraints 
at A and B, we have 

064 = —(Ma/a) O08 = —(Mz/8) 
where a and £6 are coefficients describing the degree of fixity at 
A and B—i.e., a = moment required to produce unit rotation at 
A,ete. If we let 
y = fa(W)/fa(W) 


Eqs. (1) may be rewritten as 


[ki + (1/a@)| Ma + = —fa(W) (2) 
+ [ki + (1/8)| Me = —vfa(W)| 
If we further define 
ky + (1/a) ky 
Then: 
vke — [ki + (1/8)] 
Ma = — }) 
k [ki + (1/a@)] 
Mp = saw) \\ 


It is evident, therefore, that, except for the nondimensional 
factor, y, the bracketed expressions in Eqs. (4) are independent of 
the lateral load W. If the load is symmetrical with respect to 
the beam center line, y = 1.0. Otherwise, the dependency is of 
the simple form shown above and readily evaluated. 

It is easily shown that 


ky 


(A cot Al/P) — (1/PI) 

(1/Pl) — (A/P sin 

= VP/EI 

Substituting the above values of k; and k2 into Eqs. (4), we 


arrive, after some algebraic manipulation, at the following ex- 
pressions for M4 and Mz: 


where 


0), 
| 
(5) 
(6) 
a 
| 
(7) 


256 JOURNAL OF THE AERO/SPACE SCIENCES 


Ma Me 


(aB/PL [f(D] + +a 
Vi(a + [f(D] + (2a8/P2) + 1f 


(5a) 


= fa(W) 


+ (aB/PL [fa] + | 
+ B)/PI [AOD] + (208/ P22) + 
(5b) 


fa(W) - 


where fi(Al) = Al cot NM — 1 
fo(Al) = (Al — sin /sin AL 
= Al[tan (Al/2) + (Al/2)] 


Except for a possible additive or multiplicative constant, it is ob- 
vious that f( 1”) may also be expressed as a function of X/. 
Hence a solution for the end moments now becomes feasible 
by representing the various functional relations in a series of plots 
where, for example, one plot of f,(A/), where » = 1, 2, 3, suffices 
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for any column loading. The principle of superposition then 
permits a limited number of plots of f(W) and 7(A/) to account 
for almost any conceivable lateral load. It will be noted that for 
the linear combination represented by Eqs. (1), f( W) becomes in- 
finite as \/ approaches 7. Therefore, while the theoretical maxi- 
mum value of X/ is 27, the plots of the above functions can be 
terminated at Ni = rz. 

The restraint coefficients a and 8 may be chosen so as to corre- 
spond to some reduced column length, /’, such that c = (//l’)?, 
where c is the so-called fixity factor and / the actual column length. 
The relation between c and a and 8 may be established from the 
condition that at the critical load the determinant given by Eq. 
(3) must vanish. <A family of curves may then be drawn giving 
P.,/P, as functions of the nondimensional parameters al/EI 
and 8//EI. Here P., is the critical column load for the effective 
length /’ and P, the Euler load for a pin-ended column of length 
l. Since the definition c = P,,/P, also holds, the correspondence 
between c and a@ and @ is thus established. In closing, it should 
also be noted that the above family of curves may be used for 
determining the critical column load for any arbitrary degree of 
end restraint represented by a and 8. For stresses above the 
proportional limit, however, the determination only holds for 
effective column lengths, with the critical load then being found 
from any suitable column curve of the material. 
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